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Abstract 

High  speed  photographic  and  t Ime-resoi ved  spectrographic  observa¬ 
tions  have  been  made  of  the  head-on  collision  of  two  hypervelocity  copper 
jets  formed  from  metal -i I ned  explosive  charges.  Streak  camera  records 
show  that  the  jets  are  moving  at  about  8  mm/^sec,  or  equivalently,  with 
a  kinetic  energy  of  21  ev/atom.  Framing  camera  observations  of  the 
collision  in  air  and  in  a  low  pressure  helium  atmosphere  indicate  that 
a  region  of  intense  luminosity  is  created  about  the  point  of  impact. 

Time- integrates  spectrograph  records  show  that  the  composition  of  the 
collision  light  over  the  visible  region  Is  that  of  a  continuum  with 
superposed  copper  spectral  lines,  A  time-resolved  spectral  record  in  the 
4400-5300A  region  of  a  collision  in  helium  shows  that  an  intense  contin¬ 
uum  first  appears  for  about  5  lisec.  The  continuum  intensity  then 
decreases  and  a  copper  emission  line  spectrum  emerges,  consisting  of 
both  neutral  and  fl rst- ionized  copper  lines.  The  ionized  copper  lines 
originate  from  energy  levels  that  are  about  25  ev  above  the  neutral  atom 
ground  state. 

Several  neutral  copper  lines  are  shown  to  possess  combined  instru¬ 
mentation  and  resonance-broadened  profiles.  Estimates  of  the  neutral 
copper  atom  density  in  the  jet  collision  plasma  are  derived  from  line 
halfwldth  measurements  resulting  from  theoretical  profile  fits  to  the 

observed  line  shapes.  Neutral  atom  densities  are  of  the  order  of 

19  -3 

4  x  10  cm  An  intermolecular  electric  field  created  by  copper  ions 
and  electrons  in  the  plasma  result  in  a  Stark  shift  of  the  observed 
Cul  4531  line.  The  copper  Ion  densities  determined  from  measured 
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wavelength  shifts  of  this  line  are  about  4  x  10  cm 

Oscillator  strengths  for  two  neutral  copper  transitions  observed  in 
the  time-resolved  spectrum  are  also  determined.  A  comparison  of  the 
resonance-broadened  halfwidths  for  the  Cut  **525  and  4643  lines  with  that 
of  the  Cut  4651  line  whose  oscillator  strength  is  known  results  In  an 
estimate  of  the  f-numbers  for  these  transitions. 

Conclusions  derived  from  existing  theoretical  models  of  hyper¬ 
velocity  impact  are  found  to  be  in  good  agreement  with  the  observations 
and  detailed  measurements  obtained  from  these  experiments.  Applications 
of  this  work  to  othjr  areas  of  interest  are  Indicated. 


Benjamin  J.  Pernick,  Author 
Stephen  J.  Lukas ik.  Thesis  Advisor 
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Part  l  Experimental  Techniques  and  Gross  Observations 


The  primary  objective  of  this  work  Is  to  Investigate  the  possi¬ 
bility  of  plasma  production  by  the  direct  collision  of  two  high  velocity 
high  density  metallic  jets.  This  appears  feasible  on  a  simple  energetic 
basis  since  the  kinetic  energy  per  atom  in  a  jet  is  well  in  excess  of  the 
energy  required  to  ionize  the  atom.  The  veiocity  of  the  jets  used  here 
corresponds  to  a  kinetic  energy  of  21  ev/atom  whereas  the  energy  required 
to  vaporize  and  ionize  an  atom  from  the  jet  material  is  of  the  order  of 
10  ev .  Conversion  to  thermal  energy  of  this  directed  kinetic  energy  upon 
collision  was  expected  to  result  in  excitation  and  ionization  of  the  jet 
material.  Primary  reliance  has  been  placed  on  the  use  of  the  time- 
resolved  spectrum  of  the  impact  light  as  a  diagnostic  aid  to  determine 
the  nature  of  the  collision  process. 


I . I  tntroduct ion 

Recent  experiments  concerned  with  the  impact  of  explosively  driven, 
high  speed  metallic  projectiles  and  targets  have  indicated  that  charged 
particles  along  with  vaporized  neutral  atoms  of  the  projectile  and  target 
materials  are  produced  upon  collision.  Qualitative  observations  of  the 
projectile-target  impact  by  Clark,  et  al  (1959)  and  Cook  and  Keyes  (1959) 
have  indicated  that  an  intense  luminosity  is  generated  upon  collision. 
Experiments  by  Allen,  et  a  I  (1959)  and  Grow,  et  al  (I960)  have  shown  that 
spectral  line  radiation  characteristic  of  neutral  atom*  of  the  projectile 
and  target  materials  was  generated  by  the  impact.  Keyes,  et  al  (I960) 
and  Fr  i  ichtenicht  and  Slattery,  (1963)  have  indicated  that  at  high 
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projectile  velocities  Ionized  atoms  of  the  projectile  and  target  materials 
are  present  in  the  Impact-generated  vapor,  tt  has  been  suggested  by 
Koslov  (19>9)  and  others,  (Allen,  et  al,  1959)  that  collision  between 
solid  materials  moving  at  sufficiently  high  relative  velocities  could 
result  in  the  production  of  plasmas  with  appreciably  high  ion  concentra¬ 
tions.  Development  of  this  technique  for  the  production  of  high  density 
plasmas  and  the  study  of  Its  nature  is  of  considerable  interest  and  value 
since  this  method  represents  a  different  procedure  in  the  present  day 
technology  of  plasma  production,  (Bishop,  1958;  Fischer  and  Mansur,  1958; 
lochte-Holtgreven,  1958) . 

Two  potential  applications  have  evolved  from  this  work  over  and 
above  the  stated  plasma  physics  interest.  It  will  be  shown  that  these 
experimental  methods  can  be  employed  in  the  dete  mination  of  atomic 
oscillator  strengths  or  f-numbers  for  the  material  constituents  of  the 
collision  plasma.  Such  results  would  add  to  the  existing  knowledge  of 
oscillator  strengths  presently  determined  with  different  experimental 
techniques,  (Allen  and  Asaad,  1957;  Corliss  and  Bozman,  1962).  Further¬ 
more  these  methods  could  be  used  in  principle  for  the  determination  of 
relative  oscillator  strengths  of  metallic  ions  for  which  information  is 
relatively  sparce  (Dickerman  1961;  Allen,  1963). 

Lukas  Ik,  et  a l  (1964),  have  shown  that  hydrodynamic  pressures  of 
the  order  of  tens  of  megabars  exist  at  the  moment  of  impact  for  the  jet 
velocities  used  in  this  work.  These  values  are  arrived  at  from  an 
assumed  equation  of  state  for  the  Jet  material  and  using  a  hydrodynamic 
mooel  of  the  collision  process.  The  subsequent  motion  of  the  jet  colli¬ 
sion  products  was  shown  to  be  dependent  upon  the  magnitude  of  this 


initial  impact  pressure.  Conversely,  experimental  observations  of  the 
motion  of  the  collision  products  can  from  this  hydrodynamic  modei  result 
in  a  measure  of  the  initial  impact  pressure.  This  procedure  Mould  be  of 
Importance  for  mate-iais  whose  equation  of  state  is  not  well  established 
or  for  which  no  such  measurements  have  been  made  at  such  high  pressures, 
in  addition,  with  the  advent  of  hypervelocity  projectiles  at  speeds  higher 
than  that  used  in  this  study  (Jameson,  1963)  it  appears  feasibly  that 
experimental  measurements  could  be  extended  into  the  hitherto  unattain¬ 
able  megabar  pressure  ranges. 

The  primary  aim  of  this  effort  is  to  quantitatively  determine  the 
degree  of  ionization  in  a  plasma  generated  by  the  impact  of  two  bodies 
moving  at  high  relative  speeds.  This  is  achieved  from  a  detailed  spec¬ 
tral  analysis  of  the  optical  radiation  from  the  collision-produced  plasma. 
It  will  be  shown  that  the  predominant  physical  mechanisms  that  affect  the 
observed  spectral  lines  are  due  to  resonance  interactions  and  Intermole- 
cular  ionic  fields.  Estimates  for  the  neutral  and  ion  particle  density 
of  the  collision  plasma  are  attained  from  a  comparison  of  predicted  and 
observed  spectrum  iine  profiles.  Collision  plasma  temperature  estimates 
ere  derived  from  the  results  of  the  analysis. 

1.2  Explosive  Shaped  Charge  Jets 

Certain  criteria  of  a  general  nature  must  be  considered  in  the 
decision  of  what  type  of  accelerator-projectile  system  should  be  used  In 
this  study.  First,  one  must  provide  for  an  efficient  conversion  of 
kinetic  to  internal  energy  of  the  collision  products  In  order  to  produce 
a  high  density  plasma.  If  the  projectile  material  density  is  that  of  a 


gas  then  the  ••colliding**  objects  will  in  effect  penetrate  each  other  with 
Httie  interaction  and  remain  relatively  uncoupled.  On  the  other  Hand  a 
strong, localized  Interaction  and  presumably  a  high  degree  of  thermal iza- 
tlon  of  the  Initial  directed  energy  Is  expected  from  the  collision  of  two 
solid  projectiles. 

in  addition,  the  projectile  velocity  must  be  sufficiently  high  so 
that  the  kinetic  energy  per  atom  is  at  least  greater  than  the  energy 
necessary  to  vaporize  and  ionize  the  atom.  Assuming  the  heat  of  vapor¬ 
ization  or  heat  of  sublimation  as  representative  of  the  necessary  energy 
to  remove  an  atom  from  a  solid  projectile  under  impact  conditions,  then 
the  minimum  requisite  kinetic  energy  sufficient  to  vaporize  and  ionize 
an  atom  Is  typically  of  the  order  of  10  ev  for  metallic  projectiles. 

Several  different  methods  are  available  for  producing  high  velocity 
projectiles.  Gas  guns  (Boyd,  et  ai,  1959;  Crews,  1959),  electrostatic 
and  electromagnetic  accelerators  (Bergstraih,  et  al,  1959),  exploding 
wires  (Chace  and  Moore,  1962),  explosively  driven  pellets  (Summers  and 
Charters,  1959;  Kineke,  1959),  explosive  shaped  charges,  and  plasma  guns 
have  been  used  to  accelerate  materials  to  high  speeds.  Although  very  high 
velocity  plasmoids  can  be  obtained  from  plasma  guns,  it  appeared  that  a 
significant  degree  of  thermal izat ion  of  the  initial  kinetic  energy  would 
be  less  likely  to  be  achieved  due  to  the  tow  material  density.  Of  the 
remaining  possibilities,  however,  only  explosive  techniques  have  resulted 
in  the  production  of  high  speed  projectiles  with  a  substantia)  density. 
This  reflects  the  fact  that  the  available  energy  density  of  a  chemical 
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explosive  is  extremely  high,  of  the  order  of  10  joules/cm  of  high  explo¬ 
sive  (Cook,  1958)* 
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Metallic  jets  formed  with  explosive  shaped  charges  were  used  as 
projectiles  In  this  work  rather  than  explosively  driven  pellets.  This 
was  desirable  since  the  velocities  attainable  with  shaped  charge  Jets 
are  much  larger  than  that  for  pellets,  typically  8  mm/^sec  as  compared 
to  2-4  mn/|Asee.  An  additional  factor  of  two  Increase  in  relative  impact 
velocity  was  achieved  by  firing  two  jets  directly  at  each  other  causing  a 
head-on  collision,  rather  than  one  jet  into  a  stationary  target.  Thus, 
the  kinetic  energy  per  atom  of  a  shaped  charge  Jet  moving  with  a  velocity 
of  8  irm/tisec  is  0.33A  ev/atom,  where  A  Is  the  atomic  weight  of  the  Jet 
material . 

The  formation  of  a  shaped  charge  jet  is  described  briefly  as 
follows.  Consider  a  block  of  high  explosive  with  a  metal  liner  embedded 
into  one  end  as  shown  In  Fig.  1(a).  The  explosive  is  initiated  at  the 
face  opposite  the  liner  and  the  detonation  front  moves  as  indicated. 

When  the  detonation  reaches  the  liner,  the  liner  is  deformed  and  collapses 
inward  as  illustrated  in  Fig.  1(b).  As  the  liner  materiel  arrives  at  the 
axial  region  a  high  speed  jet  is  developed  and  moves  in  the  indicated 
direction.  A  comprehensive  review  of  the  theory  of  shaped  charge  jet 
formation  and  experimental  verification  have  been  presented  by  Birkoff, 
et  al  (1948) ;  Pugh,  et  a)  (1952);  Eichelberger  and  Pugh  (1952);  Walsh, 
et  a  I  (1953);  Eichelberger  (1955);  end  Cook  (1958). 

Cylindrical  explosive  shaped  charges  with  a  conical  liner  insert  of 
copper  were  used  In  this  experimental  work.  The  choice  of  copper  as  a 
liner  material  was  based  on  its  ability  to  form  reproducible  and  predict¬ 
able  jets  and  on  the  availability  of  standard  explosive  charges.  Equally 
important  aside  from  the  charge  technology  is  the  fact  that  spectral  line 
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and  energy  level  information;  ionization  potential,  and  equation  of  state 
data  are  known  for  copper.  Specifications  of  the  explosive  design  are 
given  In  Table  !.  These  shaped  charges  produce  a  jet  whose  tip  travels 
at  about  8  mm/jAsec  or  equivalently  with  a  kinetic  energy  of  about  21  ev 
per  copper  atom.  This  latter  figure  is  of  interest  since  it  is  repre- 
sentatlve  of  the  explosively  derived  energy  directly  available  to  the 
impact-generated  plasma. 


1.3  Experimental  Configurations 

Two  types  of  explosive  charge  mountings  were  employed  in  these  jet 
collision  experiments.  For  rirings  in  air,  both  shaped  charges  were 
simply  mounted  onto  a  wooden  support  frame  and  secured  in  position. 
Alignment  of  the  two  opposing  shaped  charges  was  made  by  inserting  between 
the  charges  a  collapsible  rod  whose  ends  were  machined  to  the  liner  cone 
angle. 

^5f 

For  one  firing  in  a  low  pressure  atmosphere,  (firing  no.  109  ),  the 
assembly  shown  in  Fig.  2  was  employed.  The  charges  were  mounted  behind 
1/2  In.  thick  aluminum  plates  which  acted  as  baffles  to  delay  the  explo¬ 
sive  reaction  products  from  obscuring  the  jet  collision  region  between 
the  plates.  A  cone  base-to-ba$e  separation  of  4  in.,  or  two  charge  diam¬ 
eters  was  used.  It  was  feit  that  any  significant  deviation  from  axial  jet 
motion  would  be  small  at  these  short  standoff  distances.  The  jet  from 
each  shaped  charge  passed  freely  through  a  hole  in  the  center  of  each 


A  tabulated  summary  of  experimental  firings  appears  in  Appendix  A. 


baffle  plate  Into  the  region  between  the  plates.  The  other  rods  and 
plates  are  supporting  structures  to  hold  the  explosive  charges  securely 
In  place. 

A  3  ft  section  of  steel  pipe,  I  ft  Inside  diameter,  1/2  In,  wall 
thickness  served  as  a  vacuum  chamber  for  this  configuration,  as  shown  In 
Fig.  3.  Three  rectangular  observation  ports  were  included  in  the  pipe 
section  to  facilitate  optical  observations  of  the  jet  collision.  In  sub¬ 
sequent  firings  expendable  vacuum  chambers  were  used.  These  consisted  of 
a  cylindrical  section  of  clear  iucite,  2  ft  in  length,  1  ft  outside  diam¬ 
eter,  and  1/4  in.  wail  thickness.  The  explosive  support  mounts  were 
similar  to  that  described  above.  All  vacuum  systems  were  operated  at 
pressures  below  1  mm  Hg. 


1.4  Instrumentation 

All  of  the  preparations  involving  explosive  charges  and  firings  were 
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performed  at  a  Plcatinny  Arsenal  explosive  facility.  A  plan-view  of  the 
firing  chamber  is  illustrated  ‘n  Fig,  4.  Optical  observation  of  the  explo¬ 
sive  system  during  a  firing  is  possible  by  means  of  three  iucite  covered 
openings  or  ports  in  the  firing  chamber  wall.  The  lines  of  sight  through 
each  port  are  coincident  at  the  position  of  the  explosive  system,  thereby 
allowing  convenient,  simultaneous  observation  of  an  explosive  event  by 
several  optical  systems.  A  firing  is  viewed  with  a  Beckman  and  Whitley 
Hodel  No,  189  Framing  Camera  located  at  one  of  the  observation  ports  in 

*A  detailed  description  of  the  faciiity  is  given  in  the  "Proceedings  of 
the  14th  Meeting  of  Picatinny  Arsenal  Scientific  Advisory  Council"  at 
Plcatinny  Arsenal  25,  2b  April  1957. 
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In  the  camera  room.  This  camera  provides  a  film  record  that  consists  of 
25  frames  on  35  ran  film.  A  time  Interval  between  frames  of  1.4  y,sec  was 
found  appropriate  for  these  experiments.  The  Individual  frame  exposure 
time  at  this  framing  rate  Is  about  0.5  usee. 

A  Hilger  Medium  Quartz  Spectrograph  mounted  at  an  adjacent  observa¬ 
tion  port  provided  a  time- Integrated  spectrum  of  a  shaped  charge  jet 
collision.  Although  a  t ime- integrated  spectral  record  is  of  limited  use 
for  a  detailed  inierpretat ion  of  rapidly  changing  luminous  events,  ft  Is 
nevertheless  of  some  value  in  its  ability  to  cover  a  wide  wavelength 
range.  Its  chief  function  therefore  is  to  aid  in  spectral  line  identifi¬ 
cation  of  jet  collision  luminosity  over  a  broad  wavelength  band.  The 
spectra  of  other  light  sources  were  recorded  on  the  same  plate  as  a  jet 
collision  spectrum  for  wavelength  calibration. 

A  Beckman  and  Whitley  Model  No.  194  Continuous  Writing  Streak 
Camera  was  used  in  earlier  experiments  principally  to  measure  the  shaped 
charge  jet  velocity  prior  to  impact.  In  subsequent  firings  the  entrance 
optics  to  the  streak  camera  were  modified  to  include  a  direct-vision 
prism,  which  led  to  a  time-resolved  spectral  record  of  a  jet  collision. 

A  detailed  discussion  of  this  optical  arrangement  5s  given  in  Appendix  B. 

The  high  emulsion  speed  Kodak  Tri-X  or  Pan-X  35  ran  film  was  used 
for  ail  streak  and  framing  camera  records;  Tri-X  alone  was  used  for  all 
time-resolved  spectral  records.  The  Hilger  spectrograph  accommodated 
4  x  10  in.  photographic  glass  plates,  Kodak  type  103-F  and,  on  one  occa¬ 
sion,  type  103-0  emulsions  were  employed.  All  photographic  records  were 
developed  In  Kodak  0-19  developer  for  7  minutes  under  normal  development 
conditions. 
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Initiation  of  detonation  achieved  with  electronic  firing  units  , 

Is  synchronized  with  the  framing  camera  such  *:hat  the  camera  Is  ready  to 
record  the  explosive  event  at  a  predetermined  time  after  the  initiation 
of  the  detonation.  Rutherford  Model  A2  Time  Delay  Generators  were 
employed  to  provide  this  timing  synchronization.  Since  the  streak  camera 
Is  continuous  writing,  no  special  time  synchronization  is  required. 

A  typical  timing  sequence  for  a  jet  collision  experiment  is  illus¬ 
trated  in  Fig.  5,  the  indicated  time  delays  are  those  for  firing  no.  109. 
The  explosive  configuration  is  to  be  viewed  by  the  framing  camera  start¬ 
ing  at  a  time  when  the  jets  first  appear  from  behind  the  baffle  plates. 
The  time  required  for  the  framing  camera  rotating  mirror  to  be  in  posi¬ 
tion  after  a  firing  pulse  is  generated  (at  zero  time)  was  153  nsec. 
However,  the  time  interval  in  which  the  jets  first  appear  after  initia¬ 
tion  of  detonation  was  found  to  be  P  1  nsec.  Thus  the  firing  pulse  must 
be  delayed  132  nsec  by  the  delay  generator  before  it  initiates  the  explo¬ 
sive  detonator  for  proper  camera  synchronization.  The  generation  of  an 
82  nsec  delayed  pulse  supplied  to  an  exploding  wire  is  of  importance  for 
the  time-resolved  spectrum  for  this  firing  and  is  described  below. 


1.5  Observation  of  a  Single  Shaped  Charge  Jet 

To  obtain  a  description  of  the  jet  characteristics  and  to  experi¬ 
mentally  verify  timing  estimates,  a  single  shaped  charge  was  fired  in  air, 
(firing  no.  56).  Time  delays  were  established  so  as  to  have  detonation 


A  description  of  these  firing  uni  .  is  given  by  Walbrecht  (1959). 


of  the  explosive  completed  when  the  framing  camera  is  in  position  to 
expose  the  first  frame.  The  framing  camera  record  for  this  firing  is 
shown  in  Fig.  6.  The  view  is  normal  to  the  charge  axis  and  hence  to  the 
direction  of  detonation  and  jet  motion.  The  edge  of  the  explosive  charge 
Is  located  at  the  extreme  left-hand  side  of  each  frame.  The  direction  of 
detonation  and  jet  travel  is  from  left  to  right  in  the  figure. 

The  base  of  the  undetonated  charge  is  just  visible  on  the  extreme 
ltft  side  of  the  first  frame  in  the  figure.  In  the  second  frame  the 
luminous  reaction  products  extending  from  the  base  of  the  charge  are 
visible.  Thus  completion  of  the  charge  detonation  occurret  in  the  ti.T^ 
Interval  between  the  first  two  frames.  The  luminous  tip  of  the  jet  is 
first  observed  to  emerge  from  behind  the  expanding  reaction  products  at 
the  fourth  frame;  the  jet  is  seen  to  lead  a  reaction  products  at  subse¬ 
quent  times.  The  jet  tip  is  luminous  and  comet- 1  ike  in  appearance,  due 
to  _  lation.  The  narrow  jet  stem,  whose  diameter  is  of  the  order  of  2  mm 
near  the  tip,  is  readily  seen  in  later  frames.  The  jet  has  traversed  a 
distance  of  one  charge  diameter  from  the  base  of  the  charge  just  prior  to 
the  ninth  frame,  11  psec  ... f ter  the  first  frame.  Since  the  charge  was 
initiated  17  psec  prior  to  camera  synchronization,  a  shaped  charge  jet  is 
formed  and  travels  a  distance  of  one  charge  diameter  from  the  base  of  the 
charge  in  a  time  interval  of  28  psec  after  initiation  of  the  explosive 
train.  An  uncertainty  of  the  order  of  one  to  two  microseconds  is  ascribed 
to  this  time  estimate,  reflecting  the  uncertainty  in  framing  camera  syn¬ 
chronization  tine  and  in  the  delay  generator  setting. 

A  measurement  of  the  jet  velocity  has  been  obtained  from  the  streak 
camera  record  of  this  firing.  A  summary  of  jet  velocity  measurements 
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obtained  from  this  and  other  streak  camera  records  Is  given  In  Table  II. 
A  nominal  jet  speed  of  7-9  nw/psec  was  found  from  these  results.  This 
value  of  jet  speed  and  the  above  time  estimates  were  used  as  a  basis  for 
establishing  the  timing  sequence  In  succeeding  experiments. 


1.6  Jet  Collision  in  Air 

A  number  of  colliding  jet  experiments  were  performed  in  air  in 
order  to  observe  the  main  features  of  Jet  Impact,  framing  and  streak 
camera  records  of  the  collision  were  obtained.  In  addition,  a  time- 
integrated  spectral  record  of  the  collision  was  obtained  for  one  of  these 
firings,  the  intent  of  which  was  to  ascertain  the  presence  of  character¬ 
istic  radiation  from  the  copper  atoms  resulting  from  the  collision, 

A  framing  camera  record  for  one  of  these  experiments,  (firing  no. 

77) t  1*  shown  In  Fig.  7.  The  direction  of  Jet  motion  as  shown  is  verti¬ 
cal.  The  approaching  Jets  first  appear  from  behind  the  expanding  reaction 
products  In  the  second  frame.  A  conventional  explosively-shocked  argon 
source  provided  back  lighting  for  the  framing  camera  racord.  As  was 
observed  before,  the  jet  tips  are  luminous  and  comet-like  In  appearance. 
The  Jet  stem  is  readily  noticeable  in  later  frames.  Collision  occurs 
between  the  thirteenth  and  fourteenth  frames,  with  a  high  luminous  colli¬ 
sion  region  first  observed  at  the  fourteenth  frame.  This  region  expands 
In  time,  and  Its  extreme  brightness  persists  for  the  duration  of  the  film 
record.  The  growing  dark  areas  apparent  In  the  last  three  frames,  are 
presumed  due  to  expansion  cooi ing  and  interaction  with  the  encroaching 
explosive  reaction  products. 


A  streak  camera  record  of  another  jet  collision  experiment,  (firing 
no,  59),  Is  shown  In  Fig.  8.  On  the  left  side  of  the  figure  appears  an 
image  of  a  marker  scale,  which  physically  Is  located  at  the  explosive 
setup.  This  feature  provides  a  calibration  between  spatial  position  at 
the  explosive  event  and  transverse  distance  on  film.  The  time  scale  cali¬ 
bration  Is  known  in  terms  , of  the  streak  camera  rotational  speed.  The 
direction  of  increasing  time  is  from  left  to  right.  The  streak  camera 
entrance  slit  Is  oriented  in  the  direction  of  jet  motion.  The  approach¬ 
ing  jets  are  represented  by  two  narrow  light  trails  that  start  at  the 
left  side  of  the  figure;  the  light  comes  from  the  luminous  jet  tips. 

Jet  collision  is  characterized  by  the  intersection  of  these  trails  and 

subsequent  intense  exposure  that  lasts  for  about  50  usee.  A  region  of 

•*  -  • 

intense  luminosity,  approximately  centered  about  the  impact  point,  is 
seen  to  remain  roughly  localized  for  the  length  of  the  record.  Aside 
from  the  initial  jet  velocity  information  contained  in  a  streak  record, 
the  extremely  high  light  intensity  of  the  collision  products  is  apparent, 
even  when  viewed  through  a  narrow  streak  camera  entrance  slit. 

In  firing  no.  99,  a  t ime- 1 ntegrated  spectrum  of  a  jet  collision  In 
air  was  obtained.  Figure  9  shows  the  observed  spectrum  along  with  a 
densitometer  scan  and  a  comparison  copper  spark  spectrum.  The  overall 
characteristic  of  the  t ime- integrated  spectrum  Is  that  of  a  continuum 
with  a  few  weak  superposed  absorption  lines  over  a  wavelength  range  from 
3600  to  5000A.  The  lack  of  radiation  beyond  5000A  is  due  to  a  decrease 
in  sensitivity  of  the  Type-0  emulsion  beyond  this  wavelength.  Glass  and 
luclte  In  the  optical  path  of  the  spectrograph  are  primarily  responsible 
for  the  decrease  in  photographic  sensitivity  below  36OOA.  Absorption 
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lines  In  the  spectrum  were  Identified  by  both  direct  association  with 
known  copper  spark  lines  on  the  spectral  record,  (Shenstone,  1936;  1948), 
and  by  measurements  of  line  position  on  the  wavelength-calibrated  densi¬ 
tometer  trace. 

Table  Hi  presents  the  absorption  lines  observed  In  this  time- 
integrated  spectrum  along  with  associated  known  lines  of  copper.  Devia¬ 
tions  in  wavelength  from  the  known  copper  lines  are  predoml nent ly  due  to 
poor  definition  of  the  center  of  these  absorption  lines.  The  energy  of 
the  lower  excitation  level  (Moore,  1952)  for  the  observed  lines  Is  also 
given  In  the  table.  Most  of  the  levels  are  of  the  order  of  5  ev  above 
the  ground  state  of  copper.  Thus,  the  t I me- Integrated  spectrum  of  a  jet 
collision  In  air  shows  that  some  of  the  jet  material  is  vaporized,  and 
furthermore  a  fraction  of  these  gaseous  copper  atoms  reside  In  excited 
levels  energetically  several  electron  volts  above  the  neutral  copper  atom 
ground  state. 

Spectral  line  tabulations  show  a  multitude  of  characteristic  copper 
lines  for  the  spectral  region  covered  in  this  experiment  that  were  not 
observed  here.  Moreover,  these  lines  come  from  energetically  similar 
excited  states  and  have  comparable  tabulated  relative  intensities,  hence 
observation  of  these  lines  was  expected.  That  these  transitions  are  not 
seen  either  in  emission  or  absorption  is  attributed  to  at  least  two 
factors.  A  time- integrated  spectrum  indiscr iminately  records  light  from 
other  phases  and  parts  of  the  collision  luminosity.  This  may  have  blotted 
out  or  made  indistinct  any  indication  as  to  the  presence  of  these  missing 
spectral  transitions.  Also,  masking  of  collision  luminosity,  due  to  the 
presence  of  an  air  shock  created  at  the  time  of  collision  and  propagated 


along  with  the  collision  products  and  quenching  due  to  energy  transfer  to 
the  surrounding  atmosphere  are  of  importance.  Spectral  characteristics 
of  the  radiation  emerging  from  the  shocked  air  envelope  can  be  substan¬ 
tially  different  from  that  of  the  enclosed  luminous  collision  products 
(Bethe,  et  al,  1947).  In  the  extreme  case  of  an  optically  thick  shock 
zone,  the  observed  spectral  radiation  would  be  character i st ic  of  the 
shock  rather  than  of  the  jet  collision  products. 


i ,7  Jet  Collisions  in  Helium  Atmosphere 

Two  modifications  were  introduced  In  subsequent  experiments  that 
were  designed  to  alleviate  the  above  mentioned  difficulties.  First,  a 
direct-view  prism  was  added  to  the  streak  camera  entrance  optics  to 
enable  time-resolved  spectroscopic  observations  to  be  made,  {see  Appen¬ 
dix  B) .  Thus,  knowing  the  time  sequence  of  events  In  a  collision  experi¬ 
ment,  a  time-resolved  spectral  record  enables  one  to  separate,  In  time, 
radiation  from  the  collision  products  and  from  other  luminous  phases  of 
the  experiment.  Secondly,  to  eliminate  the  influence  of  air  shock  lumin¬ 
osity  and  shock  opacity,  experiments  were  performed  in  a  low  pressure 
helium  atmosphere.  Such  an  environment  should  be  transparent  to  collision 
radiation.  Optical  transitions  for  helium  in  the  visible  region  are 
associated  with  helium  energy  levels  that  are  about  20  ev  above  its  ground 
state.  In  addition,  no  appreciable  population  of  these  highly  excited 
levels  is  expected  during  the  course  of  the  experiment  since  the  density 
of  helium  atoms  et  the  residual  pressures  is  low,  of  the  order  of  lo'^  cm" 
The  vacuum  systems  employed  for  these  firings  were  described  above. 


(a)  Framing  Camera  and  Spectrograph  Observations 

in  firing  no.  i09,  two  shaped  charge  Jets  collided  In  helium  at  a 
residual  pressure  of  0.7  mm  of  Hg.  The  framing  camera  record  for  this 
experiment  is  shown  in  Fig,  JO.  The  view  is  through  one  of  the  rectangu¬ 
lar  ports  on  the  side  of  the  vacuum  chamber.  An  outline  of  the  port  is 
readily  discernable  in  iater  frames  of  the  record.  Initial  jet  motion  is 
horizontal  as  presented  in  the  figure.  The  jets  are  not  visible  prior  to 
impact  since  no  backlighting  was  used  in  this  firing. 

In  the  first  two  frames  only  light  from  the  explosive  reaction 
products  and  reflections  from  portions  of  the  support  structure  is 
observed.  Jet  impact  Is  first  seen  in  the  third  frame,  as  evidenced  by 
the  small,  centrally  located,  luminous  region.  A  brilliant  flash  of  light 
Is  seen  that  persists  for  about  four  frames  or  6  usee  after  first  impact 
light.  After  this  there  is  a  decrease  in  intensity  over  the  collision 
region  lasting  for  a  few  frames.  An  apparent  relighting  occurs  between 
the  tenth  and  eleventh  frame,  about  JO  usee  after  collision.  The  intens¬ 
ity  remains  high  throughout  the  remainder  of  the  record  although  the  view 
becomes  complicated  by  detonation  products  entering  the  collision  region 
and  possible  mechanical  failure  of  the  vacuum  chamber  at  these  later 
times . 

A  vertical,  luminous  line  structure  at  the  position  of  first  impact 
light  is  observed  on  the  sixth  frame.  The  light  intensity  and  definition 
of  the  strip  are  enhanced  in  following  frames.  Lukaslk,  et  el  (1964)  have 
presented  a  hydrodynamic  model  for  the  impact  of  two  high  speed  projec¬ 
tiles  in  which  they  suggest  that  a  radial ly-symmetric  expending  sheet, 
traveling  outward  from  the  collision  point,  is  formed  at  Impact.  Other 
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evidence  for  the  formation  of  a  sbeet-llke  structure  In  a  jet  collision 
can  be  seen  In  the  framing  camera  record  of  firing  no,  124,  shown  In 
Fig,  II  where  two  shaped  charge  jets  collided  In  an  expendable  luclte 
vacuum  chamber.  As  before,  Initial  jet  motion  Is  horizontal  as  presented 
In  the  figure.  Jet  Impact  Is  first  seen  In  the  sixth  frame,  as  Indicated 
by  the  small,  intensely  luminous  region  centered  between  the  explosive 
reaction  products.  As  early  as  the  second  frame  after  Initial  Impact  a 
vertical,  luminous  strip  is  apparent,  originating  at  the  point  of  Impact. 
The  Intensity  of  the  light  from  this  sheet  increases  rapidly  and  becomes 
the  most  prominent  feature  of  the  record  In  later  frames.  Further  evi¬ 
dence  of  a  mechanical  nature  for  the  formation  of  a  sheet-like  structure 
was  found  in  the  remains  of  the  support  structure  recovered  after  firing 
no.  125.  Figure  12  shows  the  recovered  structure.  Jet  collision  was  not 
head-on  and  so  the  sheet  was  not  normal  to  the  direction  of  jet  travel. 

It  appears  that  the  expanding  sheet  collided  with  one  of  the  support 
plates;  the  narrow  region  of  Impact  Is  most  obvious  In  the  photograph. 

A  time- integrated  spectral  record  of  a  jet  collision  in  helium, 
(firing  no.  109),  together  with  a  densitometer  trace  and  copper  spark 
reference  spectrum  is  shown  In  Fig.  13.  The  spectral  record  contains  many 
pronounced  absorption  lines  and  a  number  of  emission  lines  superposed  upon 
a  continuous  background.  The  correspondence  between  lines  in  the  colli¬ 
sion  spectrum  and  in  the  copper  reference  spectrum  is  readily  seen. 

A  type  103-F  photographic  emulsion  was  used  here,  extending  the  wavelength 
range  covered  to  about  700QA.  Many  more  distinct  spectral  lines  are  seen 
in  this  record  then  in  the  previously  shown  time- integrated  spectrum  for 
a  collision  in  elr.  The  comparatively  rich  spectral  content  of  this 
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record  Indicates  that  the  choice  of  helium  as  a  residual  gas  in  the  vacuum 
chamber  was  reasonable.  No  detailed  line  identification  was  made  from 
this  record  or  other  simiiar  t ime- integrated  records,  instead,  major 
emphasis  was  directed  towards  an  extensive  analysis  of  the  time-resolved 
spectral  record  of  firing  no.  109,  presented  below. 

(b)  Time-Resolved  Spectrum  of  a  Jet  Collision 

A  time-resolved  spectrum  of  a  jet  collision  in  helium  Is  shown  In 
Fig,  14.  The  wavelength  range  extends  from  about  4400  to  5300A.  This 
range  was  chosen  for  two  reasons.  A  number  of  prominent  neutral  copper 
lines  appear  in  this  interval,  (Shenstone,  1948)  which  encompass  all  of 
the  salient  features  of  the  energy  level  structure  of  neutral  copper. 

In  addition,  a  number  of  Ionized  copper  lines  also  oppear  in  this  range 
(Shenstone,  1936).  In  particular,  there  are  a  large  number  of  ionized 
copper  lines  in  the  4900  to  5000A  region  with  no  Interfering  neutral 
copper  lines.  Thus,  this  spectral  region  serves  as  an  excellent  Indi¬ 
cator  for  the  presence  of  singly  ionized  copper  in  the  collision  plasma. 

Included  on  the  film  record  are  zinc  and  helium  calibration  lines, 
a  copper  spark  reference  spectrum,  and  an  exploding  gold  wire  spectrum. 

The  exploding  wire  spectrum  served  as  a  time  fiducial  mark.  Its  purpose 
was  to  enable  one  to  correlate  spectral  information  on  the  record  with 
the  overall  time  sequence  of  events  in  the  experiment.  (The  timing 
sequence  for  firing  no.  109  was  previously  shown  in  Fig.  5.)  The  gold 
wire  was  fired  50  ^sec  before  initiation  of  detonation,  a  sufficient  time 
interval  in  advance  of  detonation  to  assure  that  light  from  the  exploding 
wire  does  not  overlap  light  from  the  jet  collision.  Allowing  for  a  28 
interval  In  which  the  shaped  charge  jets  are  formed  and  travel  e  distance 
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of  on a  charge  diameter,  a  total  estimated  time  between  initiation  of  the 
exploding  wire  end  Initial  jet  impact  is  78  psec.  The  observed  time 
interval  of  79  psee  on  the  spectral  record  is  in  agreement  with  these 
known  time  delays.  Thus,  one  can  be  certain  that  the  observed  time- 
resolved  spectrum  is  that  of  the  impact-generated  light. 

The  rate  at  which  the  entrance  slit  image  traverses  the  film  plane 
is  I  mm/psec.  Thus,  the  time  required  for  the  silt  image  to  traverse  its 
full  length  of  1.8  mm  is  1.8  psec.  Since  any  given  point  on  the  film  is 
exposed  only  in  the  Interval  during  which  the  entrance  slit  image  crosses 
that  point,  the  time  resolution  for  this  record  is  1.8  psec  (Smith,  1925; 
Seay,  et  al,  1961).  The  entrance  slit  of  the  streak  camera  was  centered 
at  the  collision  point  and  oriented  with  its  length  perpendicular  to  the 
direction  of  initial  Jet  motion,  i .e.  parallel  and  in  line  with  the 
expanding  collision  sheet.  A  projection  of  the  entrance  slit  at  the 
collision  point  extended  over  a  I .0  by  0. 1  In.  area  and  is  Indicative  of 
the  spatial  resolution  of  the  collected  optica)  radiation. 

The  Initial  portion  of  the  collision  spectrum  is  seen  to  be  a  very 
intense  continuum  that  lasts  for  some  5  to  6  psec.  Recall  that  the  fram¬ 
ing  camera  record  for  this  experiment,  described  above,  shows  an  initial 
brilliant  flash  of  light  for  about  the  sane  time  interval.  The  contl nuun. 
then  appears  to  fade  and  a  rich  emission  spectrum  is  observed.  At  about 
10  psec  a  relighting  phase  commences,  consisting  of  a  continuum  of  lower 
Intensity  than  the  initial  continuum.  At  about  13  psec  after  impact  the 
continuum  fully  extends  over  the  wavelength  range  covered.  Finally, 
amission  lines  are  seen  to  revert  to  absorption  lines  at  later  times. 

The  emission  lines  seen  after  6  psec  are  shown  below  to  be  mainly 


lines  of  neutral  copper  atoms.  There  Is,  In  add!  conclusive  evidence 

of  ionized  copper  in  the  spectral  range  between  about  4900  and  495QA. 

A  narrow  band  of  light  In  this  region  is  seen  to  emerge  from  the  Initial 
continuum.  The  Intensity  of  this  band  increases  at  later  times  and  sub¬ 
sequently  blends  into  the  relight  continuum.  There  are  about  twenty-five 
ionized  copper  lines  in  this  spectral  range  (Shenstone,  1948).  Since  the 
wavelength  resolution  of  the  time- resolved  spectrum  does  not  allow  these 
lines  to  be  seen  individual ly,  they  appear  as  an  unresolved  group  90  to 
6QA  wide.  Further  evidence  Indicating  the  presence  of  ionized  copper  in 
the  jet  collision  products  has  also  been  obtained  and  1$  presented  below. 

Thus,  high  speed  Jet  impact  results  in  a  varied  and  rich  spectral 
history.  Qualitative  evidence  exists  for  the  presence  of  ionized  and 
excited  neutral  atoms  In  the  Jet  collision  plasma.  A  detailed  analysis 
of  this  time-resolved  spectrum  is  pi asented  In  the  following  sections. 


Part  ft  Analysis  of  Time-Resolved  Spectrum 


A  study  of  the  Jet  collision  time-resolved  spectrum  of  firing  no. 

109  was  made.  The  spectrum  Is  rich  In  detail  and  contains  important 
Information  regarding  the  nature  of  the  jet  collision  process.  In  partic¬ 
ular  an  examination  of  the  shapes  and  shift,  of  several  spectral  lines  In 
a  time  interval  of  Importance  was  undertaken.  The  results  of  this  study 
lead  to  a  measurement  of  the  time-varying  neutral  and  ion  particle  dens¬ 
ity  in  the  plasma. 

Only  a  time  interval  from  about  6  to  13  psec  after  impact  was  exten¬ 
sively  investigated,  tittle  information  could  be  derived  from  the  early 
continuum  phase  of  the  spectrum  before  6  psec  since  the  film  record  was 
greatly  overexposed.  The  later  portions  of  the  record  after  13  psec  were 
not  considered  since  as  will  be  shown  later  the  collision  of  the  expand¬ 
ing  plasma  with  the  vacuum  chamber  wall  would  occur  at  these  times  and 
hence  would  greatly  complicate  the  analysis  of  the  observed  spectrum. 

Radiation  observed  by  the  time-resolved  spectrograph  originates 
from  a  particular  volume  element  in  the  collision  event  as  will  be  shown 
later.  The  projection  of  the  spectrograph  entrance  silt  onto  the  object 
plane  at  the  impact  region  determines  the  cross-sectional  dimensions  of 
this  volume.  The  effective  depth  of  the  radiating  volume  Is  small  in 
comparison  with  the  overall  size  of  the  collision  region.  Light  from 
Interior  points  of  the  collision  region  is  strongly  absorbed  upon  passage 
through  the  intervening  material  and  does  not  contribute  significantly  to 
the  emerging  radiation.  Only  a  thin  section  located  near  the  front  of 
the  expanding  products  Is  effective  in  producing  the  observed  light.  Thus 
the  particle  density  measurements  derived  from  the  time-resolved  spectral 


data  represent  spatial  mean  values  for  this  thin  region  under  observation. 
The  above  mentioned  time  resolution  of  the  spectrograph  will  introduce  a 
further  modification;  that  is,  the  numerical  density  estimates  will  be 
averaged  in  time  as  well  as  in  spatial  position. 


2.1  Spectral  Line  identification 

Densitometer  scans  of  the  time-resolved  spectrum  at  two  time  inter¬ 
vals  after  initial  impact,  and  of  the  copper  reference  spectrum  are  shown 
In  Fig.  15.  Each  scan  of  the  collision  spectrum  is  seen  to  consist  of  a 
background  conti nuum  with  superposed  emission  lines.  The  increased  intens¬ 
ity  of  the  continuum  at  11.5  psec  corresponds  to  the  relight  phase  of  the 
spectral  record.  About  sixteen  prominent  emission  peaks  can  be  distin¬ 
guished  at  both  time  intervals,  all  of  which  have  rather  broad  profiles, 
typically  of  the  order  of  10  to  20A  wide.  Comparison  of  these  with  the 
densitometer  scan  for  the  copper  reference  spectrum  illustrates  the  rela¬ 
tively  large  breadth  of  these  lines.  Both  the  background  continuum  and 
emission  lines  persist  throughout  the  time  interval  from  6  to  13  psec. 

With  the  aid  of  the  calibration  spectra  on  the  film  record  these 
emission  lines  have  been  identif  ed  mainly  as  neutral  copper  transitions 
and  are  tabulated  in  Table  IV.  All  known  prominent  transitions  of  neutral 
copper  for  the  spectral  range  covered  in  this  experiment  were  observed. 

Due  to  the  width  of  the  spectrum  lines  and  the  limited  wavelength  resolu¬ 
tion  of  the  system,  howe  _r,  some  of  the  expected  lines  are  not  completely 
separated.  This  is  reflected  in  Table  IV  by  showing  a  number  of  neutral 
copper  transitions  in  close  proximity  to  the  observed  lines.  The  upper 


level  designation  and  excitation  energy  (Moore,  1952)  for  the  spectral 
lines  with  the  highest  tabulated  intensities  (Shenstone,  19^8)  are  also 
shown  In  the  table. 

The  presence  and  growth  of  ionized  copper  1 ines  is  most  effectively 
i i lustrated  In  the  spectral  region  between  4900  and  4950A,  as  noted  previ¬ 
ously.  At  6.7  psec  there  Is  little  evidence  for  any  characteristic  Ion 
transitions  in  this  spectral  region,  whereas  at  11. 5  psec  a  pronounced 
increase  in  intensity  above  the  background  continuum  is  most  obvious. 

This  increased  intensity  represents  the  growth  of  an  unresolved  group  of 
about  twenty-five  copper  Ion  lines  in  the  50A  wide  region.  Also  the  simi¬ 
lar  appearance  and  growth  of  other  Ion  transitions  can  be  seen  in  the 
vicinity  of  the  Cul  4651  transition.  At  11.5  psec  a  distinct  hump  In 

this  line  is  apparent,  as  indicated  in  Fig.  15.  This  satellite  growth  is 

the  result  of  three  known  copper  ion  lines  at  about  4660A.  Spectral  lines 
profiles  of  this  emission  line  presented  later  will  show  the  growth  of 
these  ionized  lines  more  distinctly. 

An  energy  level  diagram  that  contains  the  more  intense  spectral 
transitions  given  in  Table  IV  is  given  in  Fig.  l6.  Also  included  are  the 
observed  groups  of  copper  ion  lines.  Several  indicated  transitions  in 
this  figure  are  of  concern  in  the  work  presented  below.  A  large  number 
of  both  neutral  and  ion  energy  levels  are  seen  to  have  been  excited.  For 
the  neutral  cnpper  lines,  some  of  these  levels  are  as  high  as  or  slightly 

above  the  copper  first  ionization  potential  at  7-7  ev  (Moore,  1952). 

Levels  that  appear  above  the  ionization  potential  are  associated  with  the 
anomalous  term  or  two-electron  excitation  configuration  of  the  copper 
atom,  (Shenstone,  1932;  White,  1934;  Herzberg,  1937).  These  anomalous 
terms  are  grouped  towards  the  right  in  the  portion  of  Fig.  16  that  corres 
ponds  to  neutral  copper.  The  normal  terms  appear  below  the  ionization 
potential.  (i.P.).  at  7.7  ev. 


For  the  observed  Ion  lines,  the  corresponding  upper  energy  levels  are 
about  25  ev  above  the  ground  state  of  neutral  copper.  Other  energy 
levels  for  both  neutral  and  ionized  copper  could  not  be  observed  since 
the  light  resulting  from  transitions  out  of  these  levels  does  not  He 
within  the  spectral  range  of  the  time-resolved  spectrograph . 


2.2  Spectral  line  Profiles  from  Time-Resolved  Spectrum 

A  number  of  profiles  for  the  prominent  emission  lines  In  the  time- 
resolved  spectrum  were  derived  from  the  densitometer  traces.  The  film 
record  was  scanned  at  1  jxsec  intervals,  from  5-7  to  12.7  psec  after  first 
impact  light.  Typical  densitometer  scans  were  previously  shown  in  Fig. 
15.  Horizontal  position  on  each  film  scan  is  related  to  the  wavelength. 
The  wavelength  calibration  was  established  for  each  record  using  the 
wavelength  i dent  5 f i cat  ions  previously  giver  in  Table  IV  Vertical  dis¬ 
placement  on  e3ch  scan  is  proportional  to  photometric  density  units. 

A  density  scale  calibration  was  established  by  inserting  calibrated  neu¬ 
tral  density  filter?  into  the  optical  train  of  the  densitometer  and 
recording  the  corresponding  response. 

To  obtain  the  profile  of  a  g’ven  emission  line  the  measured  densi¬ 
ties  were  first  converted  to  intensity  values  by  means  of  a  photometric 
film  emulsion  calibration  described  in  Appendix  C.  The  background  con¬ 
tinuum  level  for  each  densitometer  scan  was  established  at  wavelength 
regions  in  which  emission  lines  were  absent  and  this  background  level, 
after  being  converted  to  intensity  was  subtracted  from  the  total  observed 
intens i ty . 
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(a)  Photograph tc  Intensity  Calibration 

Several  Hurter  and  Orlffeld  (H-0)  emulsion  calibration  eurves  relat¬ 
ing  density  with  relative  Intensity  at  certain  wavelengths  were  experi¬ 
mentally  determined  for  the  Tri-X  emulsion  and  the  developing  process  used 
to  obtain  the  time-resolved  film  records.  A  description  of  the  calibra¬ 
tion  technique  and  experimental  results  are  given  in  Appendix  0.  A  rotat¬ 
ing  step-sector  wheel  method  was  employed.  A  light  source  and  step-sector 
wheel  were  mounted  in  front  of  the  direct-vision  prism  and  streak  camera 
entrance  optics.  The  interrupted  light  was  recorded  with  the  streak 
camera  after  being  dispersed  by  the  direct-vision  prism.  Zinc,  mercury- 
cadmium,  and  helium  spectral  lamps  were  used  as  light  sources.  Emulsion 
calibration  curves  were  constructed  at  discrete  wavelengths  for  the  vari¬ 
ous  spectral  lamp  sources.  The  results  were  checked  for  self-consistency 
as  suggested  by  Sawyer  and  Vincent,  (1943).  These  emulsion  calibration 
curves  were  found  to  be  in  agreement  with  similar  measurements  of  other 
investigators,  (Kodak,  1962;  Pressman,  1962).  Values  for  the  slope  of 
the  emulsion  calibration  curve,  y  ,  found  from  repeated  photometric  mea¬ 
surements  at  several  wavelengths  are  given  In  Table  V.  The  average  value 
of  y  Is  1.0  -  0.1  for  the  wavelength  range  studied.  This  value  was  used 
in  all  spectral  line  profile  calculations. 

(b)  Measurement  of  Spectral  Line  Profiles 

Horizontal  and  vertical  position  measurements  on  a  densitometer 
record  were  obtained  with  the  use  of  a  Gerber  Scanner  and  punched  onto 
cards  for  subsequent  analysis.  At  a  given  wavelength  position  vertical 
displacements  corresponding  to  both  the  total  density  and  to  the  back¬ 
ground  continuum  were  measured.  Approximately  fifty  to  one  hundred  points 
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were  lifted  to  reproduce  a  line  profile.  In  addition,  the  vertical  dis¬ 
placements  of  the  photographic  density  calibration  steps  on  the  densi¬ 
tometer  reccrd  were  recorded.  This  information  was  used  to  establish  a 
vertical  height  to  density  scale  calibration.  For  a  wavelength  scale 
calibration,  the  horizontal  positions  at  the  maximum  points  of  each 
prominent  emission  line  were  also  recorded. 

Spectral  line  profiles  were  constructed  from  these  densitometer 
measurements.  A  listing  of  the  IBM  1620  Fortran  program  together  with 
a  brief  description  and  the  results  of  a  typical  calculation  are  given 
In  Appendix  C.  It  was  found  in  all  cases  that  the  observed  density  values 
fell  with!'  the  linear  portion  of  the  emulsion  calibration  curve.  Hence 
a  density  value,  0  ,  and  its  corresponding  relative  Intensity,  J  ,  are 
related  by,  (Mees,  1946;  Sawyer,  1951) 

D  «  const  +  y  log  J  (2-1) 

where  y  Is  the  slope  of  the  emulsion  calibration  curve  In  the  linear 
region.  Equation  (2-1)  was  then  used  to  evaluate  the  relative  background 
and  total  intensities,  Jg(X)  and  Jj(\)  respectively,  at  the  given 


wavelength,  X  .  That  is, 

Jg(x)  -  const  x 

°o/Y 

10  B 

(2-2) 

J^.(X)  *  const  x 
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10  1 
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where  0D  and  D_  are  the  corresponding  density  values.  The  spectral 

D  T 

line  intensity,  J(X)  ,  is  simply  the  difference  between  the  total  and 
background  intensities 
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It  was  found  convenient  to  normalize  the  line  intensities  to  the 

maximum  computed  value  J  at  a  wavelength  \  for  the  spectral  line 

o  o 

under  study,  and  also  to  deal  with  wavelength  deviations,  <  ,  from  the 
wavelength  of  the  maximum,  viz. 

e  -  X  -  \e  (2-5) 

Thus  the  spectral  line  profile  program  yields  the  relative  line  intens¬ 
ity,  J(«)/JQ  ,  at  a  wavelength  deviation,  e  ,  for  each  emission  line  on 
a  given  densitometer  trace. 


2.3  Line  Broadening  Mechanisms 

In  general,  a  number  of  causes  result  in  the  broadening  of  a  spec¬ 
tral  line.  They  can  be  classified  as  follows,  (Mitchell  and  Zemansky, 
1934;  White,  1934;  Aller,  1 953 *  Penner,  1959;  Breene,  1961). 

A.  Doppler  broadening 

B.  Foreign  atom  collision  broadening 

C.  Natural  line  breadth 

D.  Resonance  broadening 

E.  Stark  broadening 

F.  Instrumentation  broadening 

The  first  two  phenomena  can  be  disregarded  as  giving  rise  to  the 
observed  broadening  for  this  work  since  their  contributions  to  the  line 
width  are  inconsequential.  Doppler  broadening  effects  due  to  the  thermal 
motion  of  the  emitting  atoms  are  usually  small  and  limited  to  the  narrow 
center  or  core  of  the  spectral  line.  The  halfwldth  of  a  Doppler-broadened 
line  is  proportional  to  the  square  root  of  the  temperature  of  the 


radiating  gas.  A  Doppler  halfwidth  of  the  order  of  several  angstroms 

observed  for  the  copper  spectrum  lines  In  the  time-resolved  spectrum 

7  8  o 

would  Imply  an  unrea 1 1st  leal 1y  high  temperature  of  1£r  -  10  °K.  Also, 

In  view  of  the  low  density  of  foreign  atoms  In  the  vicinity  of  the  jet 
collision,  (the  residual  helium  gas  In  the  vacuum  chamber),  foreign  atom 
collision  broadening  Is  extremely  small  and  Is  thus  ineffective  In  pro¬ 
ducing  the  large  observed  line  widths. 

The  natural  line  breadths  for  most  of  the  observed  copper  transi¬ 
tions  have  also  been  shown  to  be  comparatively  small,  (Allen,  1932). 
However  several  copper  lines  are  known  to  possess  unusually  large  natural 
widths,  of  the  order  of  1A,  due  to  autoion Izat Ion  of  the  upper  energy 
level  of  the  spectral  transition,  (Shenstone,  1932).  A  few  of  these  lines 
appear  in  the  time-resolved  spectrum  and  their  observed  halfwidths  may  be 
due  In  part  to  their  large  natural  line  width. 

In  the  presence  of  atoms  that  ar*.  identical  to  the  emitting  atoms 
a  spectral  line  will  be  resonance  or  self-broadened.  This  type  of  broad¬ 
ening,  due  to  the  strong  interaction  between  an  emitting  atom  and  its 
nearby  perturbing  atoms,  can  result  in  line  breadths  of  the  order  of  the 
observed  breadths  if  the  density  of  identical  atoms  is  sufficiently  high. 

Similarly,  Stark  broadening,  due  to  the  i ntermolecular  electric 
fields  of  electrons  and  ions  in  the  collision  plasma,  must  also  be  con¬ 
sidered  since  Stark-broadened  lines  are  also  known  to  possess  broad  pro¬ 
files.  In  addition,  a  wavelength  displacement  of  a  spectral  line  can 
occur  due  to  the  Stark  splitting  and  shift  of  the  energy  levels  of  the 
transition  under  the  influence  of  the  i ntermolecular  electric  fields. 

Finally,  distortions  or  modif Icat Ions  of  a  spectral  line  profile 
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due  to  the  spectrograph  1 c  Instrumentation  must  be  taken  Into  account. 

It  will  be  shown  In  the  following  sections  that  the  principal  mech¬ 
anism  giving  rise  to  the  observed  spectral  line  shapes  Is  resonance-broad¬ 
ening.  This  is  achieved  by  matching  resonance-broadened  line  profiles  to 
several  neutral  copper  transitions  In  the  time-resolved  spectrum  of  firing 
no.  109,  where  modification  of  the  form  of  the  resonance  profile  due  to 
spectrograph  I c  Instrumentation  has  been  included. 

(a)  Resonance-Broadened  Line  Profiles 

The  profile  of  a  resonance-broadened  spectral  line  as  derived  from 
either  line  broadening  impact  theory  or  statistical  theory  Is  given  by 
the  Lorentz  formula,  (Margenau  and  Watson,  1936;  Breene,  1961) 

J(v)  -  const/((v-vo)2  +  (A/2)2)  (2-6) 

where  vq  Is  the  frequency  at  the  center  of  the  symmetrically-broadened 
line  and  A  is  the  halfwldth  of  the  line  in  frequency  units.  Normal¬ 
izing  the  Intensity  J(v)  to  Its  maximum  vaiue  -Kv0) 

J(v)/J(vo)  -  (A/2)2/((v-vo)2  +  (A/2)2)  (2-7) 

The  halfwidth  is  linearly  related  co  the  number  density  of  perturb¬ 
ing  atoms  of  the  same  kind  as  the  emitter,  N  ,  and  to  the  oscillator 
strength,  f  ,  of  the  spectral  transition  as, 

A  «  const  (e2f/mv  )N  (2-8) 

o 

where  e,m  are  the  electronic  cnarge  ana  mass  respectively.  ihe  constant 
In  eq  (2-8)  has  been  evaluated  from  both  statistical  and  impact  broadening 

The  term  haifwidth  as  used  here  refers  to  the  full  width  of  the  line  at 
half  of  Its  maximum  intensity. 
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theory.  Comparative  results  given  by  Breene  are  shown  in  Table  VI  and 
are  seen  to  be  in  dose  agreement. 

(b)  instrumentation  Broadening 

Before  attempting  to  fit  resonance-broadened  profiles  to  the  emis¬ 
sion  lines  in  the  time-resolved  spectrum  of  firing  no.  109  modifications 
due  to  the  spectrographic  instrumentation  must  be  taken  into  account. 

In  particular  if  the  image  of  the  entrance  slit  Is  comparable  in  size  to 
the  actual  spectral  line  width  then  the  resultant  line  profile  will  be 
influenced  by  the  characteristics  of  the  spectrograph.  Consider  a  spec¬ 
trograph  entrance  slit  of  finite  width  irradiated  by  a  light  source  that 
contains  several  discrete  frequency  components  Vj  .  If  the  slit  Is  of 
sufficient  width  so  that  diffraction  effects  can  be  neglected  then  silt 
images  of  width  Aj  ,  will  appear  in  the  image  plane  of  the  spectrograph 
at  positions  corresponding  to  the  frequencies  Vj  .  The  light  intensity 
of  the  Vj-component  will  be  zero  outside  of  each  strip  Aj  and  of  magni¬ 
tude  J(Vj)  within  the  strip,  (Baly,  1905;  Slater,  1925;  Forsythe,  1932; 
Clark,  l 960 ;  Brugel,  1962).  if  there  is  no  overlapping  of  slit  images 
then  the  distribution  of  light  intensity  along  the  image  plane  is  simply 
described  by  the  array  of  strips  of  width  Aj  and  intensity  J(v.)  . 

When  the  several  slit  images  overlap,  the  component  intensities  are 
superposed.  The  resultant  light  intensity  at  a  frequency  v  in  the 
regie.,  wf  overlap  is 

t(v)  -  Z  J(v.)H(A,/2  +  (v,-v))H(A,/2  -  (v,-v))  (2-9) 

where  H  is  the  unit  step  function.  The  summation  extends  over  the 
number  of  intensity  components  whose  slit  image  width  includes  the  fre¬ 
quency  position  v  .  In  the  limit  as  the  number  of  spectral  components 
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Increases  Indefinitely,  the  resultant  Intensity  becomes  (Oalton,  t963) 


vhA/2 


*<v) 


J 


J(v')  dv' 


v-A/2 


(2-10) 


Thus  the  resultant  light  intensity  at  a  given  frequency  Is  simply  the 
Integrated  or  average  intensity  over  a  narrow  frequency  band. 

The  modification  of  a  resonance-broadened  profile  Is  obtained  by 
evaluating  eq  (2-10)  where  the  intensity  distribution  of  the  light  source 
J(v')  Is  given  by  the  Lorentz  formula,  eq  (2-7) .  The  resultant  Inte¬ 
grated  Lorentzian  function  Is 


4{v) 


r  i  rv“v«  +  A/ 2, 
J(v0)(A/2)  (tan’  {— 7575 - } 


tan 


•'  Pbi^} } 


(2-11) 


Equation  (2-10)  can  be  expressed  in  a  more  general  form  as 


*(v) 


J* S(v,v')  J(v’)  dv 


(2- 10a) 


The  slit  function  S(v,v')  represents  the  weighted  contribution  to  the 
resultant  intensity  at  the  frequency  v  from  radiation  of  frequency  v', 
(Dennison,  1928;  Nielsen,  et  al,  1944;  King  and  Emslie,  1951;  Brodersen, 
1953;  Penner,  1959).  Triangular  and  trapezoidal  slit  function  shapes 
are  contnonly  used  for  spectrographs  with  both  entrance  and  exit  slits. 
Gaussian  shapes  are  used  to  approximate  diffraction  effects  associated 
with  small  slits.  A  rectangular  slit  function  shape  representing  a  limit¬ 
ing  form  of  a  trapezoidal  shape  is  useful  for  spectrographs  with  only  an 
entrance  slit.  A  rectangular  slit  function 

S(v,v')  *  1  v-A/2  <  v'  <  v+A/2 

$(v,v')  »  0  otherwise 

has  been  used  in  this  work. 
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Normalizing  #(v)  to  Its  maximum  value  #(v  )  one  has 


-1  {'"wo  **/Z]  .  -I  lv-vo  -  ^ 

t(v)  .  tan  l'"A/8 - J  ~  tan  l"  A/2  J 


»<V 


2  tan 


-I 


(2-12) 


The  halfwidth  of  the  integrated  Lorentzian  profile  5  ,  evaluated  directly 
from  eq  (2-12),  Is  found  to  be 

»  A2  +  A2  (2-13) 

Thus  the  use  of  a  finite-width  entrance  slit  will  alter  the  shape  of  a 
resonance-broadened  profile  and  further  increase  the  observed  line  half¬ 
width. 


If  two  or  more  resonance  broadened  lines  are  not  completely  resolved 
the  resultant  profile  is  described  by 

v-v.  ♦  A/ 2-* 


*N)  -  SJ(v,)(A,/8)  {tan  {— JJ— } 


-  } 


(2-14) 


where  is  the  halfwidth,  the  center  frequency  of  the  ith  par¬ 

tially  resolved  line,  and  J(vj)  represents  the  relative  intensity  of 
the  ith  component.  The  summation  includes  all  the  partially  resolved 
broadened  lines  that  can  contribute  to  the  resultant  line  profile. 

The  width  of  the  slit  image  is  related  to  the  entrance  slit  width 
and  the  optical  parameters  of  the  spectrograph.  Evaluations  of  the  slit 
image  width  have  been  made  for  the  direct  vision  prism-streak  camera 
combination,  details  of  which  are  presented  in  Appendix  E.  Essentially 
what  is  done  is  to  construct  an  image  of  the  entrance  slit  by  following 
the  path  of  several  principal  rays  through  the  optical  system.  This  is 
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repeated  at  several  wavelengths  In  the  spectral  range  covered  by  the 
time-resolved  spectrum  to  obtain  the  variation  of  slit  image  width  with 
wavelength.  The  results  of  the  calculations  are  illustrated  in  Fig.  17- 
Two  families  of  curves  are  shown  which  correspond  to  slightly  different 
values  of  the  entrance  slit  width.  Changes  in  the  slit  image  width  as  a 
result  of  small  deviations  about  the  nominal  angle  of  incidence  are  also 
i nc luded . 

Experimental  values  of  the  slit  image  width  have  been  independently 
determined  from  integrated  Lorentzian  profile  fits  to  several  lines  in 
the  time-resolved  spectral  record  as  described  below,  and  are  indicated 
in  the  figure.  The  accuracy  of  the  slit  image  widths  determined  in  this 
manner  is  limited  by  the  sensitivity  of  the  line  profile  fitting  process, 
(see  section  2.5).  Uncertainties  in  slit  width  values  are  indicated  by 
the  error  flags  in  Fig.  17*  In  view  of  the  fact  that  nominal  values  for 
the  optical  parameters  of  the  direct  vision  prism-streak  camera  combina¬ 
tion  were  used  in  the  calculations,  the  agreement  between  the  calculated 
and  experimental  values  is  good. 

(c)  tine  Profile  Fits  to  Selected  Transitions 

Integrated  Lorentzian  profiles  were  matched  to  several  of  the 
observed  neutral  copper  lines  In  the  time-resolved  spectrum.  This  estab¬ 
lished  that  resonance  broadening  is  the  principal  mechanism  responsible 
for  the  breadth  and  shape  of  the  spectral  lines.  These  profiles  were 
computed  from  either  eq  (2-12)  or  eq  (2-14).  Fortran  listings  of  the 
programs  for  both  the  single  component  and  multicomponent  profiles  are 
given  in  Appendix  F,  together  with  a  brief  description  of  the  programs. 


JJ 


The  procedure  for  calculating  a  single-component  line  profile  is  as 
follows.  A  time  in  the  spectral  record  for  which  the  observed  line  half¬ 
width  is  smallest  is  first  examined.  The  resonance  halfwidth,  A  ,  is 
then  computed  from  eq  (2-13)  for  an  assumed  silt  width  A  ,  and  an  assumed 
total  halfwidth  6  .  A  line  profile  described  by  eq  (2-12)  is  then  con¬ 
structed,  By  appropriately  varying  both  the  slit  width  and  total  half¬ 
width  values,  a  profile  that  best  matches  the  observed  line  shape  is 
generated.  This  results  in  a  measure  of  the  slit  width  at  the  wavelength 
under  consideration  and  an  evaluation  of  the  true  resonance  halfwidth  of 
the  line  at  the  particular  time  studied.  Line  profile  fits  for  the  same 
line  at  other  times  are  then  achieved  by  varying  just  the  total  halfwidth, 
keeping  the  slit  width  fixed,  until  a  best  match  with  the  observed  line 
shape  is  obtained. 

!f  contributions  to  the  observed  line  shape  due  to  partially 
resolved  neighboring  lines  are  taken  into  account,  a  mu l t i component  line 
profile  given  by  eq  (2-14)  is  employed.  In  this  case  values  for  the  res¬ 
onance  halfwidth  A.  ,  and  relative  Intensity  J(vj)  ,  of  each  spectral 
line  that  contributes  to  the  observer!  line  are  needed.  First  estimates 
for  these  quantities  are  established  by  matching  a  single-component  line 
profile  to  the  most  prominent  spectral  lire.  The  difference  between  the 
observed  line  shape  and  the  single-component  profile  fit  indicates  the 
relative  amplitude  and  width  of  the  nearby  partially-resolved  line  or 
lines.  A  line  profile  is  then  constructed  with  these  data.  Appropriate 
variations  in  all  the  resonance  halfwidth  and  relative  intensity  values 
are  then  made  in  order  to  generate  a  profile  that  best  matches  the 
observed  line  shape. 
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Single-Component  Profiles 

Integrated  Lorentzian  fits  to  the  observed  transitions  Cul  4587, 

4651  and  5153  were  generated  in  the  above  manner,  at  1  psec  intervals  over 
a  time  period  from  about  6  to  13  |J.sec  after  initial  jet  Impact.  The  Cul 
4587  line  was  specifically  chosen  since  there  are  no  nearby  transitions  of 
neutral  or  ionized  copper  in  the  wavelength  region  of  the  observed  broad¬ 
ened  line.  Hence  a  s I ng i e-component  integrated  Lorentzian  profile  fit  is 
m  1st  appropriate  for  this  line.  However,  the  natural  width  of  the  Cul  4587 
transition  is  known  to  be  relatively  large  due  to  autoionization  of  its 
upper  level,  (Allen,  1932;  Shenstone,  1932).  Thus,  the  observed  line  pro¬ 
file,  aside  from  instrumentation  broadening,  is  the  result  of  a  combina¬ 
tion  of  natural  braodentng  and  resonance  broadening  interactions. 

The  Cul  4651  line  was  considered  next  since  it  is  a  component  of  the 
same  multiplet  as  the  Cul  4587  transition,  but  is  associated  with  an  upper 
energy  level  that  does  not  exhibit  autoi on i zat ion .  Thus,  natural  broaden¬ 
ing  of  this  line  has  little  influence  cn  the  overall  observed  line  shape. 
Although  several  partially  resolved  copper  lines  are  present  and  spaced 
about  IOA  on  both  sides  of  the  Cul  465'  line,  only  single-component  inte¬ 
grated  Lorentzian  profile  fits  were  performed.  In  determining  the  best 
match  to  the  observed  line  shape  in  this  case,  attention  was  directed 
towards  the  fitting  over  a  wavelength  region  that  encompassed  no  neighbor¬ 
ing  lines  rather  tuan  the  overall  line  shape. 

Integrated  lorentzian  profiles  to  the  Cu!  5153  line  were  also  made. 
This  transition  was  specifically  chosen  since  it  originates  from  a  lower 
energy  level  and  involves  a  configuration  in  which  only  one  electron  is 
Involved  in  an  optical  transition  as  distinct  f>uii<  the  two-electron 
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configuration  of  the  previous  two  lines.  As  with  the  Cul  4651  line,  the 
best  fit  determination  was  made  over  wavelength  regions  that  contained  no 
nearby  part lal ly-resolved  spectral  lines.  Fo‘  all  three  lines  the  slit 
width  contribution  to  the  observed  line  shapes  was  first  established  as 
previously  indicated.  Values  of  the  slit  widths  determined  and  used  in 
these  profile  fits  were  shown  in  Fig,  17. 

The  integrated  Lorentzian  fits  to  these  lines  are  presented  in 
Figs.  18,  19,  and  90.  Observed  line  shapes  in  the  time-resolved  spectrum 
are  represented  by  the  data  points  while  the  solid  curves  show  the  fitted 
integrated  Lorentzian  profiles.  The  corresponding  time  after  Initial  jet 
impact  is  shown  in  each  figure.  The  relative  position  and  amplitudes  of 
known  spectral  transitions  of  copper  that  fall  within  the  spectral  range 
studied  are  also  indicated. 

The  Integrated  Lorentzian  profiles  shown  in  Fig.  18  adequately  rep¬ 
resent  the  observed  Cul  4587  line  shape  over  the  time  interval  studied. 
Irregularities  in  the  observed  line  shapes  at  the  later  time  intervals 
are  believed  to  be  the  result  of  impact  of  impact  between  the  expanding 
jet  collision  sheet  and  the  vacuum  chamber  window. 

The  observed  Cul  465 1  line  shapes  and  corresponding  integrated 
Lorentzian  profiles  shown  in  Fig.  19  are  also  seen  to  be  in  good  agreement 
over  a  substantial  portion  of  the  theoretical  line  profile.  Of  particu¬ 
lar  importance  is  the  appearance  and  growth  of  three  ionized  copper  trans¬ 
itions,  (Cull  4660,  466 1 ,  4663),  which  are  present  as  an  unresolved  cluster 
In  the  observed  line  shapes.  The  relative  intensity  of  these  unresolved 
lines  Is  seen  to  increase  rapidly  at  the  later  times.  The  appearance  cf 
the  partially  resolved  Cul  4643  line  during  early  times  is  also  of  interest. 
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Contr {buttons  from  this  transition  to  the  observed  line  profiles  are  con¬ 
sidered  In  the  following  section. 

The  observed  line  shapes  for  the  Cul  5153  transition  shown  in  Fig. 
20  are  also  seen  to  be  adequately  represented  by  an  integrated  Lorentzian 
profile.  Contributions  to  the  line  shape  from  nearby  broadened  neutral 
and  Ion  spectral  lines  are  not  as  pronounced  as  In  the  case  of  the  Cul 
46$1  line. 

Part  (a)  of  TabJe  VII  lists  values  of  the  resonance-broadened  half¬ 
width  A  ,  obtained  from  the  above  integrated  Lorentzian  profile  fits. 
Accuracies  for  the  halfwidths  found  by  this  procedure  are  limited  by  the 
sensitivity  of  the  line  fitting  process  and  are  of  the  order  of  -10  per¬ 
cent,  (see  section  2.5). 

Multicomponent  Profiles 

Single-component  profile  fits  to  the  Cul  4651  transition  shown  in 
Fig.  19  indicate  that  the  partially  resolved  Cul  4643  line  noticeably 
contributes  to  the  observed  line  shape  at  early  time  intervals.  Fits  to 
the  observed  line  shapes  with  a  multicomponent  profile  fit  that  included 
both  the  Cul  4643  and  4651  transitions  were  made.  Values  for  the  reso¬ 
nance  halfwidth  of  the  Cul  4651  line  determined  from  the  above  single¬ 
component  profile  results  were  initially  used.  Estimates  for  the  reson¬ 
ance  halfwidth  and  relative  intensity  of  the  Cul  4643  line  were  first 
obtained  from  a  comparison  of  the  s i ng I e- component  profiles  and  the 
observed  line  shapes.  By  appropriate  variations  of  the  resonance-broad¬ 
ened  halfwidth  values  and  relative  intensity  a  best  match  to  the  observed 
line  shape  over  an  extended  wavelength  range  was  achieved.  Multicompo¬ 
nent  fits  at  5.7>  6.7,  and  7-7  psec  were  made,  the  resulting  fit  at 
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6.7  psec  is  shown  in  Fig.  21.  Agreement  between  the  observed  line  shapes 
and  the  fitted  curves  is  good  in  a t !  three  cases.  The  resonance-broadened 
halfwidth  values  and  relative  intensities  for  the  Cui  4643  and  4651  trans¬ 
itions  used  to  obtain  the  best  multicomponent  profile  match  are  listed  in 
Table  VII,  part  (b)  and  Table  VIII  respectively.  The  intensity  of  the 
Cut  4651  line  was  arbitrarily  chosen  as  unity  in  all  three  profile  fits. 

Resonance-broadened  halfwidths  for  several  other  observed  lines  in 
the  time-resolved  spectrum  have  also  been  obtained  from  further  multicom¬ 
ponent  line  fits.  Halfwidth  values  and  relative  intensities  for  the 
partially  resolved  components  Cut  4525,  4531,  and  4540  at  three  times  are 
aiso  shown  In  part  (b)  of  Table  VII  and  in  Table  VI!  1.  The  relative 
Intensity  of  the  principal  Cul  4540  line  was  taken  as  unity.  The  result¬ 
ant  multicomponent  fits  to  the  observed  line  shapes  are  shown  in  Figs.  22 
and  23.  The  study  of  these  lines  is  of  importance  in  subsequent  work  and 
is  discussed  in  the  following  sections. 

In  conclusion,  an  integrated  Lorentzian  profile  has  been  shown  to 
adequately  represent  the  shape  of  a  number  of  broad  and  overlapping  lines 
in  the  time-resolved  spectral  record.  When  the  natural  line  breadth  can 
be  neglected,  the  observed  line  breadths  are  the  results  of  a  combination 
of  instrument  broadening  and  resonance  broadening.  Slit  image  widths 
obtained  from  these  line  fits  have  been  previously  shown  to  be  in  agree¬ 
ment  with  calculated  slit  width  values. 

(d)  Influence  of  Plasma  Ions  on  Spectral  Characteristics 

The  presence  of  charge  particles  in  the  jet  collision  plasma  will 
Influence  the  nature  or  the  spectral  line  radiation  in  general  as  a  broad¬ 
ening  and  shift  of  a  spectral  line.  These  effects  are  related  to  the 
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Stark  splitting  and  displacement  of  the  energy  levels  of  the  spectral 
transition.  The  splitting  and  shift  are  In  turn  a  consequence  of  the 
fluctuating  Intermolecular  electric  fields  created  by  the  charged  par¬ 
ticles  in  the  vicinity  of  the  emitting  atom.  The  Stark  effect  for  several 
neutral  copper  transitions  in  the  wavelength  region  under  consideration  is 
known.  The  pattern  and  magnitude  of  the  line  splitting  with  an  applied 
electric  field  has  been  measured  by  TakamSne  (1919),  some  results  of  which 
are  given  in  Table  IX.  Three  of  the  lines  are  split  into  two  components 
asymmetrically  displaced  with  respect  to  the  original  wavelength. 

An  atom  in  an  excited  state  is  influenced  simultaneously  by  colli¬ 
sions  with  both  identical  neutral  atoms  and  charged  particles  in  the  jet 
collision  plasma.  Thus  the  spectral  line  profile  for  the.  transition  from 
this  excited  state  will  in  general  reflect  both  broadening  effects  and 
would  be  described  in  a  complicated  way  by  the  resonance  and  charged  par- 
tical  interactions.  Previous  results  however  have  shown  that  instrumenta¬ 
tion  and  resonance-broadening  are  the  predominent  mechanisms  affecting 
those  lines  studied; in  particular,  line  profile  fits  t »  the  Cut  5153 
transition,  which  exhibits  a  Stark  splitting  and  thus  could  be  expected 
to  show  Stark  broadening,  were  given  in  fig,  20.  isence  superposed  Stark 
broadening  contributions  in  the  time-resolved  spectrum  are  not  apparent. 
This  result  is  not  surprising  in  view  of  the  magnitude  of  the  Stark  split¬ 
ting  for  this  spectral  line.  Comparative  Stark  splitting  for  hydrogen  for 
which  Stark  broadening  effects  have  been  well  established.  (Edels  and 
Craggs,  1951;  Griem,  1 95^+ ;  Henkel,  195^;  Margenau  and  Lewis,  1959)  are 
several  orders  of  magnitude  greater  than  the  Stark  splitting  for  copper. 
Aiyamovskii  and  Kitaeva  (1959),  have  experimentally  generated  ionic 
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Stark- broadened  H_  lines  with  a  haifwidth  of  about  50A  correspond ing  to 

P 

an  Intermoiecuiar  electric  field  of  the  order  of  100  kv/cm.  Similarly 
Holtsmark  and  Trumpy  ( 1 925) >  have  found  Stark-broadened  halfwidths  of 
about  1A  with  an  Intermoiecuiar  field  of  the  order  of  25  kv/cm  for  the 
Cui  4023  and  4062  transitions.  The  magnitude  of  the  Stark  splitting  for 
these  copper  lines  is  larger  than  that  for  the  lines  shown  in  Table  IX  by 
approximately  a  factor  of  ten.  Hence  line  broadening  contributions  due 
to  Ions  and  electrons  In  the  jet  collision  plasma  are  expected  to  be  neg¬ 
ligibly  small  for  all  the  observed  lines  compared  to  the  predominant 
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resonance  and  instrumentation  broadening. 

The  asymmetric  splitting  of  the  transitions  given  in  Table  IX  sug¬ 
gest  that  line  shifts  may  be  present  In  the  observed  spectrum.  Since 
these  lines  are  instrument  and  resonance-broadened  the  individual  Stark 
components  are  not  resolved.  However  superposition  of  the  shifted  Stark 
components  will  result  in  a  spectral  line  profile  that  has  a  net  wave¬ 
length  displacement  in  the  direction  of  the  Stark  component  shifts.  The 
net  spectral  line  shift  can  then  be  related  to  the  component  splitting  „ 
and  hence  to  an  estimate  for  the  intermoiecuiar  electric  field  that  is 
responsible  for  the  observed  shift.  Evidence  for  a  net  wavelength  shift 
affords  an  excellent  example  of  the  influence  of  charged  particles  on  the 
observed  jet  collision  spectrum  since  resonance-broadened  lines  exhibit 
no  wavelength  shift  (Margenau  and  Watson,  1936;  Foley,  1946).  Thus  in 
contrast  to  the  situation  for  line  broadening,  resonance  interactions  and 
charged  particle  effects  are  decoupled  in  a  line  shift  measurement. 

Further  evidence  for  this  conclusion  Is  contained  in  the  results  of  the 
next  section. 
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Of  the  several  lines  given  above  for  which  Stark  effect  information 
is  known  the  Cut  4531  1 ine  appeared  to  be  the  best  choice  for  further 
study.  This  line  occurs  as  a  partially  resolved  component  of  the  Cul  4540 
transition  along  with  another  nearby  line,  Cui  4525.  ho  Stark  effect  data 
has  been  given  by  Takamlne  for  these  companion  lines.  It  is  presumed  that 
the  Stark  splitting  for  these  lines  is  much  smaller  than  that  for  the  Cul 
4531  line,  and  hence  was  not  observed  In  Takamine's  experiments.  With 
the  assumption  that  these  companion  lines  exhibit  negligible  shift  in  the 
time-resolved  spectrum,  the  position  and  shift  of  the  resonance  broadened 
Cul  4531  line  with  respect  to  these  two  adjacent  lines  can  be  measured. 
Furthermore,  due  to  the  close  proximity  of  these  three  lines  a  line  shift 
measurement  would  be  uneffected  by  the  uncertainty  and  error  in  the  wave¬ 
length  scale  associated  with  the  observed  spectral  records. 

For  the  Cul  5153  transition  no  nearby  comparison  lines  were  found 
in  the  spectral  data.  Hence  a  line  shift  determination  over  extended 
wavelength  intervals  would  be  necessary  and  consequently  significant  wave¬ 
length  scale  errors  would  be  included  in  a  line  shift  measurement.  The 
Cul  5218,  5220  transitions  which  appear  as  one  unresolved  emission  line  on 
the  spectra!  record,  were  not  considered  since  Stark  effect  data  for  Cul 
5220  appears  incomplete. 

The  resonance-broadened  Cul  453!  line  was  split  into  its  two  Stark 
components,  both  of  which  were  shifted  in  the  direction  of  higher  wave¬ 
lengths.  If  Xp  and  X^  represent  the  two  Stark  component  wavelengths 

A 

Other  Stark-effect  experiments  at  slightly  higher  electric  fie1ds  were 
performed  by  Nagaoka  and  Suglura  (1924),  Fujioka  and  Nakamura  (1927), 
but  were  restricted  to  wavelengths  below  about  4100A. 
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and  Xq  the  original  wavelength,  then  assuming  a  linear  Stark  shift  for 
this  transition 


X  +  C  F 
o  p 


X  »  X  +  C  F 
non 


(2-15) 


where  are  the  respective  line  shifts  per  unit  electric  field 

and  F  is  the  electric  field  at  the  emitting  atom. 

The  mean  displacement  can  be  conveniently  expressed  as 

5s  *  (Xp  +  Xn)/2  '  Ko  (2-16) 


or 

6s  *  (Cp  +  Cr )J/2 

The  separation  of  the  Stark  components  is  simply 

d  ■  X  -  X 
p  n 

or 

d  -  (CP  -  Cn)F 

Hence  the  ratio  of  component  separation  to  mean  displacement  Is 
d/6s  -  2(Cp  -  Cn)/(Cp  ♦  C„) 

For  a  quadratic  Stark  shift  of  the  Cul  4531  line  one  has 


(2-17) 


(2-18) 

-/ 

(2-19) 


(S3-  20) 


X  -  X  +  K  F‘ 

pop 


X  »  X  +  K  F 
non 


(2- !5o) 


where  Kp  ,  are  the  respective  line  shifts  per  unit  electric  field 
squared  for  the  p-  and  n-  Stark  components.  The  resultant  ratio  of  com¬ 
ponent  separation  to  mean  displacement  is  of  a  similar  form  to  that  of 
eq  (2-20),  namely 
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A 


d/6.  -  2(K  -  K  ) / (K  +  K  )*  (2- 20a) 

s  p  n  p  n 

In  fitting  an  observed  spectral  profile,  the  Cul  4531  line  was  split 
and  shifted  In  accord  with  the  above  relations.  The  wavelength  positions 
of  the  Cul  4  525  and  4540  transitions  were  unaltered.  A  mu  1 1 (component 
Lorentzlan  profile  match  to  the  observed  profile  at  several  times  was 
made.  The  procedure  used  In  fitting  the  line  profiles  was  as  follows, 

A  single-component  fit  was  first  matched  to  the  principal  Cul  4540  line 
over  a  wavelength  region  that  Included  no  partially  resolved  components. 
This  resulted  in  a  first  estimate  for  the  Cul  4540  halfwidth.  Comparison 
of  the  observed  line  shape  and  this  s ingie -component  fit  indicated  the 
contributions  of  the  partially  resolved  Cul  4525  and  4531  components  to 
the  observed  profile.  Initial  estimates  of  the  relative  intensities  and 
halfwidths  for  these  two  components  were  then  made  from  this  comparison. 

A  multicomponent  fit  to  the  observed  line  profile  was  made  using  these 
initial  estimates.  The  Cul  4531  line  was  shifted  and  split  Into  two  com¬ 
ponents,  each  having  the  same  halfwidth  as  determined  above.  The  relative 

Intensities  of  the  \  -  and  \  lines  were  taken  in  the  ratio  of  two  to 

n  p 

one  since  a  Stark  split  n-component  has  double  the  statistical  weight  of 
a  p-component,  (Greenstein,  I960). 

By  varying  the  resonance-broadened  halfwidth  values  for  all  the 
lines,  the  relative  intensities  of  the  Cul  4525  and  4531  components,  and 
the  line  splitting  and  shift  of  the  Cul  4531  transition,  a  best  fit  to 
the  observed  line  shape  was  achieved.  Resonance  halfwidth  values  and 

*From  Takamine's  data  in  Table  IX  this  ratio  is  for  the  Cul  4531  line 
numerically  d/6s  *  2/7  . 
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relative  intensities  for  these  spectral  lines  used  to  obtain  the  best 
prof  lie  match  were  prev ious ly  g ; ven  in  Table  VII,  part  (b)  and  Table  Vli I 
respectively.  Final  values  of  the  mean  shift  6s  ,  and  component  separa¬ 
tion  d  ,  for  the  Cul  4531  line  used  in  the  profile  fits  are  shown  in 
Table  X.  The  resultant  multicomponent  profiles  at  5.7  and  6.7  p.sec  were 
shown  above  in  Figs.  22  and  23.  For  comparison,  spectral  line  profiles 
In  which  the  Cul  4531  line  was  not  split  or  shifted  In  wavelength  are 
presented  In  Figs.  24  and  25.  The  shift  in  the  observed  Cul  4531  line  at 
each  time  Interval  is  most  obvious  In  these  comparative  profiles.  Thus 
the  observed  wavelength  displacements  can  be  attributed  to  a  Stark  effect 
that  results  from  an  Intermolecular  electric  field  created  by  charged 
particles  in  the  jet  collision  plasma. 


2.4  Neutral  Atom  and  Ion  Density  Estimates 

Estimates  of  the  neutral  copper  atom  density  in  the  jet  collision 
plasma  can  be  obtained  from  the  resonance-broadened  line-halfwidth 
measurements.  The  linear  dependence  of  a  line  halfwidth  on  the  neutral 
particle  density  discussed  previously  will  be  used  in  the  calculations. 
Observed  spectral  lines  whose  upper  energy  level  is  effected  by  autoion¬ 
ization  and  hence  exhibit  a  large  natural  line  width  will  not  be  con¬ 
sidered,  In  addition,  oscillator  strengths  for  two  of  the  spectral  lines 
studied  will  be  obtained  from  their  measured  halfwidths.  This  is  of  par¬ 
ticular  importance  since  experimentally  determined  oscillator  strengths 
for  excited  transitions  of  copper  are  relatively  scarce. 

Observed  Stark  wavelength  shifts  lead  to  estimates  of  the  density 


’■  -  ■■  *  «*  /  v 
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of  copper  Ions  in  the  plasma.  Results  of  Stark  broadening  theory  will  be 
used  to  obtain  the  number  density  of  copper  ions  and  electrons  in  the  jet 
cc! 1  is  ion  plasma. 

(a)  Calculation  of  Neutral  Atom  Density 

Halfwidth  values  of  the  Cul  5153,.  4651,  and  4531  spectral  lines  were 
used  to  calculate  the  neutral  copper  atom  density  in  the  collision  plasma. 
The  relation  between  the  neutral  particle  density  and  a  resonance-broad  - 
ened  line-halfwidth  was  previously  given  in  eq  (2-8),  with  numerical 
values  of  the  constant  in  the  equation  shown  in  Table  VI.  Although  there 
are  small  differences  between  the  values  of  the  constant,  Hargenau  and 
Watson's  statistical  result  was  considered  to  be  appropriate  for  further 
calculations.  This  choice  is  preferable  since  in  evaluating  the  constant 
a  linear  halfwidth-neutral  particle  density  relation  was  shown  to  be  a 
natural  consequence  of  a  statistical  theory  and  requires  no  assumption  of 
an  optical  collision  diameter,  as  in  the  impact  theories.  Furthermore, 
in  view  of  the  high  particle  densities  that  will  be  shown  to  result  from 
the  application  of  eq  (2-8),  frequent  close  encounters  between  emitter 
and  perturber  particles  occur  and  hence  a  statistical  line  theory  is 
appropriate. 

A  resonance-broadened  line-halfwidth  and  neutrai  atom  density  are 

“it 

from  eq  (2-8)  related  as 

A  ■  j  (e2f/mvQ)N  »  j  (e2fKo/mc)N  (2-21) 

For  the  calculations  reported  here  the  consistent  set  of  oscillator 
strengths  given  by  Corliss  and  Bozman  (1962)  for  the  Cul  5153,  4651,  and 
4$31  lines  presented  in  Table  Xi  were  used.  Halfwidth  values  for  both 


Cut  *4651  and  4531  at  5.7,  6.7,  and  7.7  nsec  after  initial  jet  Impact  were 
taken  frcm  the  multicomponent  profile  matches,  part  (b)  of  Table  VII.  At 
later  times  halfwidth  values  from  single  component  fits  were  used  for  the 
number  density  calculations. 

The  resultant  neutral  copper  atom  densities  obtained  for  each  of 

the  three  spectral  lines  a^e  listed  in  Table  XII  and  also  illustrated  in 

Fig.  26.  Error  estimates  included  in  the  figure  indicate  the  sensitivity 

of  the  line  profile  match  in  the  determination  of  a  halfwidth  and  do  not 

reflect  possible  systematic  errors  in  the  measurements  and  data  reduction 

techniques.  Variations  due  to  the  fitting  sensitivity  are  of  the  order 

of  -15  percent.  Further  discussion  of  the  uncertainties  in  these  results 

Is  presented  below  in  section  2.5.  The  overall  agreement  In  the  density 

estimates  from  these  three  spectral  lines  is  good.  Results  derived  from 

the  Cul  5153  and  4651  halfwidths  are  especially  close  throughout  most  of 

the  time  interval  studied.  The  neutral  atom  density  in  the  jet  collision 

19  -3 

plasma  is  seen  to  be  of  the  order  of  4  x  10  c:n  , 

The  structure  and  radiative  transport  processes  in  the  expanding 
Jet  collision  sheet  have  been  studied  by  Lukas ik,  et  al  (1964).  The  time 
variation  of  material  density  along  the  radius  of  the  expanding  sheet  has 
been  calculated  from  a  fluid-dynamic  model  of  the  jet  collision.  Typical 


For  A  in  angstrom  units  one  has  numerically 

A  «  (4.68  I0"6)  fx’  N  (2-21a) 

where  A(angular  frequency  units)  ■  2rr  A  (sec  *)  and 
A(angstrom)  »  10®  A(cm)  ■  \^/c  *  10®  x  A  (sec  ^) 
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profiles  show  a  moderate  increase  In  material  density  in  a  radially  out¬ 
ward  direction  from  the  point  of  jet  impact,  i.e.  the  center  of  the 
expanding  sheet.  After  reaching  a  maximum  value  at  about  50  percent  of 
the  radial  distance  from  the  center  to  the  front  of  the  sheet,  the  density 
profile  fails  rapidly  to  zero  at  the  front.  Numerical  solutions  of  the 
radiative  transport  equations  using  these  density  profiles  and  calculated 
optical  opacity  values  were  performed  by  Lukas ik,  et  al.  The  results  of 
these  calculations  indicate  that  the  observed  spectra!  radiation  comes 
from  a  spatially  narrow  region  between  the  expanding  front  of  the  colli¬ 
sion  sheet  and  the  radial  position  of  peak  density.  Thus  the  neutral 
copper  atom  densities  measured  from  resonance-broadened  haifwidths  are 
indicative  of  conditions  in  this  region. 

A  variation  of  the  neutral  atom  density  with  time  is  apparent  in 
Fig.  2 6.  The  particle  density  initially  decreases  until  about  9  psec 
after  jet  impact  and  then  increases.  This  feature  of  the  spectral  record 
is  in  accord  with  results  of  the  fluid-dynamic  model .  In  particular  one 
can  show  that  the  front  of  the  expanding  sheet  moves  with  a  constant 
velocity.  The  calculated  time  of  arrival  of  the  expanding  sheet  at  the 
vacuum  chamber  wall  is  8.8  psec  from  jet  impact.  This  is  in  excellent 
agreement  with  the  observed  time  for  a  minimum  atom  density  of  about 
9  psec.  What  is  thought  to  occur  physically  is  that  after  the  expanding 
sheet  strikes  the  chamber  wall  a  compressive  shock  is  reflected  radially 
inward  along  the  sheet.  This  inward-moving  shock  compresses  the  material 
in  the  sheet.  The  observed  increase  in  neutral  atom  density  after  9  psec 
is  in  agreement  with  this  qualitative  prediction.  Furthermore,  a  wall 
collision  will  result  in  radiative  recombinations  of  the  excited  jet 


materia!  and  heating  of  the  chamber  wall,  which  would  correspond  to  the 
onset  of  a  background  continuum  during  the  relight  phase.  At  the  same 
time  the  Increase  in  atom  density  leads  to  the  copper  absorption  lines  in 
the  time-resolved  spectrum  (Fig.  14). 

Since  the  observed  halfwidtns  are  governed  by  eg  (2-21)  coupling 
halfwidth  measurements  of  other  lines  with  the  now  known  number  densities 
can  be  used  to  determine  transition  probabilities.  Oscillator  strengths 
for  the  Cul  4643  and  4525  transitions  have  been  determined  in  this  manner. 
At  a  given  time  after  jet  Impact  the  ratio  of  halfwidths  for  two  resonance 


broadened  lines  X^  ,  Xj  can  be  expressed  as,  (see  eq  (2-21a)) 
A./Aj  *  f j/f j  x  (X./X.)3 


(2-22) 


If  the  oscillator  strength  of  one  of  the  lines  Is  known  then  the  oscilla¬ 
tor  strength  of  the  second  line  can  be  obtained  from  this  halfwidth  ratio. 
The  Cul  4651  transition  was  chosen  as  a  reference  line  since  its  oscilla¬ 
tor  strength  is  known  and  its  wavelength  is  close  to  the  Cul  4643  and 
4525  lines.  Comparison  line  halfwidths  were  taken  from  part  (b)  of  Table 
VII;  the  Cul  4651  oscillator  strength  was  given  in  Table  XI.  Halfwidth 
ratios  ^4543^4551  and  *4525^4651  ’  were  calculated  at  the  three  times  for 
which  multicomponent  profile  fits  were  made.  The  oscillator  strengths 
obtained  for  these  lines  are  given  in  Table  XIII.  The  three  values  Indi¬ 
cated  for  each  line  are  derived  from  the  measured  halfwidths  at  the  three 
early  times.  A  weighted  mean  oscillator  strength  for  each  spectral  line 
Is  also  shown  in  the  table.  The  relative  intensities  of  the  observed 
spectral  lines  given  In  Table  VIII,  were  used  to  weight  the  Individual 
f-number  values  in  determining  the  mean.  This  reflects  the  sensitivity 
of  the  fitting  process  in  e  determination  of  parameters  that  define  a 
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resonance-broadened  Hne  profile, 

(b)  Calculation  of  Ion  Density 

The  observed  wavelength  displacements  of  the  Cul  A53I  transition 
are  due  to  a  Stark  splitting  and  shift  of  this  line.,  a  consequence  of  the 
intermolecular  electric  field  that  exists  in  the  jet  collision  plasma. 

The  resultant  line  shifts  are  a  measure  of  this  field  strength  which  in 
turn  Is  established  by  the  density  of  charged  particles  in  the  plasma. 
Copper  Ions  and  hence  electrons  are  known  to  exist  in  the  jet  collision 
products  as  evidenced  by  the  appearance  of  character i st lc  ion  lines  in 
the  time-resolved  spectrum.  Contributions  to  the  intermolecular  field 
from  the  low  pressure  helium  gas  in  the  vacuum  chamber  which  may  have 
been  partially  ioniied  by  collisions  with  the  expanding  jet  plasma  are 
extremely  small  in  view  of  the  low  helium  density  and  the  high  copper  ion 
density  found  to  exist  in  the  plasma. 

Consider  an  emitting  (neutral)  atom  In  the  presence  of  a  number  of 
perturbing  ions  and  electrons.  The  net  electric  field  at  the  position  of 
the  emitter  Is  the  sum  of  individual  contributions  from  all  particles. 

For  every  charged  particle  configuration  there  will  be  a  corresponding 
net  electric  field  strength.  By  virtue  of  the  motion  of  these  charges 
their  positions  will  change  with  time  and  so  also  will  the  value  of  the 
electric  field  strength.  As  a  result  there  is  a  probability  distribution 
for  the  field  at  the  emitter;  Stark  splitting  and  shifts  will  be  governed 
by  this  distribution.  Each  Stark  component  of  a  spectral  line  will  be 
broadened  and  displaced  In  accord  with  the  local  field  distribution.  In 
the  absence  of  other  broadening  mechanisms  tfe  sum  of  the  contributions 
from  each  Stark  component  will  form  the  total  line  profile.  If  the  Stark 


49 


component  lines  are  symmetrically  displaced  about  the  center  wavelength, 

(as  Is  the  case  for  hydrogen),  then  the  line  profile  will  only  be  broad¬ 
ened  by  the  presence  of  charged  particles.  On  the  other  hand  if  the 
pattern  of  the  Stark  splitting  is  asymmetric  about  the  center  wavelength 
the  total  line  shape  will  indicate  in  addition  a  wavelength  displacement 
determined  by  the  form  of  the  electric  field  distribution. 

Holtsmark  Theory-Linear  Stark  Effect 

Holtsmark  (1919*  1924)  first  derived  a  probability  distribution  from 

a  statistical  theory  for  the  electric  field  strength  at  the  emitter  due  to 

* 

the  presence  of  nearby  ions.  If  F  represents  the  magnitude  of  the 
field  and  Fq  an  average  field  strength  given  by 

Fq  -  2.61  e  Nj273  (2-2?) 

where  e  is  the  electronic  charge  in  cgs  units  and  Nj  the  number  dens¬ 
ity  of  ions  then  the  Holtsmark  distribution,  usually  represented  as  W(0), 
can  be  expressed  as  a  function  of  the  dimensionless  ratio, 

fl  *  F/Fo  (2-24) 

Tables,  graphs  and  asymptotic  representations  of  the  function  W(0)  can 

be  found  in  the  literature  (Margenau  and  Lewis,  1959;  Greenstein,  I960; 

Breene,  1961) .  Under  the  influence  of  this  electric  field  distribution 

each  Stark  component  of  a  spectral  transition  will  be  shifted  to  various 

frequencies  In  accord  with  the  probability  function  W(3)  .  If  Avq  is 

the  shift  due  to  the  average  field  F  then  for  a  linear  Stark  effect 

o 

For  most  physically  realizable  plasma  densities  and  temperatures,  includ¬ 
ing  that  plasma  reported  here,  electron  effects  cannot  be  treated  by  a 
statistical  theory  but  are  described  by  an  impact  theory,  (Hargenau,  (1999)). 


where 


one  has 


Av  -  C.  F 

i 


Av/Avq  »  F/Fo  “  3 


(2-25) 


(2-26) 


Av  is  the  frequency  difference  from  the  normal  position  of  the  line  and 
C.  is  the  Stark  displacement  per  unit  electric  field  for  the  ith  Stark 
component  of  a  line.  The  spectral  line  intensity  distribution  J(Av)  , 
is  given  in  terms  of  the  Holtsmark  function  as  (Henkel,  1954;  Griem, 
et  al,  1959;  Penner,  1959;  Greenstein,  I960) 

S. 


J(AV)  -  S  c7  W(^-) 
i  <  I 


(2-27) 


where  S}  Is  the  relative  intensity  of  the  ith  Stark  component.  The  peak 
of  each  component  will  be  displaced  by  an  amount 

where  8  is  the  value  of  0  for  which  W(p)  is  a  maximum.  Thus  for 

IftaX 

a  spectral  line  with  only  two  Stark  components  the  mean  line  shift  is 

S5  *  1  {(AV)P  +  (AV)„} 


or 


ft  -  -5  (C  ♦  C  )  F  B 
s  2  p  n  o  max 


(2-29) 


A  value  for  the  average  electric  field  can  be  obtained  from  a 
measured  shift  6s  .  The  ion  density  corresponding  to  this  field  strength 
can  then  be  calculated  from  eq  (2-23).  Copper  ion  densities  in  the  jet 
collision  plasma  were  computed  in  this  manner  from  the  measured  line  shifts 
of  the  Cul  4531  line  given  In  Table  X.  The  results  shown  in  Table  XIV 


A  value  of  8 


max 


of  3.5  was  used  for  the  calculations  (Breene,  1957) 


«?  • 
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Indicate  that  ion  densities  of  the  order  of  5  x  10  cm  *  are  obtained 
from  such  a  treatment  of  linear  Stark  d isplaeements  using  Holtsmark's 
theory. 

Folded  Distribution-Linear  Stark  Effect 

In  certain  situations  the  Holtsmark  theory  has  proven  to  be  quite 

adequate  In  the  description  of  the  Influence  of  ionic  fields  on  spectral 

lines  (Edels  and  Craggs,  1951;  Olsen  and  Huxford,  1952;  Lochte-Hol tgreven, 

1958).  Recent  Investigations  have  shown  that  a  reliance  on  the  Holtsmark 

theory  is  possible  for  low  charged  particle  densities,  of  the  order  of 
1 6  m  3 

10  cm”  or  less;  however  for  higher  densities  electron  effects  should 
be  Included  Into  a  more  comprehensive  theory.  Combined  Ion  and  electron 
effects  at  high  densities  on  the  Stark  broadening  of  hydrogen  lines  have 
recently  been  studied  by  Greim  (1954),  Kolb  (1956),  Baranger  (1958),  Kolb 
and  Griem  (1958),  Griem,  Kolb  and  Shen  (1959).  Comparative  experiments 
performed  by  Henkel  (1954),  Turner  (1956),  and  Bogen  (1957)  have  shown 
that  the  inclusion  of  electronic  contributions  to  hydrogen  line  profiles 
result  in  a  substantially  improved  fit  than  just  the  Holtsmark  profile 
alone  and  that  the  added  electron  broadening  is  comparable  In  magnitude 
to  the  ionic  broadening  In  the  wings  of  a  spectral  line.  Their  Stark 
broadened  hydrogen  lines  were  found  to  be  well  represented  by  an  expres¬ 
sion  of  the  form 

CD 

J(v)  -  const  x  j  W(F )  J (v,F)  dF  (2-30) 

vo 

where  W(F)  is  the  Holtsmark  distribution  function  and  J(v,F)  is  the 
electron  contribution  to  the  line  broadening.  From  impact  theory  the 
electron  broadened  profile  is  given  by  (Hargenau  and  Lewis,  1959), 


(2-3  0 


J(v,F)  *  const/((v-vo  -  CF)c  +  y2) 

where  y  )s  the  halfwidth  of  the  eiectron-broadened  profile  and  C  is 
the  Stark  displacement  per  unit  field  for  a  Stark  component  of  a  line, 

A  complete  line  profile  is  generated  by  evaluating  the  folded  distribu¬ 
tion,  eq  (2-30)  for  all  Stark  components. 

Direct  evidence  for  spectral  line  shifts  due  to  the  presence  of  a 
high  election  density  Is  not  apparent  in  their  experiments  since  the 
Mnear  Stark  splitting  for  hydrogen  is  syrsnetric  about  the  central  wave¬ 
length.  Line  shifts  In  helium  and  argon  which  exhibit  an  asymmetric  quad¬ 
ratic  Stark  pattern  have  been  experimentally  studied  by  Seay,  et  al  (196!) 
and  Petschek,  et  al  (1955).  Seay,  et  a!  investigated  several  helium  lines 
and  compared  their  experimentally-determined  line  profiles  and  shifts  with 
computed  profiles  and  shifts,  using  quantum  mechanical  electron  halfwidth 
and  shift  calculations  of  Kivel  (1955).  The  uncertainties  in  the  neces¬ 
sary  matrix  elements  for  argon  made  theoretical  comparisons  with  the 
measurements  of  Petschek,  et  al  somewhat  difficult.  Nevertheless,  it  was 
shown  in  both  cases  that  the  observed  line  shifts  were  greater  than  those 
predicated  from  Holtsmark's  theory  alone  and  indicated  that  high  electron 
densities  can  also  contribute  to  a  spectra!  line  shift. 

Estimates  of  line  shifts  in  the  time-resolved  spectrum  due  to  the 
simultaneous  action  of  copper  ions  and  electrons  present  in  the  jet  colli¬ 
sion  plasma  were  also  made.  The  calculations  consisted  of  a  determination 
of  the  combined  electron-ion  broadened  line  profile  for  both  Stark  split 
components  of  the  Cul  4531  transition,  by  application  of  eqs  (2-30)  and 
(2-31).  For  a  given  cl.jrged  particle  density*  the  shift  of  the  peak 

Equal  Ion  and  electron  densities  were  assumed. 


position  of  the  line  were  found,  Kolb's  expression  for  the  electron 
broadened  halfwidth  y  was  used  (see  Appendix  G) .  The  variation  of  Cul 
4531  to an  line  shift  6g  ,  with  particle  density  found  in  this  manner  is 
given  in  Fig.  27.  For  comparison,  line  shifts  predicted  by  Holtsmark*s 
theory  (eq  (2-29)),  are  also  shown.  Thus  the  measured  line  shifts  given 
In  Table  X  lead  to  the  revised  mane r i ca I  estimates  for  the  ion  end  elec¬ 
tron  density  shown  in  Table  XIV.  Chargeu  particle  densities  in  the  jet 
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collision  plasma  of  the  order  of  4  x  10  cm  3  are  Indicated  and  are  about 
15  percent  smaller  than  the  corresponding  density  estimates  found  using 
Holtsmark's  theory  alone  as  was  illustrated  in  Fig.  27.  Descriptions  of 
the  Fortran  programs  employed  in  these  computations  are  given  in  Appen¬ 
dix  G. 

Holtsmark  Theory  -Quadratic  Stark  Effect 

The  nature  of  the  dependence  of  the  Cul  4S31  line  shift  with  applied 
electric  field  bears  further  investigation.  The  displacement  of  an  energy 
level  under  the  influence  of  an  electric  field  depends  upon  the  proximity 
of  other  nearby  perturbing  energy  levels  with  opposite  parity  to  that  of 
the  emitting  level  (Condon  and  Short  ley,  1959;  Minnhagen,  1964).  If  these 
perturbing  levels  are  energetically  close  to  the  excited  level  of  concern 
then  the  spectral  transition  originating  from  this  level  will  exhibit  a 
linear  Stark  shift  with  applied  field.  If  on  the  other  hand  the  excited 
level  Is  far  removed  from  other  perturbing  levels,  then  higher  order  line 

*The  Stark  displacements  of  lower  energy  states  are  in  general  much  smaller 
than  those  of  the  upper  states  so  that  the  observed  shifts  are  for  the 
most  part  indicative  of  the  Stark  displacement  of  the  upper  level  of  a 
trans It  Ion. 


shifts  are  prevalent.  For  example,  hydrogen  spectra*  Sines  exhibit  large 
linear  Stark  displacements  at  low  field  strengths  4 ye  to  the  degeneracy 
of  excited  levels  In  the  absence  of  an  electric  field  As  the  field 
strength  Is  increased  this  degeneracy  is  removed.  I.e.,  the  initially 
overlapping  levels  are  further  separated.  As  a  result  of  this  Increased 
level  separation  a  quadratic  variation  of  the  line  shift  with  electric 
field  becomes  important  (Kolb,  1956).  Alternately,  the  NaD-lines  show 
only  a  quadratic  Stark  effect  since  the  upper  levels  for  these  transi- 
*■  ons  are  far  removed  from  other  perturbing  states  (White,  193*0. 

Cul  energy  level  tabulations  (Moore,  1952)  indicate  that  the  near- 
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perturbing  states  of  opposite  parity  to  the  6s  3  upper  level  for  the 
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Cut  4531  line  are  the  6p  P  states,  some  2000  cm  removed.  Whether  this 
energy  separation  is  sufficient  to  Indicate  a  clear  choice  between  a 
linear  or  higher  order  Stark  effect  for  the  Cul  4$3l  line  is  not  certain. 
However,  the  following  considerations  strongly  suggest  that  the  Cul  **531 
Une  exhibits  a  quadratic  Stark  effect. 

Hydrogen  term  separations  for  which  a  linear  Stark  effect  predomi¬ 
nates  are  extremely  small,  of  the  order  of  0.05  cm  V  The  7d^D  levels  of 
helium  interacting  with  nearby  levels  at  a  separation  of  about  1 50  cm"  * 
produce  a  first  order  Stark  effect  in  fields  of  about  100  kv/cm.  Some¬ 
what  closer  to  the  copper  level  spacing  under  consideration  are  the 
Interact Ing  Ss8*,  5p2P°,  and  5d^D  states  associated  with  the  principal 
series  of  potassium.  Typical  energy  separations  are  of  “he  order  of 
4000-5000  cm'*.  These  levels  show  quadratic  Stark  displacements  in  fields 

of  the  order  of  100  kv/cm.  At  the  uppe-  extreme  perturbing  level  spacings 
P  o 

about  the  3p  r  states  for  the  NaD-lines  which  show  a  quadratic  Stark 


effect  are  in  excess  of  50,000  cm 

Thus,  these  comparisons  indicate  that  the  Cui  ^53  I  tine  exhibits  a 
quadratic  Stark  shift  with  electric  field.  Although  both  Holtsmark  and 
Trumpy  (1925),  and  Margenau  and  Watson  (1936)  have  considered  Takamlne's 
Stark  effect  measurements  to  Imply  a  linear  variation  with  applied  elec¬ 
tric  field,  Fujioka  and  Nakamura's  measurements  for  several  spectral  lines 
of  silver  show  a  mixture  of  linear  and  non-linear  Stark  component  dis¬ 
placements  under  experimental  conditions  similar  to  those  for  which  the 
copper  line  shifts  were  measured. 

Previous  Ion  density  estimates  must  be  rev!sed  to  reflect  the  quad¬ 
ratic  nature  of  the  Cui  -453 f  Stark  splitting  and  shift.  in  analogy  to 
eq  (2-29)  the  mean  line  shift  is  now  given  by 


5  «  ~  (K  +  K  )  (F  3  )2 

s  2  p  n'  v  o  Mnax' 


(2-32) 


where  Kp,  Kn  are  the  line  shifts  per  unit  electric  field  squared  for 

the  p-  and  n-  Stark  components  respectively.  Copper  Ion  densities  in  the 

jet  collision  plasma  were  computed  from  the  measured  line  shifts  of  the 

Cut  453 1  line  given  in  Table  X  and  from  the  relation  between  the  average 

intermolecular  field  strength  and  the  Ion  density  by  means  of  eqs  (2-32) 
and  (2-23)  respectively.  The  results  given  in  Table  XIV  Indicate  that 
Ion  densities  of  the  order  of  4  x  10  cm  are  present  In  the  jet  colli¬ 
sion  plasma.  These  estimates  are  lower  by  a  factor  of  ten  from  previous 
results.  The  variation  of  the  Cui  453 1  mean  line  shift  ,  with  ion 
density  found  in  this  manner  is  also  given  in  Fig.  27.  Estimates  of  the 
electron  contribution  to  the  observed  line  shift,  obtained  from  an  eval¬ 
uation  of  the  folded  distribution  for  a  quadratic  Stark  effect,  were  not 
made  since  the  electronic  contributions  were  previously  found  to  be  small. 
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2,5  Error  and  Uncertainty  Estimates 

A  number  of  systematic  errors  in  the  measurements  and  data  reductton 
procedures  are  inherent  in  the  above  results.  The  extent  to  which  these 
errors  influence  the  overall  accuracy  of  these  resuits  is  discussed  here. 
No  attempt  has  been  made  to  ascertain  limits  of  error  in  each  detailed 
measurement  and  calculation",  instead  uncertainty  estimates  arrived  at  by 
studying  typical  situations  were  considered  to  be  applicable  to  all  other 
similar  measurements. 

(a)  Film  Sens i t  i  v  i  ty 

The  halfwidth  of  a  spectral  line  as  determined  from  a  line  profile 
will  depend  upon  the  value  of  the  slope  of  the  photometric  emulsion  cali¬ 
bration  v  used  to  obtain  the  line  shape.  A  variation  In  y  will  In 
general  alter  this  ^roflle  and  hence  the  observed  halfwidth  of  the  spec¬ 
tral  'ine.  Photometric  calibration  studies  summarized  in  Table  V  indicate 
that  the  uncertainty  in  y  i s  of  the  order  of  0.1  density  unit  about  a 
mean  value  of  unity.  To  ascertain  the  sensitivity  of  a  measured  line 
halfwidth  to  this  y-uncertai nty  three  profiles  of  the  tut  5153  line  at 
8.7  tisec  were  constructed  using  values  for  y  of  0.9,  1.0,  and  1.1. 

The  total  observed  halfwidth  of  this  line  is  about  14A.  The  maximum  devi¬ 
ation  In  the  total  halfwidth  were  found  to  be  of  the  order  of  0.3A  or 
about  2  percent  of  the  total  halfwidth.  This  estimate  is  considered  to 
be  representative  of  the  error  in  all  halfwidth  measurements  since  the 
line  profiles  studied  extended  over  comparable  ranges  in  film  density. 
Thus,  the  uncertainty  In  a  total  halfwidth  determination  is  of  the  order 
of  -2  percent  over  the  total  range  of  observed  line  halfwldths. 
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(b)  Wavelength  Calibration  Effect  on  Line  Width 

Errors  In  the  wavelength  scale  calibration  of  a  densitometer  trace 
will  also  distort  the  observed  profile  and  affect  the  accuracy  of  the 
halfwldth  measurement.  The  transformation  from  a  horizontal  position 
measurement  to  wavelength  on  each  densitometer  record  was  attained  with  a 
dispersion  curve  of  the  form 

v  *  a  +  bx  +cx^  (2-33) 

where  v  is  the  wavenumber  of  reciprocal  wavelength  corresponding  to  the 
position  x  ,  The  coefficients  a,b,c  were  determined  from  a  least  square 
fit  of  eq  (2-33)  using  the  measured  positions  of  known  spectral  lines  on 
each  densitometer  trace,  (see  Appendix  C) .  The  error  In  a  given  wave- 
number  measurement  can  be  expressed  in  terms  of  the  uncertainties  in  the 
numerical  coefficients  and  the  error  In  a  position  measurement,  (Beers, 
1958).  Position  errors  were  negligibly  small  as  evidenced  by  the  resolu¬ 
tion  of  the  densitometer  record  scanning  device  (Gerber  model  S-10-C)  and 
the  repeatability  of  successive  position  measurements.  It  was  found  that 
the  quadratic  term  in  eq  (2-33)  was  comparatively  small  and  could  also  be 
neglected  in  estimating  errors.  From  the  approximate  relation 

v  **  a  +  bx  (2-31*) 

the  contribution  to  the  overall  accuracy  in  wavelength  due  to  uncertain¬ 
ties  in  the  coefficient  c  were  of  the  order  of  IA.  Error  contributions 
due  to  uncertainties  In  b  were  much  smaller,  of  the  order  of  0.05  to 
0.10A  for  the  range  of  observed  line  widths.  This  latter  estimate  is  of 
Importance  here  since  wavelength  differences  from  a  center  wavelength 
value  were  used  for  all  data  points  describing  a  line  profile.  Thus  the 
distortions  of  an  observed  line  profile  In  the  wavelength  direction  and 
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the  consequent  uncertainties  in  the  total  line  halfwidth  are  of  the  order 
of  0.05  to  0.10A.  This  amounts  to  about  l  percent  of  a  typical  line  half¬ 
width. 

(c)  Spectrograph  SHt  Width 

For  a  given  line  shape  errors  In  the  slit  width  value  at  the  wave¬ 
length  of  the  observed  line  contribute  to  the  error  In  a  measured  reson¬ 
ance  halfwldth  as  determined  from  a  Lorentzlan  profile  fit.  Fits  to  the 
sample  Cul  5153  line  were  made  In  which  ten  percent  variations  were  intro¬ 
duced  for  both  the  assigned  slit  width  and  resonance  halfwldth  values. 

All  three  profiles  appeared  to  adequately  represent  the  observed  line 
shape.  Larger  variations  resulted  In  a  degradation  of  the  quality  of  the 
profile  fit.  Experience  In  the  determination  of  a  number  of  single  com¬ 
ponent  Lorentzlan  line  profiles  have  indicated  In  general  that  the  con¬ 
struction  of  a  given  profile  with  approximately  ten  percent  variations  in 
either  the  slit  width  or  resonance  halfwidths  can  result  In  an  acceptable 
profile  match.  For  multicomponent  Lorentzlan  line  profiles  small  changes 
of  the  order  of  ten  to  fifteen  percent  in  the  slit  width,  resonance  half¬ 
widths,  and  relative  intensities  of  the  spectral  component  were  also  found 
to  yield  acceptable  line  fits. 

(d)  Densitometer  Slit  Width 

The  finite  width  of  the  densitometer  entrance  slit  will  modify  the 
shape  and  hence  the  breadth  of  a  spectral  line.  This  comes  about  due  to 
an  averaging  of  the  film  density  across  the  densitometer  entrance  slit. 
Estimates  of  the  densitometer  slit  broadening  given  In  Appendix  H  indicate 
that  the  error  In  the  calculated  profile  # (\>)  ,  depends  upon  the  square 
of  the  dimensionless  ratio  of  the  densitometer  silt  width  (in  wavelength 
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units)  to  resonance  halfwidth  A  .  For  the  slit  width  used  in  obtaining 
the  densitometer  records  and  typical  resonance  halfwtdths  as  given  In 
Table  VII  the  maximum  error  In  the  calculated  profile  Is  less  than  one 
percent . 

(e)  Absolute  Oscillator  Strength 

The  accuracy  of  the  neutral  atom  density  estimates  Is  in  part  given 
by  the  above  variations  in  resonance  halfwidths.  Errors  In  the  oscillator 
strength  values  used  in  the  density  calculations  can  also  affort  signifi¬ 
cant  contributions.  Corliss  and  Bozman  (1962)  have  assessed  the  accuracy 
of  their  oscillator  strength  determinations  and  have  found  that  their 
absolute  f-values  can  be  In  error  by  as  much  as  a  factor  of  two  due  to 
uncertainties  In  their  measurements  and  the  associated  data  reduction. 

This  Is  by  far  the  largest  source  of  error  in  the  neutral  atom  number 
density  evaluations. 

(f)  Self-Absorption 

In  view  of  the  relatively  high  particle  densities  found  to  exist  in 
the  Jet  collision  plasma,  the  effects  of  self-absorption  on  the  collision 
luminosity  must  be  considered.  The  intensity  of  the  spectral  line  radia¬ 
tion  from  Interior  portions  of  the  plasma  will,  in  general,  be  diminished 
upon  passage  through  the  absorbing  exterior  regions.  This  will  result  in 
a  distortion  of  the  observed  line  profile  and  consequently  an  error  in  a 
measured  line  halfwidth.  The  overall  agreement  in  the  neutral  atom  number 
density  estimates  for  the  three  lines  studied  offers  indirect  evidence 
that  absorption  effects  have  not  significantly  altered  the  observed  line 
shapes.  In  particular,  the  close  agreement  In  particle  density  values 
derived  from  halfwidth  measurements  of  the  Cul  4651  and  5153  lines  is 
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encouraging  since  these  lines  represent  transitions  with  lower  energy 
levels  that  are  about  2ev  apart  and  have  oscillator  strengths  that  differ 
by  a  factor  of  two.  Thus,  these  lines  could  be  expected  to  exhibit  dif¬ 
ferences  due  to  varying  amounts  of  seif-absorption.  Nevertheless,  it  is 
of  importance  to  further  investigate  the  extent  to  which  so! f -absorpt Ion 
could  modify  the  observed  spectral  line  radiation,  A  decrease  in  intensity 
of  a  spectra!  profile  about  the  line  center  will  lead  to  an  error  in  a 
line  halfwidth  measurement.  This  will  result  In  an  overestimate  of  the 
true  line  halfwidth  and  consequently  of  the  neutral  particle  density  value 
derived  from  this  estimate. 

In  order  to  study  the  effects  or  se 1 f-absorpt ion  consider  an  array 
of  plane-parallel  absorbing  regions  which  are  traversed  by  radiation  with 
an  initially-known  profile  JQ(v)  .  Let  the  number  of  absorbing  particles 
and  temperature  of  each  region  differ  from  that  of  adjacent  regions.  After 
passage  through  the  first  absorbing  section,  the  profile  of  the  emerging 
radiation  will  be  given  by  (Chandrasekhar,  1939;  Aller,  1953) 

-T,(v) 

J j (v)  =  JQ(v)  e  (2-35) 

where  Tj(v)  Is  the  optical  depth  and  Jj(v)  the  resultant,  intensity. 

The  optical  depth  Is  frequency-dependent,  I.e.,  the  spectral  line  Is 
selectively  absorbed  This  Intensity  profile  will  in  turn  be  further 
altered  upon  passage  through  the  second  absorbing  section  as 

-Tp(v) 

J2(v)  -  J,(v)  e 
In  general,  for  the  j th  slab 

-Tj  (V) 


Jj(v)  -  Jj_,(v)  e 


(2-36) 


and  hence  for  the  radiation  exiting  from  the  regions  of  absorption  after 
traversing  a  distance  R  from  the  entrance  face 
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JR(v)  *  JQ(v)  • 


-V"> 


(2-37) 


The  optical  depth  of  the  jth  region,  Xj (v)  ,  Is  known  In  terms  of 
the  atomic  absorption  coefficient  Oj (v)  and  the  number  density  of  absorb¬ 
ing  atoms  N  ,  as 
a 


Xj  (v) 


(v)  dr 


(2-38) 


J-l 


The  atomic  absorption  coefficient  Is  given  os  (Penner,  1959;  Kulper,  1962; 
Allen,  1963) 


of(v)  -  -*5 


abs 


mc  brr2  (v-v  )2  +  (y/krr)2 


(2-39) 


where  y  Is  the  collision  damping  constant  or  halfwldth  due  to  broadening 

•fa 

collisions,  the  absorption  oscillator  strength,  vQ  the  center 

frequency  of  the  line,  and  e,m,c  have  th  1 1  usual  meanings. 

Self-absorption  effects  on  spectral  line  radiation  from  the  interior 
of  the  jet  collision  plasma  upon  traversal  through  the  outer  fringes  of  the 
plasma  were  derived  from  this  one-d 1 mens  Iona  1  representation.  Particle 
density  and  temperature  profiles  obtained  from  the  fluid-dynamic  model 


The  magnitude  of  the  emission  and  absorption  oscillator  strengths  are 
related  by 


9u  fem  "  9l  *abs 

where  gy  ,  are  the  statistical  weights  of  the  upper  and  lower  states 
of  the  transition. 
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of  a  jet  collision  were  used  In  calculating  absorption  coefficients  and 
optical  depths.  Previous  results  (section  2.4}  have  indicated  that  the 
observed  radiation  originates  from  a  spatially  narrow  region  between  the 
expanding  front  of  the  collision  sheet  and  the  position  of  peak  density. 
THet  is,  the  observed  spectral  line  profile  is  fully  developed  in  thrs 
narrow  region.  Thus  the  radiation  incident  upon  the  absorbing  layers  was 
assumed  to  have  a  resonance-broadened  profile  whose  haifwidth  was  deter¬ 
mined  from  the  measured  neutral  atom  density.  The  entrance  face  of  the 
exterior  absorbing  layers  was  at  the  radial  position  at  which  the  measured 
atom  density  was  equal  to  the  calculated  density.  This  model  appeared 
reasonable  since  the  absorbing  layers  under  consideration  corresponded  to 
a  region  of  the  calculated  profiles  where  both  the  particle  density  and 
temperatures  were  rapidly  decreasing  as  one  progressed  towards  the  exit 
face  or  front  of  the  plasma.  Furthermore,  the  observed  profile  would  be 
a  close  reproduction  of  the  incident  profile  if  self-absorption  effects 
are  small,  as  Is  indicated  by  the  above  mentioned  agreement  in  neutral 
atom  density  estimates  derived  from  Cut  4651  and  5153  halfwidth  measure¬ 
ments. 

Numerical  calculations  for  the  optical  depth  and  emerging  spectral 
line  profile  determined  by  eqs  (2-36),  (2-38),  and  (2-39)  were  performed. 
Details  of  the  computer  program  used  and  typical  results  are  contained  in 
Appendix  I.  The  computed  spectral  profiles  of  the  Cul  4651  and  5153 
transitions  for  two  values  of  the  time  after  impact  are  shown  in  the 

•ft 

This  assumption  will  tend  to  overestimate  the  amount  of  absorption  since 
one  neglects  the  decrease  in  the  number  of  absorbing  particles  due  to 
Ionization.  The  effect  is  small,  however,  for  the  degree  of  ionization 
found  to  exist  in  the  collision  plasma. 


figure*  of  this  appendix.  The  magnitude  of  the  temperature  of  the 
material  at  the  entrance  face  of  the  absorbing  layers  appeared  at  a 
parameter  In  the  calculation  and  is  indicated  In  the  figures.  Each  line 
profile  was  normalized  to  the  maximum  intensity  value  appearing  as  a 
rasult  of  the  computation  to  facilitate  comparison  with  previously  dis¬ 
played  line  profiles. 

For  the  Cul  4651  line  absorption  effects  are  negligible  below  a 
temperature  of  about  4CCC°A;  the  line  center  has  been  diminished  by  at 
most  a  few  percent.  In  the  temperature  range  between  4000  and  5000°A 
absorption  becomes  pronounced  as  indicated  by  the  Increased  width  of  the 
line  and  the  appearance  of  a  characteristic  dip  or  self-reversal  at  the 
line  core.  Above  5000°A  the  spectra!  line  Is  strongly  attenuated  and 
distorted. 

Similarly  for  the  Cul  5153  transition,  profile  distortions  due  to 
self-absorpt Ion  are  small  for  temperatures  below  3000°A.  Between  3000 
and  4000°  absorption  effects  are  noticeable,  again  as  Indicated  by  the 
Increased  width  of  the  line  and  self-reversal  at  the  line  center.  Above 
4000°A  the  Cul  5153  line  is  strongly  attenuated  and  distorted. 

Estimates  of  the  degree  of  self-absorpt I on  In  the  observed  spectre) 
lines  cen  be  obtained  by  a  comparison  of  the  observed  line  shapes  and  tha 
computed  profiles  In  the  vicinity  of  the  line  center,  under  the  assump¬ 
tion  that  obsarvad  distortions  from  the  fitted  curves  are  attributable  to 
absorption  affects.  A  few  of  the  Cul  4651  line  shapes  In  Fig.  19  exhibit 
a  decrease  In  Intensity  at  the  line  core,  •  which  the  most  pronounced  dip 
appears  with  tha  6.7  psee  result  (saa  elsu  g.  21).  Tha  two  side  peaks 
separated  by  about  IA  for  this  profile.  Keeping  In  mind  that  tha 


spectrograph  entrance  slit  width  has  essentially  smoothed  variations  In 
the  observed  line  profiles,  but  that  side  peak  separations  are  relatively 
unchanged  by  entrance  modifications,  comparison  with  appropriate  computed 
absorption  profiles  suggests  that  errors  In  a  halfwidth  determination  may 
be  of  the  order  of  40  to  50  percent  in  this  case.  For  the  Cul  5153  line 
profiles  In  Fig,  SO  similar  absorption  effects  are  seen  to  occur.  In 
particular  the  line  profile  at  7.7  psec  shows  a  dip  the  line  center 
and  a  side  peak  separation  of  about  3A.  Comparison  with  self-absorbed 
profiles  again  indicates  a  40  to  50  percent  error  in  halfwidth.  Since 
radiation  source  contributions  from  the  nearby  exterior  plasma  layers 
have  not  Been  Included  in  these  absorption  profile  calculations  these 
errors  are  at  best  crude  and  reflect  extreme  upper  limits  on  the  measured 
halfwidt  uncertainty,  i ,e.  line  halfwidth  errors  due  to  selective  absorp¬ 
tion  of  an  observed  spectral  line  are  significantly  less  than  50  percent. 

(g)  Wavelength  Calibration  Effect  on  Line  Shift 

There  are  several  systematic  uncertainties  that  can  influence  the 
determination  of  the  Ion  density  In  the  Jet  collision  plasma.  In  order 
for  the  observed  line  shifts  of  the  Cul  4531  transition  to  b’e  meaningful, 
the  magnitude  of  the  shift  must  be  significantly  greater  than  the  error 
In  a  wavelength  position  measurement.  It  has  been  previously  Indicated 
that  uncertainties  in  the  wavelength  differences  from  a  center  wavelength 
value  were  small  for  the  data  points  describing  an  observed  line  profile. 
In  the  vicinity  of  the  Cul  4531  line,  wavelength  value  errors  are  of  the 
order  of  0.05A.  Such  an  error  is  small  compared  to  the  observed  line 
shifts  which  are  about  1A  (see  Table  X) . 


U4.  jump? 


(h)  of  Fitting  Sensitivity  on  Line  Shift 

Furthermore,  It  was  found  that  uncertainties  in  the  measurement  of 
•  line  shift  due  to  the  line  profile  fitting  procedure  were  also  small, 
approximately  0.1A.  That  Is,  variations  of  the  shift  of  the  observed  Cut 
4531  line  of  this  magnitude  could  be  tolerated  In  the  determination  of  a 
best  match  line  profile.  Thus,  the  combined  error  in  a  line  shift  measure¬ 
ment  due  to  the  mechanics  of  the  data  reduction  and  profile  fitting  are 
about  15  percent.  The  correspond ing  uncertainty  In  an  Ion  density  value 
determined  from  a  measured  shift  is  of  the  order  of  15  to  2$  percent  In 
view  of  the  power  law  relations  between  line  shift  and  Ion  density. 

(1)  Precision  of  Stark  Coefficients 

Another  source  of  error  appears  In  the  values  of  the  Stark  coeffic¬ 
ients  for  the  spectral  line  used  In  the  ion  density  calculations.  In  the 
laboratory  measurements  of  these  coefficients  the  wavelength  splitting 
and  shifts  are  very  small  (c.f.  Table  IX) .  In  addition  the  simultaneous 
electric  field  measurements  introduce  another  element  of  uncertainty. 

Hence  the  Stark  coefficients  can  be  subject  to  relatively  large  errors. 

To  ascertain  the  magnitude  of  the  variations  In  a  measured  Stark  coeffic¬ 
ient  a  comparison  of  the  available  Stark  effect  data  for  several  copper 
lines  obtained  by  different  investigators  was  made.  The  results  presented 
as  linear  Stark  coefficients  are  given  in  Table  XV.  For  each  of  the  tab¬ 
ulated  spectral  lines  there  is  agreement  as  to  the  direction  of  the  line 
shift  with  applied  electric  field  for  both  normal  and  parallel  Stark  com¬ 
ponents.  However,  factor  of  two  variations  in  the  magnitude  of  the  line 
shift  per  unit  field  are  seen  to  exist  in  several  comparisons.  It  Is 
unfortunate  that  only  one  of  the  Investigators  (Takamlne)  has  measured 


th*  Stark  splitting  and  shift  for  the  Cul  4531  transition.  Thus  a  direct 
Indication  of  the  reliability  of  the  Stark  coefficients  for  this  spectral 
line  Is  lacking.  If  one  assumes  that  the  variations  In  the  Cul  4531  Stark 
data  used  In  the  Ion  density  determinat ions  are  at  worse  in  error  by  a 
factor  of  two,  then  the  resultant  Ion  density  values  will  be  in  error  by 
approximately  a  factor  of  three, 

(j)  Validity  of  Stork  Coefficient  Measurements 

There  is  another  systematic  uncertainty  in  the  ion  density  calcula¬ 
tions  that  deserve  mention  but  for  lack  of  sufficient  data  cannot  be 
quantitatively  evaluated.  At  such  high  ion  densities  the  mean  intermoie- 
cular  electric  field,  i.e.  Fo  in  eq  (2-23),  is  very  large.  Typically 
for  an  ion  density  of  10  cm  J  this  field  is  about  375  kv/cm.  This  is 
much  greater  than  the  electric  fields  of  the  order  of  50  kv/cm  used  to 
experimentally  study  the  Stark  effect  for  copper  transitions.  Thus  there 
is  a  possibility  that  further  line  splitting  and  shifts  may  result  with 
the  extremely  high  fields  that  exist  in  the  jet  collision  plasma,  Takamlne 
has  observed  the  appearance  of  such  additional  Stark  components  at  high 
field  strengths  for  several  other  copper  spectral  lines.  The  net  effect 
of  these  added  Stark  components  on  the  mean  shift  of  the  observed  Cul  4531 
line  will  depend  in  detail  upon  the  magnitude  and  wavelength  direction  in 
which  the  presumed  line  splitting  progresses. 

(k)  Validity  of  Assumption  of  Negligible  Shift  of  Central  Component 

In  an  effort  to  avoid  large  errors  in  the  determination  of  the  Cul 
4531  line  shift  due  to  the  limited  wavelength  resolution  of  the  time- 
resolved  spectrum,  shifts  relative  to  the  nearby  Cut  4540  transition  were 
measured.  This  implies  that  the  Cul  4540  line  exhibits  a  negligible  shift 


under  the  physical  conditions  existing  in  the  Jet  collision  plasma,  Con¬ 
sequently  the  adequacy  of  a  Cul  4531  line  shift  measurement  as  a  good 
indicator  of  the  copper  Ion  density  in  this  experiment  depends  upon  the 
degree  of  validity  of  this  assumption.  There  are  two  factors  which  indi¬ 
cate  that  Cui  4540  line  shifts  are  indeed  negligible  in  this  experiment. 

The  nearest  perturbing  states  of  the  e  D  upper  level  for  the  Cul  4540 
transition  are  the  x^P°,  D°,  F°  levels  in  the  anomalous  term  configura¬ 
tion  for  copper.  Typical  energy  separations  are  about  7000  un"'  between 
these  perturbing  states  and  the  e  D  level  of  concern.  Since  the  Stark 
shift  Is  inversely  proportional  to  the  separation  of  Interacting  levels 
one  can  therefore  expect  any  shift  in  the  observed  Cul  4$40  line  to  be 
smaller  than  the  corresponding  Cul  4531  shift  for  which  the  energy  sepa¬ 
ration  of  the  nearest  perturbing  level  is,  as  pointed  out  previously, 

2000  cm"1. 

Furthermore  there  is  evidence  for  a  negligible  line  shift  contained 
In  the  previous  analysis  of  the  t ime-resolved  spectral  data.  The  Cul  4651 
transition,  studied  In  conjunction  with  the  Cut  4643  line,  is  a  component 
of  the  same  multi  pi et  as  the  Cui  4540  line.  Hence  the  magnitude  and  direc¬ 
tion  of  Its  wavelength  displacement  in  an  electric  field  should  be  compar¬ 
able  to  that  for  the  Cul  4540  line  (lochte-Hol tgreven,  1938),  Multicom¬ 
ponent  lorentzian  profile  fits  to  the  Cui  4651-4643  pair  were  shown  to 

match  the  observed  line  profiles.  No  line  shifts  of  either  component  were 

p 

Introduced  Into  the  profile  match.  Since  the  e  D  upper  level  of  the  Cul 
4643  transition  is  a  component  of  a  different  multiplet  complex  the  wave¬ 
length  direction  of  its  Stark  splitting  and  shift  pattern  could  be  differ¬ 
ent  from  that  of  either  the  Cul  4651  or  4540  line.  Also  this  level  is 
about  7000  cm  *  removed  from  potential  Interacting  states  and  hence  the 
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magnitude  of  a  Stark  shift  for  this  line  wouid  be  comparable  to  that  of 
the  Cut.  4$40  shift.  Since  no  relative  line  shift  of  the  Cul  4643  line 
was  necessary  to  achieve  a  good  match  to  the  observed  Cut  46$1  profile, 
the  assumption  of  a  negligible  Cui  4540  line  shift  appears  justifiable. 


2.6  Plasma  Temperature  Estimates 

Estimates  of  the  temperatures  that  exist  in  the  collision  plasma 
can  be  obtained  from  previous  results  in  several  Independent  ways.  If 
the  neutral  atoms  are  in  thermal  equilibrium  then  the  relative  popula- 
tions  of  excited  states  are  governed  by  Boltzmann  factors.  The  ratio  of 
the  total  light  intensity  of  two  spectral  lines  of  the  neutral  atom  orig¬ 
inating  from  two  upper  energy  levels  in  thermal  equilibrium  with  energies 
E.,Ej  can  be  written  as 

1  j  /  *  I  =  (vjgjfj)/(v]gjf .)  exp(E  j-Ej/kT)  (2-40) 

where  I  represents  the  total  line  intensity,  v  the  center  frequency, 
f  the  oscillator  strength  for  the  transition,  g  the  statistical  weight 
of  the  upper  energy  level,  k  the  Boltzmann  constant,  and  T  the  abso¬ 
lute  temperature,  (tochte-Ho I tgreven,  1958;  Donohue  and  Majkowski,  1962; 
Griem,  1964).  The  total  light  intensity  is  given  by  the  integral  of  the 
spectral  profile  over  the  full  extent  of  the  line,  i.e. 

1  "  J  dV  (2-41) 

For  a  resonance -broadened  fine  whose  profile  is  described  by  eq  (2-7) 

the  total  light  intensity  is  determined  by  the  product  A  J(v  )  where 

o 

A  is  the  resonance  haifwidth  and  j(vQ)  the  peak  intensity  at  the  center 
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frequency  •  The  ratio  of  total  light  Intensities  for  two  resonance- 
broadened  lines  is  then  given  by 


y*i  ■  »jj,<V/Vi<vi) 

where  vj»vj  are  cen{:er  frequencies  of  the  two  I  lines* 
oscillator  strengths  for  the  spectral  transitions  are  known 
ture  T  can  be  calculated  from  the  intensity  ratios  of  eqs 


(2-42) 

Thus  »  I f  the 
then  a  tempera- 
(2-40)  and 


(2-42) * 


Values  of  the  temperature  at  5 - 7 ►  6*7.  and  7-7  psec  after  initial 
jet  impact  have  been  obtained  from  the  measured  resonance  halfwidths  and 
relative  peak  intensities  of  the  observed  Cuf  4651-4643  and  Cui  453l“4525 
pairs-  The  resultant  estimates  given  in  Table  XVI  and  Fig*  28  indicate 
that  the  temperature  of  the  region  from  which  the  radiation  originates 

If  Q 

is  of  the  order  of  10  K  and  decreases  with  time- 

The  degree  of  ionization  in  a  plasma  in  thermal  equilibrium  is 

governed  by  the  Saha  equation*  Estimates  of  the  temperature  can  also  be 

obtained  from  the  Saha  equation  relating  the  equilibrium  distributions 

★ 

of  ions,  electrons,  and  neutral  atoms  as 


N? 

> 

N 


2  U,(T) 
U  (T) 


h 


(2-43) 


where  Nj  is  the  number  density  of  ions,  N  the  number  density  of 
neutral  atoms,  Uj(T)  and  U(T)  the  partition  function  for  the  ion  and 
neutral  atom,  E  the  ionization  potential,  T  the  absolute  temperature, 
and  the  constants  m,k,  and  h  have  their  usual  meanings*  The  ionization 


Equai  ion  and  electron  densities  are  assumed* 
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potentla,  Is  given  by  the  first  Ionization  potential  of  the  free  atom 
reduced  by  the  Coulomb  Interaction  energy  of  the  plasma  AEw  where 

,2  f4n  Tl/3  -7  1/3 

ttm  **  (y  N|  }  «  7-0  x  10  7{N})  (2-44) 

o 

for  N,  In  an  3  and  AE  In  ev,  (Allen,  1963)* 

Temperature  estimates  have  been  obtained  from  the  Saha  equation 
using  the  measured  jet  collision  plasma  ion  and  neutral  atom  densities, 
and  are  also  given  In  Table  XVI  and  Fig.  28  at  several  times  after 
Initial  impact.  Temperatures  are  seen  to  be  of  the  order  of  6  x  103  °K 
and  show  a  moderate  decrease  with  time  after  impact- 

The  shape  of  the  background  continuum  in  the  time- resolved  spectrum 
can  also  yield  a  temperature  estimate-  For  temperatures  between  5000  and 
7000°K  a  black  body  spectrum  has  a  peak  Intensity  that  falls  within  the 
4400-5300A  range  of  the  time-resolved  spectrum.  In  particular  the  wave¬ 
length  at  peak  intensity  Xffl  ,  is  given  by  Wien's  law 

X  T  -  0.2896  (2-45) 

m 

for  X  in  cm  and  T  in  °K  (Joos,  1951)-  Ascribing  a  black  body  shape 
m 

to  the  observed  continuum  in  the  t ime-reso I ved  spectrum  the  peak  intensity 

occurred  at  a  wavelength  in  the  vicinity  of  5000A  at  a  time  of  6.7  jisec 

3  o 

after  impact.  The  as^ciated  black  body  temperature  is  about  5-8  x  )0J  K 
and  Is  In  close  agreement  with  previous  temperature  estimates  at  this  time- 
Jet  collision  pias~ia  temperatures  have  also  been  obtained  from  the 
fluid-dynamic  collision  model-  The  variation  with  time  of  the  tempera¬ 
ture  at  the  center  of  impact  and  at  the  spatial  region  in  the  expanding 
sheet  from  which  the  observed  radiation  originates  is  illustrated  In 


Fig,  28-  Both  curves  show  the  expected  decrease  In  temperature  with  In¬ 
creasing  time,  but  the  calculated  temperature  values  are  about  a  factor 
of  seven  lower  In  magnitude  that  the  temperature  estimates  previously 
described.  As  will  be  shown  this  discrepar  /  can  be  expected  since  the 
fluid-dynamic  calculations  infer  heavy  particle  temperatures,  whereas 
the  estimates  derived  from  spectral  observations  are  electron  temperatures • 

The  temperature  of  a  gas  determines  the  distribution  of  Internal 
energy  for  all  the  different  constituents  present  in  the  gas  In  accord 
with  a  Boltzmann  distribution.  The  internal  energy  Is  composed  of  the 
kinetic  energy  of  the  free  electrons,  Ions,  and  neutral  atoms,  and  of  the 
excitation  energy  of  the  ions  and  neutral  atoms*  The  partition  of  the 
internal  energy  is  determined  by  collisional  and  radiative  processes  in 
the  plasma.  If  collisions  are  the  predomlnent  mechanism  in  determining 
the  state  of  the  gas  then  thermal  equilibrium  can  be  expected,  (Dewan, 
1961). 

Ambartsumyan  (1958)  has  calculated  the  ratio  of  the  rate  of  Ioni¬ 
zation  by  electron  Impact  to  the  rate  of  ionization  due  to  photo-absorp¬ 
tion*  For  the  electron  densities  and  temperatures  found  to  exist  in  a 
jet  collision  plasma,  the  electron- induced  ionization  rate  exceeded  that 
due  to  radiative  absorption  by  about  a  factor  of  ten*  Due  to  their 
relatively  low  velocities  and  hence  collision  frequencies,  heavy  particle 
col  1  is  Ion- ion i zat ion  rates  are  also  negligible*  Furthermore  Griem  (I96h) 
has  indicated  that  if  the  relative  populations  of  the  neutral  atom  ground 
state  and  first  excited  level  are  governed  by  a  Boltzmann  factor  contain¬ 
ing  the  electron  temperature,  then  the  relative  pupulations  of  all  higher 
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levels  and  the  degree  of  ionization  are  also  determined  from  appropriate 
thermal  eQuinbriim  relations  in  terms  of  the  electron  temperature- 
The  condition  that  the  col  1 1 s lonal -exci tat  ion  rate  of  the  first  excited 
level  exceed  the  ratiative-population  rate  leads  to  a  minimum  value  of 
the  electron  density  at  a  given  temperature  that  is  required  to  assure 
thermal  equilibrium  of  aii  the  excited  states  of  the  atom-  The  electron 
densities  and  temperatures  in  a  jet  collision  plasma  easily  fulfill  this 
minimum  requirement-  Thus  the  electrons.  Ions  and  atoms  of  the  collision 
plasma  are  in  thermal  equilibrium  and  the  above  temperature  estimates 
derived  from  relative  line  intensity  measurements  and  neutral  atom-ion 
equilibrium  relations  are  to  be  interpreted  as  electron  temperatures- 
Electron-electron  collisions  will  establish  a  Maxwellian  velocity 
distribution  for  the  electrons  while  the  heavier  particle  collisions 
will  also  yield  a  velocity  distribution  for  the  ions  and  atoms  but  at 
a  temperature  that  may  differ  from  the  electron  temperature-  The 
electrons  will  reach  an  equilibrium  distribution  more  rapidly  than  the 
heavier  atoms  due  to  their  greater  velocity  and  collision  frequency. 
Finally,  energy  equlpart » t ion  between  all  the  particles  comes  about  by 
electron-ion  and  electron-atom  collisions-  The  time  scale  in  which 
energy  equipart i t ion  corns*  about  is  much  greater  than  the  time  scale 
in  which  electrons  and  heavy  particles  establish  their  Maxwellian  dis¬ 
tribution  due  to  the  small  energy  transfer  in  an  electron-heavy  particle 
encounter,  a  consequence  of  the  large  particle  to  electron  mass  ratio¬ 
using  Griem’s  estimate  of  the  equilibration  time  for  the  electrons  to 
attain  their  distribution  via  electron-electron  collisions  and  the 
electron  densities  and  temperatures  found  in  this  work,  the  estimated 


T 


sec.  To  permit 


electron  equilibrium  time  Is  of  the  order  of  10 
colllslonal  excitation  of  the  atoms  and  ionization  to  become  established 
one  should  allow  a  somewhat  longer  time,  about  10  ^  sec-  Collisions 
of  the  Ions  and  neutral  atoms  with  electrons  results  In  an  energy  equi- 

partstion  but  at  much  later  times-  The  equilibrium  time  for  electrons 

-8 

and  Ions  Is  found  from  Gr Jem's  estimates  to  be  about  2  x  10  sec- 
whereas  for  atoms  It  Is  about  2  x  10  ^  sec-  Since  these  time  estimates 
are  in  the  nature  of  a  time  constant  one  would  expect  electron- ion 

-7 

equilibrium  over  a  time  span  of  10  sec,  and  electron-atom  equilibrium 

-5 

in  10  sec- 

Thus,  within  the  time-resolution  of  these  experiments  and  for  the 
time  period  after  jet  impact  studied,  the  electrons  and  Ions  have  the 
same  kinetic  temperature  indicated  In  the  above  measurements.  The 
neutral  atoms  on  the  other  hand  have  not  fully  achieved  the  measured 
electron  temperatures  In  this  time  period-  Collision  temperatures  cal¬ 
culated  from  the  fluid-dynamic  collision  model  pertain  to  the  neutral 
atoms  In  th~  plasma  and  were  previously  seen  to  be  lower  than  the  cor¬ 
responding  electron  temperatures  and  decreasing  with  time*  Both  neutral 
atom  and  electron  temperatures  however  appear  to  converge  at  somewhat 
later  times  in  the  experiment  as  indicated  in  Fig-  28. 


Part  Mi  Discussion  of  Results 


An  experiment*!  study  of  the  collision  of  two  Impinging  shaped- 
charge  jets  has  led  to  some  understanding  of  a  high  velocity  impact.  The 
gross  mechanical  features  of  such  a  collision  have  been  clearly  Illus¬ 
trated  and  estimates  of  the  neutral  atom  and  Ion  densities  In  the  colli¬ 
sion-generated  plasma  we$e  derived.  Relative  oscillator  strengths  for 
two  neutral  copper  lines  have  been  determined  from  line  halfwidths  in  the 
time-resolved  spectrum.  The  results  of  these  experiments  are  to  be  com¬ 
pared  with  conclusions  of  other  investigators. 


3.1  Pertinent  Hyperveioc ity  Impact  Theory 

For  the  jet  velocities  in  these  experiments  the  jet  material  will 
act  dynamically  as  a  fluid  upon  impact.  Goldsmith  (I960)  has  briefly 
considered  the  analogous  behavior  of  a  cylinder  striking  a  plane  rigid 
target  at  lower  speeds  than  those  involved  in  the  work  presented  here. 

As  a  result  of  this  collision  two  regions  are  formed  in  the  cylinder. 
Adjacent  to  the  rigid  wall  the-materlal  is  fluid-like  and  flows  outward 
along  the  target  Interface,  thereby  simulating  the  formation  of  a  radially 
expanding  sheet.  Upstream  this  fluid  region  is  separated  by  a  shock  wave 
from  the  incoming  undisturbed  material.  Goldsmith  has  indicated  that  the 
entire  body  can  be  treated  as  a  fluid  and  the  usual  hydrodynamic  shock 
relations  (Courant  and  Friedrichs,  19^8)  used  to  describe  material  condi¬ 
tions  In  the  vicinity  of  the  shock  wave  reflected  from  the  rigid  wall  Into 
the  advancing  material. 

GUvarry  and  Hill  (1956)  have  applied  the  hydrodynamic  shock 


equations  to  their  mode!  of  a  high  speed  projectile  collision.  They  have 
calculated  pressure,  density,  and  temperature  behind  the  reflected  shocks 
created  by  the  collision  of  two  plane-para  1  lei  solid  regions  using  an 
equation  of  state  obtained  from  the  statistical  Thomas-Fermi  model  of  the 
atom.  Iron  and  silicon  were  chosen  for  these  investigations.  This  one- 
dlmensionai  shock  model  is  applicable  to  a  jet  collision  for  a  suffici¬ 
ently  short  time  after  impact  before  rarefaction  waves  from  the  cylind¬ 
rical  jet  boundaries  are  experienced  at  the  original  impact  point.  Sue 
to  these  rarefaction  waves  the  shock  wave  progressing  up  the  jet  stem  is 
no  longer  planar  and  a  one-d imens iona 1  collision  model  becomes  unrealis¬ 
tic,  At  impact  velocities  in  the  range  of  10  mm/psec  Gilvarry  and  Hill 
estimate  the  pressure  behind  the  reflected  shock  wave  to  be  of  the  order 
of  10  megabar,  temperature  somewhat  greater  than  10  K,  and  the  density 
compression  ratio  about  1,5.  Furthermore  the  resultant  shock  tempera¬ 
tures  were  shown  to  be  in  excess  of  the  estimated  melting  temperatures  at 
such  pressures  and  densities  so  that  at  least  melting  of  the  impacting 
material  is  to  be  expected.  These  extreme  conditions  representative  of 
the  earlier  moments  of  a  jet  collision  are  soon  reduced  upon  arrival  of 
rarefaction  waves  and  radial  expansion  of  the  impact  material. 

Lukas ik,  et  al  (196k)  have  created  a  fluid-dynamic  mode!  descrip¬ 
tive  of  a  head-on  jet  collision.  The  early  stage  of  the  collision  was 
represented  by  the  impact  of  a  plane-parallel  solid  region  and  a  station¬ 
ary  rigid  wall  In  a  similar  fashion  to  the  work  of  Gilvarry  and  Hill. 

Upon  impact  a  shock  wave  progresses  into  the  solid  region  and  the  material 
between  this  shock  front  and  t^e  rigid  wall  comes  to  rest.  For  compara¬ 
tive  purposes  several  equation  of  state  representations  were  used  in  their 


calculations  of  material  properties  In  the  shocked  region.  Tfllotson's 
equations  of  state  were  considered  to  be  most  applicable  for  this  work 
since  they  are  based  on  both  shock  wave  data  up  to  about  one  megabar  and 
theoretical  results  of  Thomas-Ferml -Dirac  theory  at  pressures  In  the 
vicinity  of  100  megabars;  thus  the  material  properties  at  the  10  megabar 
pressures  encountered  here  were  probably  correctly  described.  For  a 
copper  jet  with  a  velocity  of  8  mm/psec  the  impacted  material  was  found 
to  hove  a  Hugoniot  pressure  of  ubout  12  megabar,  temperature  of 
1,3  x  10**  °K,  end  a  density  compression  ratio  of  1,9,  Since  the  shaped 
charge  jet  tip  has  a  small  diameter  (of  the  order  of,  2  mm)  these  condi¬ 
tions  exist  only  for  a  short  period  of  time,  typically  less  than  0.5  psec 
after  which  the  ensuing  expansion  of  the  collision  products  will  termi¬ 
nate  this  early  phase  of  the  collision  process. 

The  second  phase  of  the  jet  impact  model  ts  concerned  with  the 
expansion  characteristics  of  the  collision,  where  results  of  the  one¬ 
dimensional  calculations  were  employed  as  Initial  conditions.  The  expan¬ 
sion  of  a  one-dimensional  slab,  of  a  cylinder,  and  a  sphere  of  radius 
comparable  to  the  initial  jet  tip  radius  were  studied.  Three  important 
results  were  derived  from  these  considerations;  namely,  the  colliding  jet 
material  will  be  vaporized  upon  expansion;  the  velocity  of  the  edge  or 
front  of  the  expanding  region  Is  the  same  for  either  geometry,  and  is 
about  17  mm/psec.  The  density,  pressure,  and  temperature  of  the  collision 
products  is  extremely  high  in  the  vicinity  of  the  impact  point;  for  the 
expanding  material  near  the  front  they  fall  off  sharply  as  (I  -  r/vt)n  , 
where  r  is  the  radial  distance  from  the  collision  center,  t  is  the 
time  after  impact,  and  v  the  velocity  of  the  front.  The  exponent  n 


has  the  values  k,  6  and  2  for  the  material  density,  pressure,  and  temper- 

'k 

ature  profiles  respectively. 


3.2  Discussion  of  Experimental  Results 

These  tatter  theoretical  conclusions  have  been  verified  fey  the 
shaped  charge  Jet  collision  experiments.  The  rich  copper  line  spectral 
content  of  both  the  t ime- Integrated  and  t Sme-resoi ved  records  show  that 
copper  vapor  Is  present  in  the  collision  products.  The  calculated  time 
for  the  expanding  products  to  reach  the  vacuum  chamber  wall  was  found  to 
be  in  excellent  agreement  with  the  time  at  which  the  neutral  atom  dens  ity 
Is  a  minimum.  A  subsequent  increase  In  the  neutral  atom  density  is 
attributed  to  the  pile-up  of  material  at  the  window  accompanying  the 
creation  of  an  Inward  moving  compressive  shock  created  by  the  wall  colli¬ 
sion.  Furthermore,  the  extremely  high  pressure,  density,  and  temperature 
conditions  that  are  predicted  to  exist  are  made  evident  by  the  intense 
continuum  radiation  that  appears  on  the  early  portion  of  the  time-resolved 
spectral  record,  the  appearance  of  an  extremely  bright  regipn  about  the 
point  of  collision  on  the  framing  camera  record,  and  the  presence  of 
Highly  excited  neutral  and  ionized  copper  atoms  with  up  to  2$  ev  excita¬ 
tion. 

The  emission  lines  of  several  spectral  transitions  in  the  time- 
resolved  spectrum  were  shown  to  possess  instrumentation  and  resonance- 
broadened  contours.  Integrated  Lorentzian  line  profiles  which  combine 

*This  density  profile  was  the  one  used  In  the  previously  discussed  spec¬ 
tral  line  absorption  studies. 


the  effect*  of  both  broadening  mechanism  were  seen  to  match  the  observed 
line  shapes  In  the  time  Interval  studied.  Values  of  the  spectrograph  slit 
width  determined  at  several  wavelengths  were  found  to  be  in  good  agreement 
with  calculated  values  obtained  from  an  optical  ray  tracing  through  the 
direct  view  prism  streak  camera  combination.  Estimates  of  the  neutral 
atom  density  for  several  microseconds  after  initial  jet  impact  were 
derived  from  resonance-ha l fwidth  measurements  of  three  observed  spectral 
lines  with  known  oscillator  strengths.  The  computed  densities  were  found 
to  be  in  close  agreement,  of  the  order  of  4  x  10  "  per  cm  over  the  time 
span  considered.  Since  these  spectral  transitions  arise  from  different 
excited  levels  and  configurations  of  the  neutral  copper  atom  and  have 
different  oscillator  strength  values,  this  agreement  indicates  that  the 
assumptions  and  methods  used  in  the  neutrai  atom  density  calculations  are 
valid. 

Lukas ik,  et  a  I  have  further  attempted  a  numerical  solution  of  the 
radiative  transport  equation  for  the  jet  collision  products  using  results 
of  their  hydrodynamic  calculations  and  theoretical  optical  opacity  esti¬ 
mates.  They  have  shown  that  the  observed  spectral  radiation  has  its 
origin  in  a  narrow  region  near  the  front  of  the  expanding  collision 
products  and  have  found  the  particle  density  in  this  region  to  be  approx¬ 
imately  two  to  three  times  larger  than  the  atom  densities  derived  from  the 
spectral  line  halfwidth  measurements.  In  view  of  the  uncertainty  in  the 
opacity  estimates  the  use  of  a  one-dimensional  mode)  neglecting  material 
expansion  normal  to  the  outward  moving  sheet,  and  the  fact  that  the  spec¬ 
tral  data  from  which  atom  densities  are  obtained  are  averaged  over  time 
and  spaca  due  to  the  finite  extent  of  the  spectrograph  entrance  slit, 


the  agreement  In  density  estimates  Is  considered  to  be  reasonably  good. 

The  magnitude  of  the  copper  Ion  density  in  the  jet  collision  plasma 
has  been  determined  from  the  measured  wavelength  shift  of  a  neutral  copper 
line  for  which  some  Stark  effect  data  is  available.  Since  the  nature  of 
the  variation  of  Stark  splitting  and  shifts  with  electric  field  strength 
is  not  well  established,  both  linear  and  quadratic  Stark  effect  character¬ 
istics  were  considered  and  their  influence  upon  the  calculated  ion  densi¬ 
ties  has  been  discussed.  In  addition,  a  comparison  of  the  results  of 
Holtsmark's  theory  for  ionic  perturbers  and  numerical  evaluations  of  the 
folded  distribution  to  include  the  additional  effect  of  electrons  was 
made,  assuming  a  linear  Stark  effect  for  the  affected  spectral  line.  Use 
of  Holtsmark's  results  was  seen  in  this  case  to  yield  a  slightly  higher 
ion  density  estimate.  Evidence  exists  however  which  indicates  that  a 
quadratic  electric  field  dependence  is  more  realistic.  With  this  assump¬ 
tion  ion  densities  were  calculated  using  Holtsmark's  theory.  Copper  ion 
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densities  of  the  order  of  4  x  10  cm  were  found  to  exist  In  the  impact 
plasma  indicating  that  the  degree  of  ionization  was  about  one  percent 
over  the  time  span  under  consideration. 

Oscillator  strength  values  for  the  Cut  4525  and  4643  spectra!  trans¬ 
itions  of  neutral  copper  have  been  determined  from  comparative  line  half¬ 
width  measurements  using  as  a  basis  the  Cul  4651  transition  whose  f-number 
has  been  previously  established.  These  results  are  to  be  compared  with 
the  only  other  known  f-number  measurements  for  these  spectral  lines 
obtained  by  Allen  and  Asaad  (1957).  Table  XV l|  contains  the  oscillator 
strength  values  for  the  Cul  4525  and  4643  transitions  found  by  Allen  and 
Asaad  and  from  this  work.  Comparative  f-values  for  Cul  4651  from  Corliss 
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and  Bczman  (1962),  and  from  Alien  and  Asaad  are  also  given.  Corliss 
(1962)  has  indicated  the  existence  of  a  large  systematic  discrepancy  In 
Allen  ar d  Asaad's  results  due  to  their  incorrect  use  of  a  constant 
normalization  factor  to  establish  absolute  oscillator  strengths  from 
relative  measurements.  Fu  ;hermore,  after  correction  for  this  system¬ 
atic  difference  a  residual  discrepancy  between  Allen  and  Asaad's  results 

and  Corliss  and  Bosnian1  s  measurements  of  about  a  factor  of  three  still 

■  sv 

remained.  In  view  of  these  uncertainties  it  is  difficult  to  assess  the 
agreement  between  comparative  oscillator  strength  determinations  given 
in  Table  XV 1 1  .  At  best  one  can  say  that  the  results  are  not  grossly 
out  of  line. 


3.3  Further  Extensions  and  Applications 

There  are  a  number  of  areas  of  physical  interest  to  which  the 
results  of  shaped  charge  jet  collision  experiments  can  be  of  value. 

A  high  velocity  jet  impact  represents  a  new  technique  for  the  production 
of  a  plasma  with  a  high  particle  density  and  moderate  degree  of  ioniza¬ 
tion.  The  possibility  of  generating  thermor uc lear  power  from  impacting 
materials  as  suggested  by  Harrison  ( I 963 )  is  of  great  interest  and  can  be 
viewed  as  a  goal  of  such  efforts.  The  feasibility  of  such  experiments 
must,  however,  await  the  development  of  still  higher  speed  projectiles. 


Alien  and  Asaad  further  estimate  their  f-numbers  to  be  accurate  to  within 
about  a  factor  of  two  due  to  their  experimental  uncertainties. 


Some  technological  progress  in  this  direction  has  been  made  hy  Walsh, 
et  al  0953)>  and  Jameson  (1963)  In  their  design  of  explosive  accelera¬ 
tors  and  by  Walsh,  et  al  and  Lukaslk,  et  al  In  the  choice  of  projectile 
mater ia is . 

A  most  obvious  contribution  to  the  field  of  atomic  spectroscopy  and 
allied  areas  in  astrophysics  is  in  the  determination  of  oscillator 
strengths  for  spectral  lines  from  excited  energy  levels  of  metallic  atoms. 
At  present  oscillator  strength  values  are  determined  either  from  electric 
arc  experiments  (Locbte-Hol tgreven,  1958)  or  by  numerical  computations  of 
questionable  ref lability  (Allen,  1963).  Spectroscopic  bbservat Ions  of  a 
hypervelocity  projectile  impact  can  offer  both  an  independent  experimental 
verification  of  existing  results  and  also  a  continuing  program  of  new 
f-number  determinations  for  selected  materials.  With  the  advent  of  higher 
speed  shaped  charge  devices  and  the  anticipated  increase  in  the  collision 
plasma  ion  density,  measurements  of  ionic  spectral  line  oscillator 
strengths  by  the  methods  contained  in  this  work  appear  possible.  This  is 
of  particular  importance  since  oscillator  strength  measurements  for  ionic 
spectral  transitions  are  difficult  to  achieve  with  present  day  techniques. 

Shaped  charge  jet  collision  experiments  may  also  serve  as  a  means 
of  obtaining  high  pressure  and  density  equation  of  state  measurements  for 
materials  that  can  be  made  to  form  a  jet  under  explosive  loading.  The 
extremely  high  pressures  and  densities  experienced  at  the  moment  of  jet 
impact  are  much  greater  than  that  obtained  in  present  day  experiments 
specifically  oriented  towards  equation  of  state  determinations,  (Walsh 
and  Christian,  1955).  With  the  advent  of  higher  velocity  projectiles 
further  extensions  in  the  range  of  these  quantities  will  be  attainable. 
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Unfortunately  direct  pressure  and  density  measurements  are  not  possible 
and  hence  conditions  In  the  vicinity  of  the  impact  center  can  only  be 
Inferred  from  a  detailed  analysis  of  spectroscopic  observations  and  appli¬ 
cation  of  a  f luld-dynamic  mode?  of  a  jet  collision.  Nevertheless,  one  can 
proceed  In  a  manner  similar  to  that  used  In  this  work  and  correlate  exper¬ 
imentally  derived  particle  densities  with  analogous  results  from  the  model 
calculations.  These  calculations  wouid  be  performed  parametr Icaily  for 
the  various  equations  of  state  under  study  in  order  to  determine  the  best 
representation.  An  alternate  procedure  would  be  to  determine  experi¬ 
mentally  the  time  after  initial  impact  when  the  material  front  from  the 
jet  collision  reaches  the  vacuum  chamber  wail.  The  constant  velocity  of 
the  front  derived  from  the  wall  collision  time  can  be  related  to  Initial 
conditions  that  exist  during  the  early  stages  of  the  collision  by  the 
equations  of  motion  descriptive  of  the  expansion  phase.  These  initial 
conditions  in  turn  are  determined  by  the  usua 1  hydrodynamic  shock  rela¬ 
tions  and  the  material  equation  of  state. 

There  is  one  area  of  current  interest  for  which  the  experimental 
and  analytical  procedures  employed  In  this  work  can  be  directly  applied. 
Extensive  research  efforts  in  exploding  wire  phenomena  have  led  to  a 
basic  understanding  of  the  important  processes  attendant  upon  a  wire 
explosion.  A  large  amount  of  theoretical  development  and  experimental 
results  pertaining  to  the  electrical  characteristics  and  dynamic  proper¬ 
ties  of  the  wire  are  available  (Chace  and  Moore,  1959,  1962).  However 
scant  attention  has  been  given  to  the  information  that  can  be  obtained 
from  spectroscopic  observations.  T ime- integrated  spectral  investigations 
and  to  a  limited  extent  time-resolved  spectral  studies  have  been  made. 
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but  the  analysts  of  these  records  has  been  restricted  to  the  Identified* 
tlcn  of  spectral  lines  characteristic  of  the  wire  material  (Anderson, 

1920;  Smith,  1925;  Relthel,  et  al,  1959)*  Exploding  wire  spectra  are 
rich  In  detail  and  the  results  derivad  from  an  analysis  of  the  spectral 
data  along  the  I  Ines  of  approach  used  In  this  work  would  be  an  Important 
contribution  to  the  knowledge  end  understanding  of  the  phenomenon. 

Present  day  plasma  physics  applications  and  astrophysics!  require* 
ments  have  placed  an  obvious  emphasis  on  obtaining  Stark  effect  measure* 
ments  for  certain  elements  such  as  hydrogen,  helium,  argon,  etc.  Determ¬ 
inations  of  the  ion  density  In  a  jet  collision  plasma,  however,  have  made 
use  of  experimentally-determined  Stark  splitting  and  shifts  of  a  selected 
spectral  line  of  copper.  The  lack  of  available  Stark  effect  data  for 
such  uncommon  plasma  const ItuantS  Is  to  be  noted.  It  Is  hoped  that  this 
and  future  work  In  hyperveloc I t y  impact  experiments  will  provide  an  Incen¬ 
tive  to  obtain  Stark  effect  measurements  for  spectral  lines  of  these 
unusual  plasma  components. 
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TABU  1 


Explosive  Shaped  Charge  Specifications 


Designation 


RC-338,  E.l.  DuPont  de  Nemours 


Material 


Explosive  RDX 
Liner  Copper 


Dimens ion 
Explosive 
Liner 


Cylindrical,  2  in.  dia,  3-3/8  In.  length 

0.063  in.  wall  thickness,  45°  total 
included  cone  angle 


Hess 


Explosive  approximately  200  grams 
Liner  approximately  44  grams 


TABLE  i! 


Summary  of  Jet  Velocity  Measurements 


Firing 

Average  Velocity 

no. 

(mm/nsec) 

56 

7.96 

57 

7.98 

59 

7.86 

99 

7-68 

Mean  7-9  -  0. I 


ram/ nsec 


TABLE  lit 


Absorption  Lines  in  the  T I me- Integrated  Spectrum 
of  a  Jet  Collision  in  Air 


Observed 

Wave  length 
(A) 

As  $ i gned 

Copper 

Wavelength 

(A) 

Lower  Energy  Level 
of  the  Transit  Ion 
(ev) 

koz 

4023 

3.78 

4067 

4063 

3.82 

4276 

4275 

4.84 

4513 

4507 

5.57 

4509 

5.24 

4580 

4587 

5.10 

4646 

4651 

5.07 

4682 

4675 

5.15 

4697 

5.24 

4703 

4697 

5.24 

4704 

5.10 

TABLE  IV 


Emission  Lines  in  the  T ime-Resolwed  Spectrum 
of  a  Jet  Collision  in  Helium 


Observed 

Cul  Line 

(A) 

Upper  Level 
Des i gnat  ion 

Upper  Energy 
Level 

(cm"*)  (ev) 

4507.35 

4509.37 

4513.19 

li 

e  u  1/2 

64472.3 

7.99 

4525.11 

4530.79 

4539-70 

f2D  1  1/2 

6s 2S  1/2 
ek0  1  1/2 

63584.6 

7.88 

4586.97 

eS>  2  1/2 

62948.3 

7.80 

4642.58 

4651.12 

e20  1  1/2 

J*0  2  1/2 

62403.3 

7,73 

4674.72 

4677.34 

e40  2  1/2 

62948,3 

7.80 

4697 .49 
4701.71 

4704.59 

e40  3  1/2 

62403.3 

7.73 

4767.49 

4776.22 

f 2D  2  1/2 

67142.7 

8.32 

4794.00 

4797.04 

e20  1  1/2 

62260. 1 

8.09 

4842.29 

4866.10 

f20  2  1/2 

67142.7 

8.32 

TABLE  IV  (continued) 


Observed 

Cul  Line 

(A) 

Upper  Level 

Oes i gnat  ion 

Upper  Energy 
Leve  1 

(cm“!)  (ev) 

5016.61 

L 

e  0  1/2 

64472.3 

7.99 

5034.36 

e4D  1  1/2 

63584.6 

7.88 

5076. 17 

e20  2  1/2 

64657.8 

8.01 

5105.54 

5111.91 

5115.49 

4p2P° 1  1/2 

30783.7 

3.82 

5144.12 

5153.24 

5158-36 

4d2D  1  1/2 

49535.2 

6.20 

5212.78 

5218.20 

5220.07 

4d2D  2  1/2 

49942.1 

6.19 

5292.52 

e40  3  1/2 

62403.3 

7.73 
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TABU  V 


Values  of  Emulsion  Calibration  Curve  Slope  Y 


Zn  4680 

Zn  4722 

Hg~Cd  4800 

0.90 

1.45 

1.13 

1.03 

1. 10 

0.84 

1.00 

0.98 

1.28 

0.99 

Zn  4810 

He  5016 

Hg-Cd  5086 

1.42 

1.01 

1.37 

0.99 

0.80 

0.72 

1.11 

1.17 

0.90 

1.04 

Heen  Value 

1.0-  0.1  density 

units 

TABLE  VI 


Values  of  the  Constant  In  Eq. 


Statistical  Theory 

Margansu  and  Watson  ( 1936} 


Interruption  Theory 

Weisskopf  (193-) 

Furssow  and  WJassow  (1936) 
Lindholm  (1933) 


(2-8) 


TT/3 


V3 

Tt/2 


TABLE  VI! 
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TABLE  V'HI 

Relative  Intensities  as  Determined 


from  Multicomponent  li 

ntegrated 

Lorentzlan 

Prof  1 le 

Fits 

Transition 

><ne 

sec) 

4651 

4643 

4540 

4531 

4525 

5.7 

1.00 

0.38 

1.00 

0.44 

0.34 

6.7 

1.00 

0.23 

1.00 

0.56 

0.035 

7.7 

1.00 

0.18 

1.00 

0.39 

0.075 

m*r 


T 


TABLE  IX 


Stark  $pi 

itting  and  Shift 

for  Selected  Neutral 

Copper  Lines* 

Wavelength 

(A) 

p-€omponent 

(A) 

n-Component 

(A) 

Maximum  Field 

(kv/cm) 

453) 

+0.04 

+0.03 

44.0 

5153 

-0.10 

-0.05 

253 

5218 

-0.09 

-0.07 

253 

5220 

__ 

-0.06 

.. 

From  Takamine,  (I919);  p-  and  n-  refer  to  direction*  parallel 
and  norma)  to  the  applied  field. 


TABLE  X 


Cul  4531  Wavelength  Shift  and  Stark  Component 
Separation  as  Determined  from  Integrated 
Lorentzian  Profile  Fits 


Time 

Mean 

Shift 

6 

C 

Component 
Separat ion 
d 

(nsec) 

(A) 

(A) 

5.7 

1.3 

0.36 

6.7 

1.5 

0.42 

7.7 

1.0 

0.29 

TABLE  XI 


Oscillator  Strengths  for  Several  Cul  Lines* 


Wavelength 

(A) 

g  f 

9 

f 

4530.78 

0.20 

2 

0.10 

4651.12 

1.8 

8 

0.23 

5153.24 

1.9 

4 

0.48 

* 

From  Corliss  and  Bozmen  (1962). 
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TABLE  XM 


Neutral  Copper  Atom  Density  In  the  Jet  Collision  Plasma 


Time 

Neutral 

Copper  Density 

(x  1 0  ^  cm  ^  ) 

(|A$ec) 

5153 

4651 

4531 

5-7 

3.6 

52 

90 

6.7 

4.1 

32 

4. 1 

77 

3-2 

2.9 

9-7* 

8.7 

2.8 

2.8 

97 

2.7 

3.2 

10.7 

3.6 

4.3 

H.7 

3.4 

56 

12.7 

4.9 

7.0 

* 

Not  shown  in  Fig.  26. 


TABLE  XII! 


Oscillator  Strengths  for  Cul  4643  and  Cul  4525 


T  rans i t i on 

f-number 

4643 

0.19 

0.22 

0.20 

Weighted  Mean 

0.20  t 

0.01 

4525 

0.31 

0.20 

0.48 

Weighted  Mean 

0.33  - 

0.08 

The  relative  intensities  of  the  observed  spectral  lines 
given  in  Table  VIII  were  used  to  weight  the  individual 


f-nuaibers . 


TA8LE  XIV 


Copper  Ion  Density 

In  the  Jet 

Collision  Plasma 

Ion  Density 

Linear  Stark  Effect 

Quadratic  Stark  Effect 

4 

Hoitsmark 

Folded 

Hoitsmark 

Time 

Theory 

D i str i but  ion 

Theory 

(nsec) 

(cm'3) 

(cm"3) 

(cm'3) 

5.7 

4.9  ID18 

4.2  I0,S 

3.3  1017 

6.7 

6.2 

5.2 

37 

7.7 

33 

2.8 

2.7 

TABU  XV 


TABLE  XVI 


Jet  Collision  Plasma  Temperatures 


Time 

(nsec) 

Relative  Line 
intens i ty 

Cul4651-4643  CuI4531-4525 

Saha 

Equation 

Black 

Body 

5-7 

12,000 

14,000 

6250 

- 

6-7 

6,000 

5.200 

6220 

5800 

77 

4,000 

6,800 

61 30 

- 

Temperatures  In  °K 


TABLE  XV  M 


Comparison  of  Oscillator  Strength  Values 


Cu! 

This 

A! len  and 

Corliss  . 

Transition 

Work 

Assad 

Bozman 

4651 

1.7 

0.23 

4643 

C.  20 

0.64 

— 

4525 

0.33 

0.32 

tmmt 

*gST"W» 
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<1  <1*4-  <3— 

DIRECTION  OF  JET  MOTION 

(b) 

FIG.  I.  EXPLOSIVE  SHAPED  CHARGE  JET  FORMATION 
(port  (b)  from  R.  J.  Eichelberger  ond  E.M.  Pugh, 

J.  Appl.  Phys.  _23,  537  ( 1952  ). ) 


FIG.  3.  VACUUM  G.  MBER  USED  IN  FIRING  NO.  109 


LOCATION  OF 
EXPLOSIVE  ASSEMBLY 


FIG.  4.  PLAN  -VIEW  OF  FIRING  CHAMBER 


FIG.  5.  TYPICAL  TIMING  SEQUENCE 


G  6.  FRAMING  CAMERA  RECORD  OF  JET  FORMED  BY  COPPER- LINED  SHAPED  CHARGE 

( 1.4/xSEC  PER  FRAME) 


O*  COLLISION  SPECTRUM 


b  *Cu  SPARK 


FIG.  9.  TIME-INTEGRATED  SPECTRUM  OF  JET  COLLISION  IN  AIR 


TIME 


FIG  10.  FRAMING  CAMERA  RECORD  OF  JET  COLLISION  IN  FIRING  NO.  109 

( 1.4 SEC  PER  FRAME) 


TIME 


FIG.  II.  FRAMING  CAMERA  RECORD  OF  JET  COLLISION  IN  FIRING  NO.  124 

(1.4 /i-SEC  PER  FRAME) 


0*  COLLISION  SPECTRUM  t>*  Cu  SPARK 


FIG.  13.  TIME- INTEGRATED  SPECTRUM  OF  JET  COLLISION  IN  HELIUM 


■  .?  ■  - 


FIG.  14.  TIME-RESOLVED  SPECTRUM  OF  JET  COLLISION  IN  HELIUM 


FIG  15.  DENSITOMETER  SCANS  OF  TIME -RESOLVED  SPECTRUM 


WAVELENGTH  DEVIATION  it) 


FIG.  18. 


WAVELENGTH  DEVIATION  (A) 


Ai.ISN3i.Nl  3NH 


FIG  20.  Cul  5<53  INTEGRATED  LORENTZSAN  FtTS 


CttI 

*.T^SEC 


A1ISN31NI  3NM 


C»l  4029,4$ »,  4940 


FIG.  22  INTEGRATED  L0RENT2  IAN  PROFILE  Co  I  4525,  4531,  4540  AT  5.7/u  SEC.  4531  SPLIT 
AND  SHIFTED 


FIG.  23.  INTEGRATED  LORENTZI AN  PROFILE  Cu  I  4525, 4531,  4540  AT  S.lfi  SEC  4531  SPLIT 
AND  SHIFTEO 


FIG.  24  INTEGRATED  LORENTZI  AN  PROFILE  Cu  I  4525,  4531,  4540  AT  5.7  SEC.  NO  SHIFT 
IN  4531 


0*1  4925,4391,4943 


WWP4IW  '.w 

*  '  *'  ^  ■■^4  j 


FIG. 27.  VARIATION  OF  CUI  4531  LINE  SHIFT  WITH  PLASMA  ION  OENSITY 
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Summary  of  Experimental  Firings 
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Summary  of  Experimental  Firings 
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APPENDIX  B 


Direct-Vision  Prl*m  and  Streak  Camara  Combination 

A  discussion  of  the  direct-vision  prism  and  Beckman  and  Whitley 
streak  camera  used  in  combination  as  a  time-resolving  spectrograph  is 
contained  In  the  following  paper. 
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Time-Resolved  Spectroscopy  of  High 
Speed  Events* 

FtANX  R.  SCHWAtTI 

Explcsitu  Restore*  Sodtem,  Picaivtny  Arsmvd,  Dour :  Sew  Jersey 

AND 

Benjamin  J.  Pibnicx 

Simau  InstikmU  ef  T*ck*Mo[y.  tfoMtm,  New  Jersey 
(Received  April  20, 1962) 


THE  time-resolved  spectral  information  relating  to 
highly  transient,  self-luminous  phenomena  over 
selective  spectra]  ranges  in  the  visible  region  can  be  readily 
obtained  with  a  high  speed  streak  camera  by  the  simple 
addition  of  a  direct-view  prism  into  the  existing  entrance 
optics  of  the  camera  system,  litis  technique  has  been 
employed  here  in  observing  spectral  variations  with  re¬ 
spect  to  time  of  such  luminous  events  as  exploding  wires, 
hyperveiodty  impact,  and  detonating  high  explosives. 

The  optical  arrangement  utilized  is  shown  in  Fig.  !. 
Light  from  a  source  is  focused  on  an  adjustable  slit.  The 
slit  is  positioned  at  the  focal  point  of  a  7  in.,//2.3  loss, 
which  in  conjunction  with  the  objective  lens  of  a  Bcckman- 
Whitley  continuous  writing  streak  camera  forms  a  collimat¬ 
ing  optical  system.  A  Leybold  20X  20  mm  direct-view 
prism,  inserted  between  the  lenses,  disperses  the  light  from 
the  first  columa  ting  lens.  A  dove  prism  is  used  with  the 
direct-view  prism  to  facilitate  preliminary  adjustments  and 
alignment  so  as  to  obtain  the  desired  spectral  range  on  the 
film  record.  Inclusion  of  the  direct-view  and  dove  prisms 
reduces  the  optical  speed  of  the  original  system  from  // 5 
to  about // 25.  The  spectral  region  that  can  be  covered  with 


the  above  optical  arrangement  extends  from  about  MOO  to 
6500  A. 

Figure  2  is  a  film  record  that  illustrates  the  spectral 
variations  in  time  of  a  gold  wire  (99.99%  pure)  0.0015  in. 
in  diameter  and  0.040  in.  long,  when  exploded  by  a  4-kV 
pulse.  The  writing  speed  of  the  camera  is  1.0  mm/#*sec. 
Calibration  lines  from  helium  and  zinc  spectra]  lamps  are 
also  included  on  the  sarrie  film  strip.  The  wavelength  region 
shown  here  is  from  4400  to  5500  A,  which  lies  within  the 
spectral  sensitivity  range  of  the  Kodak  35-mm  Tri-X  film 
used.  The  maximum  spectral  resolution  within  this  range 
is  about  2  A. 

The  calibration  lines  on  the  film  strip  are  images  of  the 
full  entrance  slit.  Thus,  for  broad  sources  that  fully  illumi¬ 
nate  the  entrance  slit,  the  observed  slit  length  corresponds 
to  a  writing  time  or  time  response  of  1.8  jmec  at  1-mm/jcsec 
writing  speed1  (i.e.,  the  slit  image  will  completely  pass 
over  a  given  point  of  the  film  in  this  time  interval). 

For  the  situation  illustrated  in  Fig.  2,  however,  the  time 
response  is  much  smaller.  The  length  of  gold  wire  is  oriented 
perpendicular  to  the  entrance  slit ;  hence,  during  the  vari¬ 
ous  stages  of  wire  explosion  the  sht  is  not  completely  il- 
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NOTES 


Fio.  2.  Hme-reaoived  ipectrucn  of  exploding  gold  wire  with  He 
and  Zn  calibration  iinca.  (Writing  rate,  1.0  mm/itttc  ) 


luminated.  The  time  response  in  this  case  is  related  to  the 
diameter  of  the  expanding  luminous  region  about  the  ex¬ 
ploding  wire  rather  than  to  the  slit  length. 

For  applications  in  which  a  spectral  resolution  of  the 


order  of  an  angstrom  over  a  visible  range  of  the  order  of 
1000  A  is  of  interest,  the  direct-view  prism  method  for 
time-resolved  spectroscopy  is  advantageous.  The  prism  can 
be  readily  added  to  the  optical  train  of  most  high-speed 
cameras  without  the  use  of  additional  optical  equipment. 
Light  transmission  losses  through  the  prism  and  reduced 
optica!  speed  of  the  camera  system  have  presented  no 
difficulties  for  the  highly  luminous  sources  employed  here. 

This  work  was  undertaken  as  part  of  a  research  program 
in  explosive  study.  We  would  like  to  thank  Dr.  S.  J.  Lukasik 
of  Stevens  Institute  of  Technology  and  Dr.  S.  Koslov  of 
Vitro  Laboratories  for  their  continued  help  ami  guidance. 
We  are  also  grateful  to  personnel  at  Picatinny  Arsenal  for 
providing  the  necessary  facilities  and  permission  to  publish 
this  note. 

*  Work  supported  by  the  Ordnance  Corps,  U.  S.  Army 

1  G.  E.  Seay,  L.  B:  Seely,  Jr.,  and  R.  G.  Fowler,  J.  App).  Phys.  32, 
2439  {19615. 
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APPENDIX  C 

Spectral  Una  Profile  Program 

A  line  profile  I*  generated  from  dens ftometer  measurements  with  the 
use  of  the  spectral  line  profile  program  listed  below.  In  brief,  the 
program  converts  horizontal  and  vertical  position  measurements  on  a  dens¬ 
itometer  scan  Into  wavelength  and  Intensity  values  respectively.  The  line 
profile,  given  by  eg  (2-4)  is  then  constructed  by  subtracting  the  back¬ 
ground  Intensity  from  the  total  observed  intensity  at  each  wavelength 
step. 

The  Input  quantity  TIME  labels  the  program  output  with  the  time 
after  first  Impact  light  at  which  the  film  record  was  scanned.  A  density 
scale  calibration  Is  established  from  a  quadratic  least  squares  fit  of 
known  density  steps  D(l),  and  the  corresponding  densitometer  displacements 
Y(l)  .  A  listing  of  the  subroutine  POLYFT  which  performs  the  least  squares 
fit  Is  also  Included  below.  As  a  check  on  the  reliability  of  the  fit, 
density  values  D£V(f),  are  computed  from  the  input  densitometer  displace¬ 
ments  Y(l),  and  compared  with  the  known  density  steps  D(l).  A  wavelength 
scale  Is  constructed  from  the  horizontal  position  measurements  POS(l),  of 
several  lines  on  the  densitometer  record  with  known  wavelength  WVLT(l), 
and  wavenumber  WVNO(I).  All  position  measurements  are  with  respect  to  an 
arbitrary  zero  reference  point  In  the  horizontal  direction.  The  subrou¬ 
tine  POLYFT  is  again  used  to  obtain  a  quadratic  least  squares  fit  between 
the  variables  POS(l)  and  WVNO(I).  Similarly  the  reliability  of  the  fit 
Is  checked  by  compering  wavenumbers  DIFF(l)  computed  from  the  Input  posi¬ 
tion  measurements  POS(l)  with  the  known  wavenumbers  WVNO(i). 

Each  spectra)  line  is  labeled  by  the  wavelength  value  of  the  line 
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center  LWAVE,  The  number  of  input  points  used  to  describe  the  Hne  pro¬ 
file  Is  ceiled  l  .  Three  variables  define  an  input  point;  namely  X(t) 
the  horizontal  position,  YC ( I )  the  vertical  position  of  the  background 
density,  and  YT { I )  the  vertical  position  of  the  total  density.  Since  the 
length  of  a  densitometer  record  exceeds  the  horizontal  scale  of  the  Gerber 
scanner,  the  scanning  device  must  be  translated  in  order  to  read  alt 
necessary  points  on  the  densitometer  record.  Thus  one  cannot  directly 
obtain  the  required  continuous  sequence  of  values  for  the  variable  X(l), 
measured  from  the  arbitrary  zero  reference  point,  over  the  full  range  of 
the  densitometer  record.  This  difficulty  does  not  occur  for  y-positlon 
measurements  since  the  width  of  the  densitometer  record  is  easily  accom¬ 
modated  by  the  scanner. 

A  uniform  X-scale  is  obtained  in  the  following  manner.  The  current 
value  of  X(l)  Is  compared  with  the  previous  value,  i ,e.  X(l)-STORE.  If 
there  has  been  no  translation  of  the  Gerber  scanning  device  X { I )  will  be 
numerically  greater  than  the  quantity  STORE.  Alternately  if  X{l)  is  less 
than  STORE,  a  translation  has  been  made.  The  amount  by  which  the  scan¬ 
ning  device  has  moved,  called  SHIFT,  is  then  added  to  the  current  value 
and  all  subsequent  values  of  X < * ) ,  l.e.  X ( I )  «  X ( I )  +  SHIFT.  Any  number 
of  translations  of  the  Gerber  scanner  can  be  handled  in  this  fashion. 
However,  this  procedure  places  restrictions  on  the  manner  in  which  the 
Gerber  translations  are  made.  First,  one  must  ascertain  that  the  last 
X ( I )  value  read  before  translating  is  greater  than  the  next  X(l)  value  to 
be  read.  Also  the  amount  of  translation  is  not  arbitrary.  The  scanning 
device  Is  relocated  so  that  the  last  X-position  read  becomes  the  starting 
point  of  subsequent  readings. 
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After  converting  sM  X(l)  values  for  a  given  spectral  line  Into  a 
continuous  sequence  the  corresponding  wavenumbers  MAVEK(i),  «nd  wave¬ 
lengths  MAVEL(t)  are  computed.  The  spectral  line  Intensity  at  each  wave¬ 
length  Is  then  calculated.  The  total  density  is  found  from  YT(l),  and 
the  background  density  from  YC(I).  For  a  given  value  of  the  photometric 
slope  GAhttA(J),  the  line-  intensity  8R!Tl(f)  is  found  in  accord  with 
eg  (£M0  .  The  maximum  value  BRITM,  of  the  quantities  BAITL(l)  Is  then 
obtained  and  alt  BRiTL(l)  values  are  normalized  by  division  with  BRITM. 
The  wavelength  at  which  the  line  intensity  is  a  maximum  Is  called  WVLHX. 
The  output  of  the  program  consists  of  the  quantities  BRITL(l)  end  OEL(I), 
where  BftiTL(l)  Is  the  normalized  spectral  line  Intensity  and  DEL(i)  the 
wavelength  difference  from  the  position  of  maximum  intensity.  These 
latter  computations  are  then  repeated  for  the  K  values  of  GAMMA (d) . 
After  completion  of  these  spectral  line  profiles  for  a  given  line  the 
program  is  ready  to  similarly  process  the  next  set  of  densitometer  input 
points . 
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SPECTRAL  LINE  PROFILE  MK6 
USES  SUtoROUTlNE  POLVF 1 T 

DIMENSION  Dt  40).  Y<  40 > .  A<3>.  DEV<40)«  POSt  40 ). 

X  4 AVEKtltoO).  WMVtL < i 50 ) *  toRlTLtibOJ.  U13). 

X  DEL  (  I  SO ) *  WVLT<40).  WVN0<40>.  XtlSO)*  YCUSO)*  YTtlSO)  * 
X  GAMMA { 5 )  ♦  COEFFtS)  .D1FF<40) 

COMMON  Y  «  O  *  COEFF 

EOUI VALENCE  t  Y  »  POS  >  ♦  <D»  4VNO  ) 

eou i valence  t  dev*  oiff.  w«ven.  del  > 
equivalence  t  X  »  WAVEL  )  .  (  YC.  WRITE  > 

1  READ  to6«  TIME 
PUNCH  to6.  TIME 
PUNCH  96 

READ  39.  DUMMY*  N*  M 
DO  2  1=1*  N 

2  READ  98  «  D<  1  >  ♦  Y  <  I  ) 

CALL  POLYFT  < Y « D . 2 *N • COEFF ) 

At  1  >  *  COEFF  t  1  ) 

A  t  2  >  *  COEFF  t 2  > 

At  3)  ^  COEFF <J) 

DO  5  l  =  1  »  N 

DEV(I)  *  Adi  +  A<2>  *Y«I>  4A(3)*Y( 1 >*Y< I J 

DEV  t 1  )  *  0(1)  -  DEV t I  > 

5  PUNCH  98.  U(I).  Y  (  1  )  «  DEV (  1  ) 

PUNCH  95*  All)*  A ( 2 ) «  A  <  3  > 

PUNCH  96 

DO  6  I  -  1  *  M 

6  READ  93*  POStl).  WVLT(I)*  WVNOt I ) 

CALL  POLYFT  tPOS.WVNO.2 *M* COEFF > 
tot  1  )  *  COEFF  t I  ) 

to ( 2  >  «  COEFF (2) 
tot  3 )  *  COEFF (3) 

OO  l 1  I  *  I •  M 

D  I FF  (  I  >  =to  <  I  )  4-  e<2>*POStl>  -f  to<3)*POS< I )*PU3t I > 

D  IFF  <  1  )  *  WVNOti)  -  OlFFCl) 

H  PUNCH  93*  POStl)*  Ib'VLT  (  1  )  *  UVNQt]  )«  DlFFW) 

PUNCH  95.  tot  1  >  «  to ( 2 ) .  to  t  3 ) 

PUNCH  96 
STORE  =  0.0 
SHIFT  *  0.0 
READ  69 «  DUMMI  «  K 

21  READ  90.  LWAVE.  E 
PUNCH  90.  Ltt«VE 
PUNCH  96 

OO  22  1=1.  L 

22  READ  to6 «  X(I).  YCtI).  YT  t  1  ) 

DO  26  I =1 »L 

irt  Xtl)  -  STORt  J  24.  24.  23 
24  SHIFT  =  SHIFT  4  STORE 

23  STORE  =  Xtl) 

Xtl)  =  X< I )  4  SHIFT 

W  AVEK  <  I  I  =  to (  1  )  4  to ( 2  )  *X  <  1  )  4  tot  3) »X<  l >*X<  I  ) 

26  wAVEL  <  I  *  *  i.O/wAVEX(I) 

PUNCH  95.  STORE  .  SHIFT 
PUNCH  96 
OO  39  J  *  1  *  x 

READ  39  *  GAMMA ( J ) 

PUNCH  69.  OMMMm  t  J  ) 

PUNCH  96 

OO  1 2  I  =  1 •  L 

12  toRITLt  I  )=EXPFt  (At l  )4Mt2>*YTt  1  )4m(3)*YT(  I  )*YT<  1  )  ) / 


X  (GAMMA! J)*u*434£9l ) 

X  EXPF4 4A4 l  M-A4£)*YC< 1  >+A<3>*YC4  l >*YC4 1 >  )  / 

X  4  GAMMA ( J )  #0 *43429 ) ) 
bftlTM  *  HR  I  TL  4  1  ) 
ttVLMX  *  MAVEL4 1 > 

DO  42  1  -  £«  L 

IF  4  BR1TM  -  ttWITLdJ  )  41*  42#  42 

41  bRITM  ■  fclRiTLU  > 

WVLMX  a  WAVLL41) 

42  CONTINUE 

DO  33  I  »  I.  L 
ttRITLdl  »  b«lTL(  I  »/dWI  TM 
DEL  4  I )*WAVEL4 I  J-WVLMX 
33  PUNCH  SS*  wMVtLil).  DEL  < 1  > ♦  6R1TL41) 

PUNCH  96 
39  PUNCH  9b 
GO  TO  £1 

4*6  FORMAT  4  3F9*3> 

69  FORMAT  4  HU.b.  £15  > 

98  FORMAT  4  F8#£*  FU.Ji  Fd*3  > 

96  FORMAT  4/) 

95  FORMAT  4  3E£Q#7  > 

93  FORMAT  4  FI  0*3*  3F12«3  > 

65  FORMAT  4  2620.4 «  fc£0*6  ) 

90  FORMAT  (  114#  15  ) 

END 


SUBROUTINE  POL YF  T  C  X  *  Y  ♦  K  *  N  *  C  ) 

C  PUT  OAT*  POINTS  AND  FINAL  COEFFICIENTS  IN  COMMON 

DIMENSION  XC  40  )  ♦  Y  (  40>*B<S>*S<4)*AC5«b)»C<5> 

COMMON  X  *  Y  *  C 
CN*N 
XI  *K+  1 
DO  2  0*2* K.I 
02*0-1 
DO  I  1*1*02 

1  Af  O*  I )«O»0 
AC 02* O) *0*0 

2  A( J* J )* i *D 
A<  1  t  1  )  *  1  *0 
tit  I  i«0*o 
DO  3  1*1 fN 

3  b< 1 >sttC 1 )4Y< I ) 
b( I >*bC i i/CN 
ALF=0*0 

DO  4  I • 1  *  N 

4  ALF*ALF+XC1) 

A(2« 1  )  s-ALF/CN 
0*0*0 

SC  I  >  *0*0 
0(2 ) *0*0 
DO  5  1*1 *N 
P*A(2 . 1 >+XC I ) 

0=0+X( l )*P*P 
SC  1  >  =  SC  1  )+P*P 

5  bC2)*bC2)+YC 1 )*P 
b<2)*bt2)/St 1 > 

ALF«0/SC l ) 
bETA*SC 1 )/CN 

I F  C K— 2 ) 1 1 *0*0 

6  DO  lO 
M1*M+1 
M2*M-1 

A  C  M I • 1 ) s-ALF* A  C  M  » 1 ) -BETA* A  t  M2  *  1  ) 

DO  7  0*2 «M 
02*0-1 

7  AC  Ml t0)*ACM«O2l-ALF*ACM» 0>-bETA*ACM2*0) 

SCM)*0*O 

0*0*0 

R*0*0 

DO  9  1*1  *N 
P*  1  *0 

DO  a  L* 1 *M 
NL  =M  1  “L 

b  P=P*X  C  i  J-fACMl  *NL) 

U*0+XC  J  >*P*P 
H*R+Y( I ) »P 
9  SC M) *SCM)+P*P 
ALF  SO/'S  CM) 
bETAsSCM) /SCM2 ) 

10  DC  Ml  ) *R/SCM ) 

11  UO  12  L=1 «X1 
CCDsO.U 

DC  12  L1*L«K1 

12  CCl  )*CCL)+t>CLl  )*ACL1  *L  ) 

Rfc  TURN 

END 
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APPENDIX  D 


Photometric  Emulsion  Calibration 


A  rotating  step-sector  wheel  method  was  used  to  obtain  the  photo¬ 
metric  emulsion  calibration  of  the  time-resolved  spectral  film  records. 
The  experimental  arrangement  employed  is  shown  schematically  In  Fig.  D-l. 
Light  from  a  spectral  lamp  source  Interrupted  by  a  rotating  step-sector 
wheel  Is  focused  on  a  rectangular  si  It.  To  provide  for  uniform  Illumina¬ 
tion  over  the  silt  the  magnification  Is  such  that  only  a  small  central 
region  of  the  extended  light  source  Is  projected.  The  silt  is  viewed  by 
the  direct  vision  prism-streak  camera  optical  train  (see  Appendix  B)  with 
the  narrow  width  spectrograph  entrance  slit  removed.  To  vary  the  total 
exposure  time  of  a  calibration  run  a  manually-operated  shutter  is 
inserted  In  the  optical  system  as  indicated  in  the  figure.  The  streak 
camera  mirror  Is  stationary  during  the  calibration  exposure  In  order  to 
obtain  static  images  of  the  silt. 

The  step-wise  variations  In  exposure  along  the  length  of  an  image 
are  determined  by  the  rotating  sector  wheel  geometry  and  the  total  expo¬ 
sure  time.  From  eq  (2-10a)  the  film  exposure  at  a  given  point  on  the 
image  can  be  expressed  as 


i 


J(v')  dv'  dt 


t  A 


(D-l) 


where  J(v')  is  the  profile  of  a  spectral  line  from  the  source  lamp, 

A  is  the  slit  Image  width  (see  Appendix  E) ,  and  t  is  an  effective 
exposure  time,  t  is  determined  by  the  time  per  revolution  of  the  sector 
wheel  during  which  a  portion  of  the  slit  is  uncovered,  and  the  length  of 
time  the  mechanical  shutter  is  kept  open.  If  t 


is  the  time  per 


02 


revolution  of  the  sector  wheel  et  a  given  step,  F  the  rotational  fre¬ 
quency,  and  t  the  total  exposure  time  of  the  manual  shutter,  then  the 

effective  exposure  time  is  simply  the  product  t  Ft 

$  £ 

Since  the  source  light  is  steady  in  time,  one  has 

#  -  tsFte  J  J(v  ')  dv*  (0-2) 

£ 

if  the  slit  Image  width  is  very  large  in  comparison  with  the  width  of  the 
spectral  fine  then  the  film  response  at  each  step  exposure  is  given  by 

#  *  I  Ftste  (0-3) 

where  *  represents  the  tote!  line  intensity,  i.e. 

m 

I  -  j  J(v')  dv' 
o 

The  total  line  intensity  and  rotational  frequency  are  constant  for  a 
given  calibration  run.  For  the  sector  wheel  used  in  the  calibrations, 
successive  values  of  tg  are  in  the  ratio  2  to  1 .  The  manually  con¬ 
trolled  exposure  times  were  chosen  so  as  to  cover  a  broad  range  of  dens¬ 
ity  values  on  the  film  record,  from  about  0.05  to  1.7  density  units. 

The  measured  variations  of  film  density  with  relative  exposure  for 
a  zinc  spectral  line  are  illustrated  in  Fig.  D-2.  Density  values  were 
obtained  from  densitometer  readings  of  the  slit  image  and  relative  expo¬ 
sure  values  derived  from  eq  (0-3).  Measurements  corresponding  to  differ¬ 
ent  total  exposure  times  t  ,  are  distinguished  in  the  figure.  Over  a 
density  range  from  about  0.1  to  1.2  the  density  values  vary  linearly  with 
the  logarithm  of  the  relative  exposure,  i.e. 

0  ■  Y  log  ♦  +  const  (0-4) 


I 


Th*  slope  Y  U  determined  front  a  least  squares  fit  to  the  data  points. 
A  summary  of  the  values  of  Y  obtained  in  this  manner  is  presented  !n 
Table  V.  Tie  density  values  for  the  observed  spectral  lines  were  within 
this  I  inear  density  range. 
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FIG  D-l.  PHOTOMETRIC  EMULSION  CALIBRATION  SETUP 
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FIG.  0-2.  PHOTOMETRIC  EMULSION  CALIBRATION 


APPEND iX  E 


Spectrograph  Slit  Image  Width  Calculations 

Numerical  estimates  of  the  width  of  the  spectrograph  entrance 
slit  image  were  obtained  from  an  optical  ray  trace  analysis  of  the 
direct  vision  prism  streak  camera  combination.  Principal  rays  origi¬ 
nating  from  the  entrance  slit  are  deviated  by  the  prism  and  focused 
at  the  image  plane  of  the  streak  camera*  thereby  defining  the  slit 
image  width.  Since  the  deviations  of  these  principal  rays  are  wave¬ 
length  dependent,  the  variation  of  slit  image  width  with  wavelength 
Is  calculated  over  the  spectral  range  encompassed  by  the  time-resolved 
spectral  records. 

Consider  the  ray  diagram  illustrated  in  Fig  E-l.  A  ray  from 
the  center  and  from  the  edge  of  the  entrance  slit  are  followed  through 
the  collimator  lens,  direct  vision  prism,  and  telescope  lens  of  the 
streak  camera.  The  total  width  of  the  entrance  slit  is  called  S  , 
the  focal  lengths  of  the  collimator  and  telescope  lens  are  fj,fj 
respectively,  and  the  total  slit  image  width  is  A  .  The  prism  angles 
VVE  *n<*  *nt«r^ac®  angles  of  importance  indicated  by  other  greek 
letter  symbols  ere  also  given  in  the  figure.  The  subscript  1  refers 
to  interface  angles  associated  with  the  ray  from  the  entrance  slit 
center,  and  similarly  the  subscript  2  is  applied  to  angles  associated 
with  the  ray  from  the  entrance  slit  edge.  The  index  of  refraction  of 
the  identical  outer  crown  glass  prisms  is  called  ,  and  for  the 
middle  flint  glass  prism,  n^  .  At  the  wavelength  for  <ero  deviation 


I 


the  silt  image  is  centered  about  the  point  located  at  O'  fn  the 
image  plane.  For  other  wave  lengths  the  sHt  image  is  located  at  a 
typical  point  x  as  illustrated.  The  angular  deviation  of  the  center 
ray  from  the  principal  axis  is  denoted  by  T  . 

The  following  angular  relations  can  be  readily  derived 


tan  9j 

-  S/2f? 

e,  -  x,-x2 

sin  Xj 

■  nc  s in  C  j 

sin  X2 

*  nc  s  ' n  ^2 

e,  * 

m  a 

c 

c2* 

m  a 

c 

sin  T]  ^ 

-  (np/nc>  X 

sin 

Sin  T}2 

■  (r>p/nc)  * 

sin 

52 

!i  *  *i 

52  *  «2 

-  “f 

<nF/nc) 

x  sin  e  j  m 

s  i  n 

<"F/nc> 

x  sin  * 

s  i  n 

*2 

«,  .V, 

m  a 

c 

*2  +  V2 

m  a 

c 

nc  sin  Y,  -  sin£j 
nc  sin  Yj  •  sin  3^ 


-  Z  +  ^  |  - 1 90° 
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The  position  of  the  slit  Image  center  Is  given  by 

x  *  f 2  ten  f 

The  edge  of  the  silt  image  It  «t  a  position 

x  ♦  A/2  »  fj  ten  ff+Gj) 

Hence  the  slit  width  A  ,  Is  given  by 

A  *  2  fj  (ten  (T+G^)  *  tan  T) 

Numerical  solutions  of  the  above  set  of  equations  were  performed 
with  the  use  of  the  Slit  Function  Estimate  program  given  below.  Vari¬ 
ables  in  the  program  are  associated  with  terms  In  the  above  equations 
as  follows. 


FI.F2 

focal  lengths  f^.f^ 

S 

entrance  slit  width,  S 

ALPC.ALPF,  SIGMA 

prism  angles,  <*c,  cr^,  £ 

XI 

angle  of  Incidence,  x( 

NEFC,  NEFF 

indices  of  refraction,  n 

i 

HAVEL 

wavelength 

2 

C 

ETA 

Tl 

SE 

? 

EPS 

t 

DELTA 

6 

GAMMA 

Y 

BETA 

a 

THET1 

•> 

X2 

*2 

E4 


PSI 

¥ 

THET2 

PHI 

¥  + 

82 

OISP 

X 

DtSPL 

x 

i/2 

SLIT 

6/ 2 

The  program  systematical 1y  evaluates  the  angles  from  the  above 
equations,  First  for  the  principal  ray  from  the  entrance  slit  center 
(f.e.  BETA!,  DELTI,  etc.),  and  then  for  the  principal  ray  from  the 
entrance  slit  edge  (l.e.  BETA2,  0ELT2,  etc.).  The  angular  deviations 
¥  and  ¥  +  0^  are  then  computed.  The  silt  Image  center  position 
OISP,  the  edge  position  D1SPL,  and  finally  the  slit  width  SLIT  at  the 
wavelength  WAVEL  are  obtained  as  part  of  the  program  output.  These 
calculations  are  repeated  at  wavelengths  for  which  refractive  index 
data  are  available. 

A  dispersion  curve  Is  constructed  from  the  above  calculated 
OISP  and  WAVEL  values  with  the  use  of  the  Least  Squares  Fitting  of 
Slit  Function  program  and  LSOF  subroutine  listed  below.  This  is  a 
necessary  step  since  the  calculated  SLIT  values  are  in  units  of  length 
In  the  image  plane  and  must  be  converted  to  wavelength  units  bv  means 
of  a  dispersion  curve.  Also  the  calculated  dispersion  curve  can  be 
compared  with  an  experimental  calibration  to  provide  a  check  on  the 
ray-tracing  calculations.  The  above  calculated  values  of  DISP,  SLIT, 
and  WAVEL  are  used  as  Input  information  to  these  fitting  programs  and 
appear  as  dimensioned  variables.  For  convenient  numerical  calculations 
a  least  squares  polynomial  fit  is  used  to  represent  the  functional 
dependence  between  wavelength  and  position  in  the  image  plane.  A 


least  squares  fit  of  order  M  of  the  quantities  FRQ{|)  and  0 1 S P ( I ) 
is  obtained,  where  FRQ(I)  is  the  wavenumber  equivalent  of  WAVEL(f). 

9y  comparing  deviations  in  calculated  and  input  wavelengths,  i.e. 

OIF  •  WL  -  WAVEL(I),  for  various  values  of  the  order  of  the  fit  M  , 
a  best  value  of  M  Is  obtained,  A  vaiue  of  M  equal  to  4  was  found 
to  yield  the  best  fit  polynomial  in  the  numerical  computations. 

The  silt  Image  width  in  wavelength  units  is  then  obtained  with 
the  us#  of  this  fourth-order  polynomial  representat ion  of  the  dis¬ 
persion  curve.  The  position  of  the  center  of  the  slit  Image  OISP(I) 
has  a  calculated  equivalent  wavelength  vaiue  WL  .  Similarly  the 
position  of  the  edge  of  the  slit  image  DiSPL  has  an  equivalent  wave¬ 
length  value  WPL  .  The  difference  WPL-  WL  gives  the  slit  image 
width  in  wavelength  units  DELW. 

Numerical  estimates  of  the  slit  image  width  were  computed 
using  refractive  index  data  supplied  by  the  prism  manufacturer,  meas¬ 
ured  prism  angles,  and  known  collimator  and  telescope  lens  focal 
lengths.  These  data  are  presented  In  Table  E-l.  The  entrance  slit 
width  $  and  the  angle  of  incidence  Xj  were  treated  as  parameters 
and  varied  about  their  nominally  known  values.  The  silt  Image  width 
was  calculated  at  those  wavelengths  for  which  refractive  indices  were 
known.  Results  of  the  computations  are  given  in  Table  E-ll  and  *llus 
trated  in  Fig.  17.  Angular  deviations  of  11°  about  a  nominal  32.5° 
incidence  angle  were  considered  representative  of  the  uncertainty  in 
the  alignment  of  the  direct  vision  prism  in  the  optical  train.  The 
slit  image  width  was  computed  for  several  values  of  the  entrance  slit 
width  at  each  value  of  the  incidence  angle.  The  calculated  position 


of  tha  canter  of  the  slit  image  y  ,  et  eech  wavelength  is  also  in¬ 
dicated  in  the  table. 

Calculated  wavelength  positions  x  ,  were  compared  with  meas¬ 
ured  positions  obtained  from  a  densitometer  trace  of  the  copper  spark 
spectrum,  originally  used  for  wavelength  identification  on  the  time- 
resolved  spectral  record.  This  would  provide  a  check  on  the  ray- 
tracing  computations.  The  calculated  position  values  given  in  Table  E-i! 
were  transformed  to  the  position  scale  of  the  densitometer  trace.  A 
linear  transformation  was  used  since  the  densitometer  scan  represents 
a  magnification  of  the  original  film  record,  A  dispersion  curve  was 
then  constructed  from  a  least  squares  fit  using  these  adjusted  position 
values.  Table  E-Ill  gives  the  calculated  position  for  various  wave¬ 
lengths  in  the  4500-5300A  range  obtained  from  the  dispersion  curve, 
and  the  corresponding  position  from  densitometer  trace  measurements. 
These  results  are  seen  to  be  in  close  agreement.  The  lack  of  system¬ 
atic  variation  between  the  calculated  and  measured  position  values  is 
to  be  noted.  The  differences  given  in  Table  E-Ill  appear  to  oscillate 
about  a  zero  mean.  This  is  further  indication  that  the  approximations 
and  numerical  estimates  employed  in  the  calculations  are  realistic  and 
that  the  slit  image  width  estimates  are  reasonably  good. 


TABLE  E-! 


Refractive  Index  Values 


s  lengt'h 
(A) 

Crown 

Flint 

4047 

1  .53022 

1.65071 

4358 

1  .52667 

1 .64206 

4800 

1.52282 

1.63312 

4861 

1.52236 

1.63210 

5461 

1.51871 

1.62410 

5876 

1.51680 

1.62004 

5893 

1.51671 

1.61992 

6438 

1-51471 

1-61583 

6563 

1.51431 

1.61504 

Prism  Angles  (deg) 


a 

c 


98.5 
132 

57.5 


Focal  Lengths 
f| 


(in.) 

7 

12 


/ABLE  E- 1 1 


Variation  of  Silt  Image  Width  with  Wavelength 


Angle  of 
Incidence 

Wavelength 

Image 

Pos  1 1 1  on 

Total  Silt  Width 

A  (mm) 

5  «  0,07 5mm  S  *  0.100mm 

*1 

X 

(deg) 

(A) 

(mm) 

(A) 

(A) 

31.5 

404? 

-37-9 

2.1 

2.8 

4358 

-20.2 

3-3 

4.4 

4800 

~  5-5 

5.0 

6.6 

4861 

-  4.0 

5-2 

7.0 

5461 

7.1 

8.5 

H.4 

5876 

12.2 

11.5 

152 

5893 

12.2 

n.5 

15.2 

6438 

17.0 

55.9 

21 . 1 

6563 

17.8 

16.7 

22.3 

32.5 

4047 

-45.6 

2. 1 

2.7 

4358 

-26.4 

3.5 

4.6 

4800 

-n.o 

5-2 

6.8 

486! 

~  9.4 

5.4 

7-2 

5461 

1.9 

8.8 

n.  7 

5876 

7.2 

1 1 .8 

15.7 

5893 

7.2 

1 1 .9 

15-9 

6438 

12,1 

16. 1 

21.5 

6563 

12.9 

17.3 

23.0 

335 

4047 

-54.7 

2.G 

2.7 

4358 

-33.3 

3.6 

4.8 

4800 

-16.9 

5.3 

71 

4861 

-15-3 

5-6 

7,4 

546 1 

-  3-5 

9-0 

12.0 

5876 

1-9 

12.  1 

16.2 

5893 

!  9 

12.  1 

16.2 

6438 

6.9 

16.8 

22.4 

6563 

7.8 

17.7 

23.7 

TABLE  E-t!l 


Calculated  and  Measured  Wavelength  Position  Comparison 


Angle  of 
incidence 

Wavelength 

Pos  i  t  ion* 

*1 

Calc. 

.  Meas. 

Difference 

(deg) 

(A) 

(arb. 

units) 

(arb.  units) 

31.5 

4500 

5.46 

5.30 

0.16 

4690 

7.73 

7.85 

-0.12 

4700 

9.78 

9-90 

-0.12 

4800 

11.62 

11.65 

-0.03 

4900 

13.27 

13.25 

0.02 

5000 

14.73 

14.65 

0.08 

5100 

16.02 

15-95 

0.07 

5200 

17.15 

17.10 

0.05 

5300 

18.14 

18.23 

-0.09 

32.5 

4500 

5.42 

5.30 

0.12 

4600 

7.72 

7.85 

-0. 13 

4700 

9.79 

990 

-0.11 

4800 

11.65 

11.65 

0.00 

4900 

13.29 

13.25 

0.04 

5000 

14.75 

14.65 

0  10 

5100 

16.02 

15.95 

-C.07 

5200 

17. 14 

17. 10 

0.04 

5300 

18.10 

18.23 

-0. 13 

33.5 

4500 

5.37 

5.30 

0.07 

4600 

7.72 

7.85 

-0.13 

4700 

9.8! 

990 

-0.09 

4800 

11.68 

11.65 

0.03 

4900 

13.33 

13.25 

v.07 

5000 

14.77 

14.65 

0.12 

5100 

16.03 

15.95 

0.06 

5200 

17. 12 

17.10 

0.01 

5300 

18.05 

18.23 

-0.18 

Pos 1 1 i on  on  a  densitometer  trace. 
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FIG.E-1.  OPTICAL  RAY  TRACE  DIAGRAM 
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c  slit  FUNCTION  ESTIMATE  mk  3 

Q  mWCSIN  FUNCTION  1NCLUOEU 

C  F|*  F3*  S  To  t>E  IN  CONSISTENT  UNITS 

c  Input  ancles  in  uegre.es 

AS  INF  (  X  )  -  1  *S7c79to3-bwRTF  <  1  • U“X  >  *  <  1 • 870  7t_8t3~ 

X  U*£I31  l44*X+U,U7436iU*X*X-U*Ol67393*X*X*X) 
l  READ  99  *  FI*  Ft 
READ  99 *S 

REAU  98*  ALPS*  mLPF  *  SIGMA 
WEAD  98*  Xl 

WEAD  97*  RLFC*  WEFF.  WAVEL 
1=0 

C  DEGREE  TO  WAD IAN  CONVERSION 

ALPC  -  ALPS*  I  .  74 
ALPF  »  ALPF*  1 • 748339E-U3 
SIuMa  *  sIgMm*  1  *  74‘jJtVt-Ot: 

XI  x  xl*  l • ?4Ss39E~U3 
X  a  Xl 

3  IF  <1*8707903  -  X)  7*4.4 

4  slNZ  =  S|NF(X)/REFL 
Z  a  AS  INF  (  SIN/  ) 

GO  TO  1 U 
7  X=  X  —  1  *87u7963 

SINZ  =  COSFtXIZREFC 
Z  a  ASINF(SINZ) 

X  =  X+  1  «b7U’Nbo 
lO  t-TA  =  ALPC  -  Z 

IF  (1.8707983  -  ETA)  18*  1  Id .  lid 

13  SlNSE  =  I«EFC*slNF (ET A ) JZWEFF 
SE  =  ASiNFt SlNSE) 

GO  TO  18 

lb  ETA  s  ETA  -  1 *8707903 

16  SINSE  =  (REFC*COSF  (ETaM/REFF 

17  SE  =  AS  1 NF  <  S I NSt ) 
t-TA  =  ETm  +  l  •  8707963 

IB  EPS  s  ALPF  -  SE 

IF  ( 1 •87G7963-EPS )  33*30.30 

30  SINUL  *  CREFF  *  S I NF < EPS ) ) ZREFC 
DELTA  =  AS  INF  l SINUL) 

GO  TO  36 

33  EPS  =  EPS  -  1*8707963 

SINUL  x  <REFF*COSF <t.PS> )/ReFC 

delta  =  asinf(sinul> 

EPS  =  EPS  +  1*8707963 
t6  GAMMA  a  ALPC  -  UELTA 

IF  (1.870796C  -  GAMMA)  31*38*36 
38  SlNb  a  REF  C  *  SiNF( GAMMA) 
ticTA  =  AS  I  NF  (SINS) 

GO  TO  33 

31  GAMMA  =  GAMMA  -  1*8707963 
SlNU  =  ReF C*COSF ( oAMMA ) 

OETA  =  ms  I  NF  (  s  1  NS  / 
vjAMMm  =  gaMMa  +  1  *8707963 

3J  IF  t  1  —  1  )  34*39*39 

34  1  =  14-1 

at  Ta  l  =  tsE  T  a 
DEL?  i  x  OeLTm 
EPS  1  a  EPS 
SE  1  =  SE 

E  TA 1  =  ETA 

Z1  =  Z 


Xl  *  X 

THtTl  =  ATANJ-  < S/(2.U*F I ) ) 
xi  a  xi  -  THtTl 
X»  X2 
uo  TO  3 

39  tofcTAg  *  ttfcTA 
OELT2  =  DELTA 
tPS£  =  EPS 
Sfc2  *  SE 
ETAS  a  fcT  A 
Z2  *  Z 
XS  *  X 

PSI  *  SIGMA  +  UtTAl  -  1.5707963 
THfcT2  =  CJt TmS  -  ufcTMl 
PHI  «  Pbl  +  TMtTS 
IF  (PSI  -  0«O)  43  * 46 ♦  46 
43  PS  I  =  -PSI 

OISP  «  -<F2*S1NF  (PSD  )/'COSF<PSI  > 
PSI =-P5l 
GO  TO  Si 

46  IF  (1*5707963  -  PSI)  49.47.47 

47  OISP  a  (F2*S1NF  C PS i ) ) ZCOSP < PS I ) 

SO  TO  Si 

49  PSI  *  Pol  -  1.5707963 

D1SP  a  -(F2  *  LOSF ( PS i ) ) XS 1 NF ( PS I ) 
PSI  a  PSI  +  1 *5707963 

51  IF  (PHI  —  O  •  U  )  52.55.55 

52  PHI  a  -PHI 

DISPL  »  -(F2*SINF(PH1 ) )/COSF(PHl ) 
PHI a-PHl 
GO  TO  60 

2S  IF  <1*5707963  -  PHI)  56.50. S6 
56  DISPL  »  (F2»SlNF(PHl > IXCOSF (PHI ) 

GO  TO  60 

56  PHI  a  PHI  -  1*5707963 

DISPL  =  -<F2*COSF<PHi ) J/SlNFCPHi > 
PHI  a  PHI  +  1*6707963 
60  SLIT  *  DISPL  -  DISP 

C  RADIAN  TO  DtGHtE  CONVERSION 

R  =  57*295766 
ALPC  «  ALPC*  R 

ALPF  =  ALPF*  « 

SIGMA  a  SIGMa*  R 
XI  *  X1*R 


PUNCH 

99. 

FI  .  F2 

PUNCH 

96. 

ALPC. 

ALPF .  S I 

GM  A 

PUNCH 

97. 

RtFC. 

RtFF 

PUNCH 

96 

PUNCH 

99* 

S .  XI. 

WAVtL* 

SLIT 

PUNCH 

96 

l F  ( sense 

Sw 1 TCH 

1  )  63 . 

1 

63  6ETAI  *  ot  TmI *R 
DELT1  a  UtLTI*R 
EPS1  a  LP51*R 
SE1  *  SE 1 *R 
fcTAl  »  ET A  I *R 
Z1  a  Zl*R 

tat  TA2  *  OE  Ta2*W 
DELT2  «  DELT2*R 
EPS2  *  tPS2*R 

oC.2  *  st«;*w 


i' 


LTA2  »  fcTA2#R 
ZZ  »  Z2*R 
XZ  «  X2*R 
PS  I  *  P6I*« 

PHI  a  PH 1*R 

PUNCH  9 6.  bLTAi  *  btT M2 
PUNCH  96.  UtLT 1 «  UhL  T2 
Punch  96.  t.psi  * 

PUNCH  96.  Sfcl.  Stc 
PUNCH  96*  LTAl*  uTAZ 
PUNCH  96 ♦  Z 1.  22 
PUNCH  98.  XI.  X2 
PUNCH  96.  PS  1  «  PHI 
PUNCH  99 .  0 l SP  »  U 1 SPL 
PUNCH  96 

PUNCH  96.  WAVLU.Ol SP.6C1 T 
GO  TO  1 

99  FORMAT  (4F9.3) 

98  FORMAT  <3F7*2> 

97  FORMAT  <2F6.S.  F9.JJ 
96  FORMAT  </> 

9b  FORMAT  (OFIO.U) 


6  NO 


n  r 


Lfc,AST  GuUAML  FITTING  OF  GL  l  T  FUNCTION  MKi 
INPUT  DATA  FROM  GLlT  FUNCTION  LGTjMATfc  MNt! 
DIMENSION  WAVLC<9).  U  1  GP  <9J  •  OL  1  T  <  9  )  *FRU(9)«  M<b> 

1  READ  98*  M 
DO  Z  I  si .9 

R£AO  99*  WAVfc.L<l>  D1GP<1).  GL  1  T  <  I  > 

2  FRG<  1  )  *1  •O/'WAVfc.Lt  1  ) 

CALL  LSOFtUIGP.  FRO.  M«  9*  A  ) 

4  DO  20  1*1.9 
FQY=A( 1 > 

MM*M+1 

DO  5  J*2.  MM 

5  FUY  *  FUY  +  A  <  U  )  *0  1  GP  <!)«■*(  J~  1  ) 

WL=  1  ♦O/'FUY 

D1F=WL-WAVLL< l > 

XPL  =  QIGP< 1 )  +  GL I T l 1 ) 

FPL=A< 1 ) 

0O  10  0*2. MM 

lO  FPL*FPL+A< UI *XPL**< 0-1 ) 

WPL* 1 • 0/FPL 
0fc.L9*WPL-WL 
PUNCH  90* 1 

punch  99.  gl i t  < i  ) 

PUNCH  99.  ttL  •  HiAVtL  i  1  >  *U  IF 
20  PUNCH  99.  WPL.DLLW 
GO  TO  1 

99  FORMAT  13F1U.J) 

90  FORMAT  ( 1 4 ) 
tND 


SUBROUTINE  LSOF  <XtY,M«N.Al 

DIMENSION*  X  <  9  >  ♦  Y  l  9  )  ♦  t)<9)«  C(5)t  At5>*  U<5»5>*  V(5> 

K.  =  2*M+1 

DO  1  L* 1  «  K 

b(L>*0*0 

DO  1  1  a  1  ♦  N 

1  b(L)*b(L) +X< I > 

•  1 1  *M  +  1 

DO  2  L* 1 »  M 1 
Cl L > aO.O 
DO  Z  1*1  »N 

2  Cl  L  )  »C<  L  )  +  Y<  I  >*X  <  1  )•»*<  L-  1  ) 

DO  3  I  «  1  .  M 1 

DO  3 
K* 1+J-l 

3  Ut  I  * J)*blK> 

DO  4  I  *  1  » M 
11=1+1 

DO  4  J* 1 1 «M1 
VI U)=UlO. 1 )/U< 1 « 1 » 

Cl J)*ClU)-VlU)*Cl 1 ) 

DO  4  K* U  «M1 

4  U(  JtM=U(  J«K>-Vl  J)*U<  1  «X) 

Ml  Ml  >*CIM1 )/UlMl »M1 ) 

DO  5  1  *  1 « M 
K=M-| +1 

0*0  •  o 

DO  6  J=K.«M 

6  0=Q+A< J+j )*UIN»U+1 » 

5  A(X)*lt(K.J-U)/g(N«K) 

RETURN 

end 


BLANK  PAGE 


n 


APPENDIX  F 

integrated  Lorentzian  Profile  Computer  Programs 
(a)  Single  Component  Profile 

7'ie  computer  program  employed  to  calculate  an  integrated  Lorentzian 
profile  for  an  Isolated  spectral  line,  described  by  eqs  (2-12)  and  (2-13) 
Is  given  below.  Input  terms  defined  In  the  program  are  identified  as 
fol lows : 

HFWIT:  total  halfwidth,  6 

WAVEL:  wavelength  of  spectral  line, 

DELWV:  increment  In  wavelength 

SLITW:  spectrograph  entrance  slit,  A/2 
All  wavelengths  are  in  angstrom  units. 

The  quantities  HFWIT  and  SUTW  are  first  converted  to  wavenumber 
equivalents,  DEFRQ  and  OSLIT.  The  resonance  halfwldth  of  the  spectral 
line  HFBRT,  is  then  calculated  in  accord  with  eq  (2-13),  and  converted 
to  angstrom  units  as  BREAT.  The  line  profile  at  selected  wavelength 
Increments  is  then  computed.  DELAM(l)  and  DELNU(I)  are  the  wavelength 
and  wavenumber  steps  measured  from  the  spectra)  line  center  wavelength 
WAVEL,  at  which  the  intensities  RELIN(l)  are  calculated.  These  intensity 
values  are  normalized  to  the  peak  value  at  the  line  center  by  division 
with  the  quantity  CONST.  The  negative  wavelength  portion  of  the  line 
profile  is  first  obtained,  represented  by  the  var tables  DELAM(l)  and 
REllN(l)  In  the  program  output.  The  positive  wavelength  portion  of  the 
Una  profile  Is  then  similarly  calculated,  again  designated  by  the  same 
variables  DELAM(i)  and  RELIN(l)  in  the  program  output. 
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C  INTEGRATED  LOREN T Z 1  AN  PROF 1 EE  MAW 

C  SINGLE  COMPONENT  PROFILE 

DIMENSION  DEL  AM  <200) «  QfcENU  (200)*  Rf-_LIN  <200) 
EQUIVALENCE  JEENU*  REE  IN ) 

I  READ  99.  HFW.  r 
READ  99*  WAVEL 
READ  99*  DEEteV 
READ  99*  SElTW 
READ  96*  ETERM 
WAVE  *  WAVEE  +  O.S*HFWlT 

OEFRQ  *  (  I *0  /  WAVE  )  -  <  1*0  /  WAVEE  ) 

WEENT  »  WAVEE  +  SLlTw 

DSEIT  *  C  1*0/  WEENT  )  -  <  1 *0  /  WAVEE  ) 

HF0RT  *  <  DEFRO*DeF RQ )  -  < USE  I T*0SE 1 T  ) 

HFBRT  «  SORTF  <  HFbRT  > 

RECIP  *  HFbRT  +  i.O  2  WAVEE 
BREAT  a  WAVEE  -  <  1.0  /  RECIP) 

PUNCH  99*  WAVEE.  HFW I T 
PUNCH  97 

PUNCH  99.  SElTW.  BREAT 
PUNCH  97 

CONST  *2.0*ATANF (USE  1 T/HFBRT > 

STORE  *  0.0 
DO  27  I  »  I  *  ETERM 
DEEAM  < 1 )  ®  ST ORE—DEEWV 
STORE  =  DEEAM(l) 

DEEAM < I )  a  DEEAM I l >  +  WAVEE 

DEENU <  I  )  *  < l .O/QEEAM! 1 ) >-! 1 .O/WAVEL > 

IF  <DELNU(1)  -OSE1T)  22.  22.  24 
22  REEIN< l ) * ATANF < < DEENU < I ) +DSL I T ) /HFBRT ) 

X  ♦  ATANF  I  IDSL  IT  -  DEENU!  I  ) ) /HFbRT) 

GO  TO  25 

24  REE 1 N< 1 >  ■  ATANF  <( DEENUt I > +DSE I T >/HF ORT ) 

X  -ATANF  i< DEENUt It  -DSEI T ) /HFbRT ) 

25  RELINU)  *  REE  I  N  (  I  )  /CONST 
DEEAM! I )  a  DEEAM! I )  -  WAVEE 

27  PUNCH  9b*  DEEAM!  I  )«  REE  IN  <  I  ) 

PUNCH  97 
STORE  «  0.0 
DO  40  1»1.  ETER.M 

DEEAM! 1)  «  STORE  +  DEEWV 
STORE  *  DEEAM! I  ) 

DEEAM! I )  *  DEEAM  <I>  +  WAVEE 

DEENU! I )  ■  < 1 .U/DEEAM! I )>-( 1 .O/WAVEE) 

DEENU! I )  ■  — UEENU! I ) 

1 F  ( DEENU ! I ) -USL IT)  35.35.07 
35  REElNt  1  )*ATANF  <  <  DEENU !  1  )  +U5E  IT)  /HF tiRT  ) 

X  + ATANF  <  (USE  1 T-DtENU ( I  )  ) /HFBRT ) 

GO  TO  3b 

37  REE  IN!  1  ) -ATANF  (  (DEENU!  1  )+D3l_i  T  )  /HFBRT  ) 

X  -ATANF  < (DEENU! 1 >-QSLIT ) /HFBRT ) 

3b  RELIN(I)  =  REE  IN!  I ) /CONST 
DELAMU)  r  DEEAM!  1)  -WAVEL 
40  PUNCH  9fi  *  DEEAM ( I  ) .  REElNil) 

PUNCH  97 
GO  TO  1 

99  FORMAT  <  2F9.3  ) 

96  FORMAT  (  F9.3.  E2U.B  ) 

97  FORMAT  (/) 

c'6  FORMAT  <  15  ) 

END 


(b)  Multicomponent  Profile 


A  listing  of  the  computer  program  used  in  the  calculation  of  a 
mult {component  Integrated  Lorentzian  line  profile  described  by  eq  (2~W 


below.  The 

following  Input  variables  are  Identified 

SLIT: 

spectrograph  entrance  slit,  A/2 

DELW: 

wavelength  increment 

WAVE (J) : 

center  wavelength  of  J th  line, 

A(J): 

resonance  halfwidth  of  J th  line,  Aj/2 

BI{J): 

relative  Intensity  of  Jth  line,  J(Vj) 

All  wavelengths  are  In  angstrom  units. 

The  variables  SLIT,  WAVE(J),  and  A(J)  are  first  converted  to  wavi 
number  equivalents  as  DSLIT,  FREQ(j) ,  and  BRET(j)  respectively.  Indlv 
dual  contributions  to  the  total  line  profile  from  each  broadened  speed 
component  are  then  computed  at  specified  wavelength  intervals  about  th< 
center  wavelength  of  the  principal  spectral  Itne  WAVE(1).  Negative  an< 
positive  wavelength  increments  are  called  DELAM(l)  and  DELWV(l)  respeo 
tlvely.  No  distinction  is  made  for  the  corresponding  wavenumber  steps 
DELNU(l).  The  spectral  component  intensities  RELIN(J)  and  SOURC(J)  an 
added  to  give  the  total  line  intensity  SUM ( I )  and  AOD(l)  about  the  cem 
wavelength  WAVE(l).  These  results  are  normalized  to  the  maximum  total 
Intensity  value  PHIZO.  The  normalized  line  intensities  SUM(l)  and  ADO 
at  wavelength  intervals  DELAM(l)  and  OELWV(l)  comprise  the  program  out| 


non 


INTEGRATED  LORENTZIAN  PROF 1 LE  MK7 

MULT  1  COMPONENT  PROFILE 

PRINCIPAL  COMPONENT  CALLED  WAVE (  l  ) 

DIMENSION  WAVE  <  5  >  •  A<5>«  B  I  (  5  )  *  FREU<S>*  BRET<b>*  DELAM < 60 ) ♦ 
X  D£LNU(60)»  RELlNtS)*  SUM<60)«  DEL WV 1 60 )  »  SOU«C<b).  ADO<60) 
EQUIVALENCE  <  REL IN*  SOURC  ) 

1  READ  99,  SLIT 
READ  99  *  DELW 
READ  96 ♦  LTERM 
READ  96 •  N 
DO  3  J*  1 «N 
READ  99*  WAVE ( J ) 

READ  99*  A ( J ) 

READ  99*  OUJ) 

3  PUNCH  99*  WAVE (  J  )  •  A<J>*  til(J) 

PUNCH  95 
PUNCH  99.  SLIT 
PUNCH  95 

SLITW  »  WAVE  < I  )  ~  SLIT 

DSL IT  »  I  l.O/SLITW  i  -  U*U/WAVE<1»  > 

DO  6  U*1 *  N 

FREQ  C  U I  *  J.O/WAVEIOJ 

6  BRET(U)  *  I  *U/  (  WAVE  I  J  )  -  A<J>  )  -  FREQ  I  J  ) 

STORE  *  DELW 

SUMA  *  0*0 

LL  *  LTERM  +1 

DO  20  1  *  1  *  LL 

DELAM ( I )  a  STORE  -  DELW 

STORE  =  DELAM l I ) 

OELAMII)  *  DELAMd)  ♦  WAVE  d  > 

DO  18  J*  I.  N 

DELNUII)  *  < 1 • O/DEL AM <  1  )  )  -  FREQ(U) 

I F  <  OELNUIl)  -  DSL 1 T  )  11*  11*  12 

11  RELIN(J)  *  ATANf*  <  ( DELNU  <  I  )  +  DSL  IT)/  BRET(U)  >  + 

X  ATANF I  (DSL  IT  -  DELNU ( i  )  >/  BRET(U)  ) 

60  TO  13 

12  RELIN(J)  *  ATANF  <  <  DELNU  (  1  )  +  DSL  I T ) /  BRET(J)  >  - 
X  ATANF  (  (DELNU(I)  -  DSL  I  T )  /  BRET(J)  ) 

13  RELIN(U)  *  BRET (  J )  *  Bi<J>*  RELIN(J) 

BLANK  *  0*0 

SUM  C I  )  a  SUMA  ♦  REL I N ( J ) 

1 8  SUMA  *  SUM  < I  ) 

20  SUMA  *  0*0 
STORE  «  0*0 
SUMB  »  0*0 
DO  35  I*  I.  LTERM 
DELWV(i)  »  STORE  +  OELW 
STORE  *  DELWV(I) 

DELW V ( I )  «  DELWV(I)  ♦  WAVE ( 1 > 

DO  34  J*1*N 

DELNU ( 1 )  »  ( 1 •O/DELWVI I )  )  -  FREQ(U» 

I F  <  DELNU (  I >  -  DSL  IT  >  27,  27,  26 

27  SOURC ( J  >  =  ATANF (  (UELNUC1)  ♦  DSL  J  T  )  /  BRET ( J )  )  ♦ 

X  ATANF (  ( DSL  I T  -  DELNU ( 1 )  >/  BRET(J)  ] 

60  TO  29 

28  SOURC <U)  *  ATANF (  ( DELNU ( I  )  +  DSL  I T ) X  BRET(J)  )  - 

X  ATANF  (  ( DELNU  (  1  1  -  DSL  IT)/  BREKU)  ) 

29  SOURC (J)  *  8RET(J)*  BI(U)*  SOURCCJ) 

ADO ( I )  *  SUMB  +  SOURC ( J ) 

34  SUMB  «  ADD < 1 ) 

35  SUMB  a  0*0 _ _ _  _ _ 


PUNCH  95 
AMAX  *  SUM( 1 ) 

BIGD  «  OELAM< 1 ) 

00  39  1  ■£  *  LL 

IFi  AMAX  -  SUMU  1  )  35*39*39 
36  AMAX  *  SUM* 1 > 

BIGD  ■  OfcUAMi 1 ) 

39  CONTINUE 

bMAx  «  AOO ( I ) 

51 GL  *  DELWV< 1 > 

00  43  1  ■  2*  LTERM 
IF  <  5MAX  -  AOO(I>  ) 43 *  43  *  43 
4£  bMAX  «  AOO  <  1 > 

5104.  ■  DELWV4  I  ) 

43  CONTINUE 

IF  4  AMAX  -  UMAX  )  45*46*46 

45  PHI  20  *  BMAX 

WVZO  *  b  t  GL  -  WAVE  I  l  ) 

GO  TO  47 

46  PH 120  *  AMAX 

WVZO  *  51G0  -  WAVE  < 1 ) 

47  PUNCH  99*  BLANK t  BLANK*  WVZO*  PHlZO 
PUNCH  95 

OO  49  1  *1*LL 

SUM  111  »  SUM 4  i  )  /  PHlZO 

OELAM (11  *  OfcLAM 4 i  )  -  WAVE ( 1 > 

49  PUNCH  99*  BLANK*  BLANK •  DEL AM ( 1 ) *  SUM(i) 
PUNCH  95 

DO  51  1  *1  *  LTERM 

AODdi  *  AOO(  1  1/PHIZ0 

DELWV4 1 1  *  DELWV(I)  -  WAVE ( 1 ) 

51  PUNCH  99*  BLANK*  BLANK*  OELWV(I).  ADOU) 
PUNCH  95 
GO  TO  I 

99  FORMAT  <  3F9*3«  E20#tt  ) 

96  FORMAT  <  15) 

95  FORMAT  (/) 

END 


APPENDIX  G 


Spectra!  Line  Shift  Calculations 
(a)  Ho!tsmark‘s  theory 

The  computer  program  used  to  calculate  the  linear  Stark  shift  due 
to  a  given  plasms  fon  density  Is  Hsted  below.  The  following  quantities 
used  in  the  program  are  related  to  terms  In  eqs  (2-21),  (2-28)#  and 
(2-29)  as: 

DION  Ion  density,  Nj 

CLIN  linear  Stark  coefficient  C} 

HAVE  spectral  tine  wavelength,  \Q 

FZERO  average  intermolecular  field,  F 

°  / 

DELTA  Stark  component  shift,  (Av)  t  A' 

/ 

DION  Is  In  units  of  10  cm”"3 ,  CLIN  In  cm  /esu  ,  andJWAVE  In  angstrom 

units.  The  program  output  FZERO  is  In  kv/cm  and  OELTa  Is  given  In 
angstrom  units.  | 

Values  of  the  net  line  shift  ,  for  the  Cul  4531  line  computed 

fof  a  range  of  Ion  densities  are  shown  In  Table  6-1.  The  Intermolecular 
field  strength  Fq  ,  and  the  Individual  p-  and  n-  component  shifts  are 
also  Included.  The  variation  of  the  Cul  4531  net  line  shift  with  Ion 
density  for  a  linear  Stark  effect  Is  Illustrated  in  Fig.  27. 

The  computer  program  used  to  calculate  the  quadratic  Stark  shift 
from  Hoitsmark's  theory  Is  given  below.  The  quantity  CQUAD  Is  the  quad¬ 
ratic  Stark  coefficient  appearing  as  Kj  in  eq  (2-34).  All  other  terms 

*The  electrostatic  unit  of  field  strength  (esu)  is  statvoit/cm. 


are  defined  as  In  tha  above  linear  Stark  shift  calculation.  CQUAD  is  In 
units  of  cm  V(esu)^.  Evaluations  of  the  net  line  shift  for  the  Cul  4531 
line  are  given  In  Table  6-11.  The  variation  of  the  net  shift  with  Ion 
density  for  a  quadratic  Stark  effect  Is  also  illustrated  in  Fig.  27 . 

(b)  Folded  distribution 

Spectral  line  shifts  resulting  from  the  combined  effect  of  nearby 
plasma  Ions  and  electrons  are  obtained  from  an  evaluation  of  the  folded 
distribution,  eq  (2-30),  where  the  electron-broadened  line  profile  is 
given  by  eq  (2-31).  Briefly  what  is  done  is  to  determine  the  electron- 
ion  broadened  line  profile  for  a  given  charged  particle  density.  The 
wavelength  position  at  which  the  resultant  profile  has  a  maximum  value 
yields  the  spectral  line  shift  corresponding  to  the  given  charged  particle 
density. 

In  terms  of  the  dimensionless  ratio  3  ,  the  electron-ion  profile 
for  a  Stark  component  of  a  spectral  transition  that  exhibits  a  linear 
Stark  effect  is  given  by 

m 

J(v)  -  const  x  f  - W(P)...<?§ - -  (G-l) 

Jc  (•  -  cJFoe)2  ♦  Yj 

where 

P  - 

<i>  ■  V  -  v 

WO)  Is  the  Holtsmark  distribution  function,  vq  the  nominal  wavenumber 
of  the  spectral  line,  Cj  a  linear  Stark  coefficient,  Fq  the  average 
Intermolccular  field  given  by  eq  (2-23) ,  and  Yj  the  halfwidth  of  the 
electron-broadened  component  profile. 


G 


Let  W(0 )  be  represented  by  a  linear  function  of  0  in  the 
Interval  *  3  s  ,  l.e, 

W(0)  -  «K  +  bK0  (G-2) 

The  values  of  the  coefficients  a^  ,  fcy  are  calculated  for  each  Interval 
In  the  following  simple  manner.  At  the  point  one  has 

W<%!  '  *K  ♦  W 


and  similarly,  at  the  point 


“«W  *  *K  +  W+l 

where  W(f^),  are  the  values  of  the  Holtsmark  distribution  func¬ 
tion  at  0^,  respectively.  Solving  these  for  the  coefficients 


®K'  bK  on®  ba#: 


*K  * 


-  \  W<W 

*Wl  ‘  *k 
w(^t)  -  W(^) 

K ’  \+\~\ 


U-3) 


(G-J*) 


Values  for  and  W(t^)  are  obtained  from  a  numerical  tabulation  of 

the  Holtsmark  distribution  function,  (Greensteln,  IJfJO)  . 

With  this  linear  approximation  for  W{0)  ,  eq  (G- 1 }  can  be  written 


as 


J(v)  »  E 

K  'J 


(a„  +  b„0)d0 


*V  (a’"CjfoP)2+YJ2 


(G-5) 


where  the  summation  is  performed  for  ail  tabulated  values  of  0  .  The 
constant  in  eq  (G-l)  has  been  set  to  unity  since  the  resultant  profile 


1 5  to  be  norma  I i zed  to  Its  max Imum  value. 


C4 


Define  th«  quant  1  tie*  as 

^<+1 

**  *K  £  (•  -  C.F 


M. 


J'o**V 


I 

L  .  b  r  _ ^ _ 

K  K\  <®'cjFop)2+Yj! 


Thus 


J(v)  »  E  (1^  -f  lk) 
K 

Upon  Integration,  one  has 


(G-6) 


*K  " 


r  k  .  -l  x  \ 

Vc'c'v  tan  77 
lcJFoYj  Tj- 


^*•1 


lk  “ 


{; 


- *  An(X2  +  Y  2)  +  - tan  • 

»«/<,>*  j  <CJFo>S  YJ 


bK” 


-1  X 


:} 


^+1 


where 


Thus 


X  -  CjFo0  -  a) 


y  c/o>'kk  *  lk> 


{<Vjfo  *  v)t,n'1 


bKY,  «  ~  >  ^K+l 


«v 


in(X2  ♦  Yj2)} 


(G-7) 


The  resultant  spectral  line  profile  is  found  from  eqs  (G ~6)  and  (G-7). 


Numerical  evaluations  of  these  equations  leading  to  the  determina¬ 
tion  of  the  line  shift  for  a  given  electron-ion  density  were  performed 
with  the  computer  program  listed  below.  The  input  quantities  are  related 
to  terms  In  the  above  equations  as  follows. 


00) 

dimensionless  field  strength,  ^ 

W(l) 

Hoitsmark  distribution  function,  W(f^) 

OION 

ion  dens ity,  N. 

CUN 

linear  Stark  coefficient,  Cj 

GAMMA 

electron-broadened  haifwidth,  y^ 

WAVE 

spectral  line  wavelength,  kQ 

OELAM 

wavelength  increment 

The  prof  I le 

calculations  were  initiated  at  an  arbitrary  wavelength 

18  •  3 

from  the  nominal  spectral  line  center,  k^  .  OiON  is  in  units  of  10  cm  , 
-1  -1 

CLIN  in  cm  /esu,  GAMMA  Jn  cm  ,  and  WAVE  in  angstrom  units. 

The  average  intermolecular  electric  field  FZERO,  is  determined  in 
accord  with  eq  (2-23)  and  appears  as  an  output  in  units  of  kv/cm.  The 
variable  OEL(J)  represents  the  wavelength  with  respect  to  the  nominal 
line  center,  k  -  kQ  ,  at  which  the  profile  is  to  be  calculated.  DNU(J) 
is  the  related  wavenumber  difference,  u>  .  At  a  particular  wavelength 
DEL(J),  the  coefficients  in  the  linear  approximation  equation  for  the 
Hoitsmark  distribution  function  are  computed  In  accord  with  aqs  (G-3)  and 
(G-4) .  is  represented  in  the  program  b/  the  term  A  ,  and  b^  by 

the  term  6B. 


Estimates  for  y 


J 


are  discussed  below  in  part  (c)  of  this  appendix. 


Gf 


+  given  by  eq  (G-7)  are  then 
obtained.  The  program  first  calculates  the  sum  of  the  Individual  con¬ 
tributions  to  the  quantity  at  the  upper  limit,  \*t  •  X,>IUS 

represents  the  value  of  the  variable  X  at  the  upper  limit  of  each  Indi¬ 
vidual  integration,  f3^+}  .  The  resulting  value  of  this  summation  is 
ADD(J).  Similarly  the  sum  of  the  individual  contributions  to 
at  the  lower  limit  ^  is  obtained  next.  XMINS  is  the  value  of  the 
variable  X  at  the  lower  limit  of  each  Integration,  .  The  value  of 
this  summation  Is  called  SUN(J).  The  quantity  ADD(J)  *  ADD(J)  -  SUM(J) 
then  determines  J(v)  in  accord  with  eq  (G-6)  at  the  wavelength  OEL(J) . 
This  procedure  is  repeated  at  each  wavelength  step  to  describe  the  line 
profile  ADO(J)  St  the  wavelength  DEL(J)  for  each  Stark  component. 

The  remainder  of  the  program  is  directed  towards  the  determination 
of  the  wavelength  at  which  the  line  profile  has  a  peak  value.  The  quan¬ 
tities  AOO(J)  are  first  normalized  to  their  maximum  value  BMAX,  and  pre¬ 
sented  as  output  along  with  the  associated  wavelengths  OEL(J).  In  order 
to  Improve  the  determination  of  the  peak  value  of  the  line  profile  a 
second-order  least  squares  fit  in  the  vicinity  of  the  above  calculated 
maximum  value  of  the  computed  profile  is  then  made.  Two  points  on  either 
side  of  this  computed  maximum  are  used  in  the  fitting  calculation.  The 
subroutine  POLYFT,  previously  listed  In  Appendix  C  is  used  to  obtain  the 
least  squares  fit.  The  wavelength  at  which  the  fitted  profile  is  a  maxi¬ 
mum  is  called  XMAX.  It  represents  the  shift  In  a  Stark  split  component 
of  a  spectral  line  under  the  combined  Influence  of  electrons  and  ions  In 
the  vicinity  of  an  emitting  atom. 

Values  of  tha  net  line  shift  ,  for  the  Cu!  4531  line  obtained 


Values  of  the  quantities 


Q7 


from  tbefolded  distribution  calculation  for  several  value*  of  Ion  dens- 
Ity  are  shown  In  Table  G-IM.  Shift*  of  the  Individual  p-  and  n-  compo¬ 
nent*  are  also  Included.  The  variation  of  the  Cul  net  wavelength 

shift  with  Ion  density  for  a  linear  Stark  effect  Is  Illustrated  In  Fig. 

27. 


(c)  Electron-broadened  halfwidth  calculations 

Kolb  0956)  has  derived  the  following  expression  for  the  electron- 
broadened  line  halfwldth  ,  for  a  spectral  line  that  exhibits  a  linear 
Stark  effect. 


I6N 


3v 


(eC j) ^(2.36  -  21n  fij) 


(G-8) 


where  is  the  electron  density,  Cj  is  a  I  inear  Stark  coefficient. 


e  the  electronic  charge  and 
v  -  (3kT/m),/2 


flj  *  2eCj/pv 
p  -  (kT/Brre2^)  1/2 

m  represents  the  mass  of  an  electron  and  T  the  temperature. 

For  the  p-  and  n-  Stark  components  of  the  Cul  **531  line  one  has 
numer leal ly 

Y  -  1.1+24(1  -  0,84?  Lnb  )n  /(T),/2 
P  P  e 

where 

5p  -  72.7 i'^)  ,/2/T 


Y  -  0.801(1  -  0.847  ln6  )n  /(T)l/2 
1 1  n  e 


and 


where 


6n  -  54.0{n  }3/2/T 
n  6 

i  g  * 

n  Is  the  electron  density  In  units  of  10  cm  ,  and  Y„»Yn  »r®  the 
e  p  n 

component  halfwtdths  In  wavenumber  units,  cm  3. 

A  program  for  computing  these  e lectron-broadened  line  halfwidths 
for  a  given  electron  density  and  temperature  is  given  below.  The  program 
input  quantities  DION  and  T  represent  the  electron  density  value  ng 
and  temperature  T  respect iveiy.  The  halfwidth  value  yp  is  called 
GAMMA,  and  y  is  called  GAMMB.  Table  G-IV  shows  the  calculated  Y-valuas 
for  several  choices  of  electron  density  and  temperature. 

These  values  for  Yj  were  employed  to  carry  out  the  folded  dis¬ 
tribution  line  shift  calculations.  Specifically  the  electron-broadened 
halfwldth  at  a  temperature  of  5000°K  for  a  given  charged  particle  density 
was  used.  To  ascertain  that  the  calculated  line  shifts  were  relatively 
Insensitive  to  the  value  of  temperature  chosen  over  a  reasonable  tempera¬ 
ture  range,  comparative  shifts  were  calculated  for  a  fixed  particle  dens¬ 
ity  and  different  temperatures.  The  numerical  results  given  in  Table  G-V 
show  that  for  the  temperature  values  considered  the  variations  in  line 


shift  are  Insignificant. 


Table  G-l 


Cut  453 1 

Linear 

Stark  Shift 

(Holtsmark  Theory) 

loo 

F 

Component  Shift 

5 

Dens  1 ty 

o 

P 

n 

s 

(cm'3} 

(kv/cm) 

(A) 

(A) 

(A) 

5.0  1016 

51 

0.07 

0.05 

0.06 

1.0  IO17 

81 

o.n 

0.08 

0.  10 

5.0  1017 

227 

0.32 

0.24 

0.28 

o 

© 

OD 

376 

0.51 

0.38 

0.45 

5.0  IO18 

1100 

1.50 

1.12 

1.31 

1.0  io’9 

1745 

2.38 

1.78 

2.08 

5.0  IO19 

5103 

6.96 

5.22 

6.09 

Table  G-t I 

Cul  4531  Quadratic  Stark  Shift  (ttoltsmark  Theory) 


Ion 

F 

Component  Shift 

6 

Dens  1 1 y 

O 

P 

n 

5 

(cm"3) 

(kv/cm) 

(A) 

(A) 

(A) 

1.0  1017 

.  81 

0.30 

0.23 

0.27 

5.0  !017 

111 

2.61 

1.96 

2.28 

1.0  10,B 

376 

6.57 

4.93 

5.75 

5.0  1018 

1100 

56.17 

42.  13 

49.15 

Table  G-lll 


Cut  453!  Linear  Stark  5hlft  (Folded  Distribution) 


ion 

Component 

Shift 

6 

Density 

P 

n 

s 

(cm*3) 

(A) 

(A) 

(A) 

i.O  1018 

0.58 

0.44 

0.51 

5.0  lo18 

1.64 

1.23 

1 .44 

f.O  1019 

2.61 

1.95 

2.28 

E lectron-Ha l fwidth  Values  for  Cui  4531  Stark  Components 


E  lectron 

Temperature 

Density 

Yp 

(cm-3) 

(°K) 

(cm'  * 

o  io’9 

1000 

0. 145 

5000 

0.092 

10000 

0.074 

0  io18 

1000 

0.572 

5000 

0.393 

J0000 

0.320 

0  io’9 

1000 

1.01  1 

5000 

0.727 

10000 

0.597 

Halfwidth 

(cm"') 

0.088 

0.055 

0.043 

0.353 
0.235 
0. 190 

0.632 

0.437 

0.356 


Table  G-V 


a. 

Variation  of  Line  Shift  with  Temperature 


Temperature 

T(°K) 

1000 

5000 

10000 


Line  Shi  ft 
fis.(A) 

1  447 
1.439 
1.436 


18  -3 

An  electron  density  of  5*10  cm 


was  assumed 
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C  ELECTRON  bROAOENING  HAL.Fl1?  £U  fH  EETJMmTeE 
I  HEAD  99*  DION 
RE AO  99 ♦  T 
A  *  SQR TF  <  0 1  ON ) 
b»  1 .G/80RTF ( T ) 

DELA  -  72 »7* A 
DELA  =  DELA/T 

GAMMA*  1  .424*  <  1  *0-(J»Q47*L0GF(0ELA>  } 
GAMMA* <  GAMMA*U I  ON ) *b 
DELti  *  54 • G*A 
OELti*OELb/T 

GAMMb  *  0.b01*( 1 *0-U.b47*L0GF<UELb> } 

GAMMti*  (  GAMMti*U  iON)*a 

PUNCH  99*  UlON*  T«  GAMMm  ,  (jAMMt) 

GO  TO  1  t 

99  FORMAT  <  2F11.4*  2tlb*8  ) 

END 


BLANK  PAGE 


HI 


APPENDIX  H 

Densitometer  Slit  Broadening  Effects 

The  operation  of  the  Baird-Atomic  densitometer  Is  schematically 
Indicated  In  Pig,  H- 1 .  A  light  beam  from  a  source  lamp  Is  chopped  by 
means  of  a  rotating  shield.  The  chopped  beam  is  first  transmitted 
through  the  film  under  study,  and  alternately  through  a  mechanically 
driven  aperture,  i.e.  the  reference  beam  path.  After  passage  through  the 
film  the  resultant  beam  Is  focused  on  the  scanning  slit  and  photomulti¬ 
plier  for  comparison  with  the  reference  beam.  Differences  in  light 
intensity  actuate  a  servo  amplifier  and  motor  which  in  turn  drives  the 
adjustable  aperture  in  the  reference  beam  path.  The  servo  system  is  at  a 
null  balance  when  both  beam  Intensities  are  equal. 

Let  Jq  be  the  source  lamp  intensity  and  Jg  the  light  intensity 
of  the  reference  beam  after  passage  through  the  adjustable  aperture. 

Thus 

Ji  -  J0  I0'D  (H-l) 

where  D  is  the  density  value  determined  by  the  opening  of. the  adjust¬ 
able  aperture.  If  D(x,y)  represents  the  density  corresponding  to  a 
given  position  x,y  ,  on  the  film,  then  the  Intensity  of  the  chopped  beam 
after  passage  through  the  film  is  given  by 

L/2  w/2 

J  -  j  \  J0  IO"D(x'y)  dx  dy  (H-2) 

-L/2" -w/2 

where  W, L  represent  the  scanning  slit  width  and  length  respectively, 

A  full  description  of  this  instrument  is  given  by  A.O.  Carpenter  and 
J.U.  White,  Analytical  Chemistry  2£,  1473  (Oct  1953). 


H2 


projected  onto  the  film  plane.  The  y-dlrectlon  corresponds  to  the 
direction  of  increasing  time  on  the  film  record,  and  the  x-dfrection  is 
related  to  a  wavenumber  scale  from  a  dispersion  relation.  Neglecting 
temporal  variations  on  the  film  record  over  the  iength  of  the  slit  L 
one  has 


where  v  Is  measured  from  v£  .  Thus 

I0-D(v)  „  IO"D(vc)  10-®''  io-*2 


H3 


where 


a  * 


V 

c 


& 


*  dv 


v 


c 

For  small  a  and  p  the  above  expression  becomes 

10‘0<V)  «  10“D^  {l  *  «v  +  ^  »2v2  -  ^ 
to  second  order  terms  In  v  .  Thus  from  eq  (H-5) 


w  10'° 


Performing  the  integration  there  finally  results 

2 

D(vc)  -  D  w  log  1  +  v2j-  (H-6) 

Thus  the  optical  density  as  measured  by  the  densitometer  deviates  from 
the  true  film  density  (at  the  center  of  the  scanning  slit)  by  the  amount 
indicated  on  the  right  hand  Side  cf  eq  (H-6). 

If  i  is  the  film  exposure  associated  with  the  density  D(v£)  and 
I  that  associated  with  D  then 

(ft 

2  l/Y 

*/*m  **  {l  +  (--,“  ?)v2}  (H-7) 

where  Y  is  the  slope  of  the  photometric  calibration  curve.  Expanding 
the  bracketed  expression  one  has  finally 


(*  -  *)/•  *  (a2  -  &)  v2/ 1 2Y 

m  m 

The  constants  a  ,  P  are  evaluated  as  follows. 


(H-8) 

From  their  defini- 


t  ion 


«  *  <1?)^  “  Y/#^v 


and  stmi  Ur  I y 


20 


<4 

dv^ 


v 


c 


where  I  is  given  by  eq  (2-11)  in  the  text.  After  some  manipulation 
one  finally  arrives  at 


(H-9) 


where 


C.  •  ten’*  X  -  tan’1  Z 


Cf>  “ 


1  +  X2  1  +  Z2 


C  .  Z  _  X 
3  (i  +  Z2)^  (I  ♦  X2)2 

X  -  (v  -v  ♦  A) /A/2 

v  O 

Z  -  (V  -v  - A) /A/2 

v  O 

and  y  •  1.0  has  been  assumed. 

A  listing  of  the  Fortran  program  used  to  numerically  evaluate 
errors  that  result  from  the  use  of  a  finite  width  scanning  slit  is  given 
below.  Several  important  quantities  named  in  the  program  are  related  to 
variables  in  the  above  equations  in  the  following  manner 
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WAVEK 

Wavelength 

of  spectral  line,  X 

BREAT 

Resonance 

halfwidth,  A/2 

SLIT 

Spectrograph  entrance  slit,  A/2 

ERROR 

Pens  1  tome ter  slit  width,  X  «  1/v 

OWAVE 

1 ncrement 

in  wavelength 

These  input  quantities  are  all  in  angstrom  units 

The  output  of  the  program  consists  of: 

DEL  wavelength  difference,  Xc  -  XQ 

ERR  relative  exposure  error,  (#  -  #  )/# 

m  m 

A  sample  calculation  of  the  error  due  to  a  finite  width  densitom¬ 
eter  silt  has  been  made  for  a  Cul  4651  profile  of  narrow  width.  At  a 
time  of  7.7  p.sec  the  resonance  halfwidth  A  ,  for  this  line  is  3.IQA; 
the  spectrograph  entrance  si  it  is  5.50A  at  this  wavelength.  The  pro¬ 
jected  width  of  the  scanning  slit  onto  the  film  plane  is  about  0.0)3  mm 
and  corresponds  to  a  wavelength  range  of  approximately  0.5A  for  the  nomi¬ 
nal  40A/mm  dispersion  at  this  wavelength.  The  calculations  indicate  that 
the  maximum  relative  exposure  error  FRR  is  less  than  0.5  percent  through¬ 
out  the  line  profile. 
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FIG.H-I.  SCHEMATIC  DIAGRAM  OF  DENSITOMETER 


c 


DENS  1  TOME  Tfc,«  SL  ST  CORRECTION  MX  2 

1  READ  99*  WAVE*  »dREAT  »  SC  1  T  4  ERROR  tDWAVE  «N 
PUNCH  99*  WAVEK  »tiREAT  «SC  l  TtERROH 
RATIO  *  ERROR/dREAT 

PUNCH99.RAT  lO 

WAVE*  WAVE K  +  dRE AT 

HFW 1 T  »  <  1  »0/WaVex ) “ <  1  • U/WMVt  ) 

WAVE*WAVEK.  +  bL  1  T 

SC  1 T*  <  1  .O/waVEK | - l  1  .O/WAVE ) 

l  =  1 

AAA=0 «U 

2  DEL=DWAVE*AAA 
WAVE  =WAVEK+UEL 

ya(  1 #0/WAVEK  J - <  1 .O/WAVt ) 

X*<  Y+SC 1 T  > /ME w i T 
Z*  <  Y-SL  J  T )/HFWl T 
1F12)3.3*4 

3  ZZ<*  -Z 

Cl «ATANF ( X  »  +aTaNF  <  ZZ  > 

CO  TO  5 

4  C l *ATANF < X > -aTanf { Z > 

5  CZ*{ I.U/I l.U+X*X) 1  -  l i .O/ ( 1 «Q+Z*Z ) > 

6  C3=(ZX(  1  .0  +  *.#Z)  MM  l  *0/<  1  *U+Z*Z)  > 

7  C3=C3“<XZ4 l *U+X*X> )#( l.U/l UU+X*X) I 

dETA  =  (C2*C£/< Cl *Cl > J  -  C3XC1 

BETA*BETA/(dREAT*dREAT ) 

ERR«  O,O833*EETA*ERR0W*ERW0R 
PUNCH9Q «  DEE  *  ERR 
IF«1-N)8t  l«i 
d  1*1+1 

aAA= AAm+ 1 • U 
CO  TO  d 

99  FORMAT  <  bF9.3«  lb  ) 
ye  FORMAT  <  F9.3«  2E2U#d  ) 

END 
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Spectra]  Line  Se 1 f -Absorpt i on  Calculations 


Numerical  estimates  of  the  effects  of  se I f-absorpt i on  on  the 
observed  spectral  lines  are  derived  from  a  one-d i mens i one  1  absorption 
model.  Distortions  of  the  line  profile  upon  passage  through  a  thin 
plane-paral lei  absorbing  region  are  given  by  egs  (2-36).  (2-38),  and 
(2-39)*  The  optical  depth  of  a  thin  slab  t  (v)  is,  in  cgs  units 


Tj(v) 


-8,448 


-3 

10 


abs 


dr 


(1-0 


where  N  ,  the  number  density  of  absorbing  atoms  and  y  ,  the  colli- 

3 

sion  damping  constant  depend  upon  the  radial  position.  The  number 
density  of  absorbing  atoms,  i.e.,  of  excited  atoms  in  the  lower  energy 
level  of  the  affected  spectral  transition,  is  related  to  the  total 
number  density  of  neutral  atoms  as 


N  JH  -  (g/U)  exp  (-E/kT) 

3 


0-2) 


assuming  local  themodynamic  equilibrium  in  each  element  of  the  slab, 
g  is  the  statistical  weight  of  the  lower  energy  level,  U  the 
partition  function  for  the  neutral  atom,  E  the  energy  of  the  lower 
level,  and  T  the  local  temperature  in  the  slab. 

The  neutral  atom  density  and  temperature  distributions  deter¬ 
mined  from  the  fluid-dynamic  model  of  a  jet  collision  are  used  in 
the  calculations  (Lukasik  et  al,  1984).  If  Nc  represents  the 


12 


neutral  atom  density  and  Tc  the  temperature  at  the  Impact  center, 
rf  the  position  of  the  front  of  the  expanding  collision-produced 
material  sheet,  and  r  a  redial  position  in  the  sheet,  then  near  the 
expanding  front  one  has 

N/Ne  -  26.0  (1  -  r/rf)4  (1-3) 

The  temperature  ratio  T/T  was  expressed  as  a  fourth-order  poly¬ 
nomial  In  r/r^  (see  program  listing).  Nc  and  r^.  are  given  as 
functions  of  the  time  after  initial  Impact  as 


N  (cm“a)  »  35-0  1020/ta 

V 

(I-1 

ff  (cm)  »  1 . 72  t 

(1-5) 

for  t  In  usee  .  T  was  treated  as  a  parameter  in  the  calculations. 
Combining  the  above  equations  one  finally  arrives  at 


With  the  substitutions 


5 

r . 

J 


r 

o 


-  r 

o 

4-  X 


+ 


rj-i  -  r  +  x 
J  o 


where  rQ  is  the  positron  for  which  the  neutral  atom  density  Is 

given  as  N  ,  then 
c 


Tj(v) 


-  («  e/5  rf 4) 


(1-8) 


The  alteration  of  the  shape  of  a  spectral  line  with  an  Incident 
intensity  profile  Jj_j(v)  upon  passage  through  the  slab  is  obtained 
from  eqs  { 2—36)  and  (1-8).  Hence,  the  change  in  a  given  Intensity 
profile  Jq(v) ,  after  passage  through  a  number  of  plane-parallel  ab¬ 
sorbing  slabs  Is  obtained  by  repeated  application  of  the  combined 
eqs.  (2-38)  and  (1-8). 

The  computer  program  used  to  numerically  evaluate  the  profile 
of  a  self-absorbed  spectral  line  is  given  below.  Several  input  terms 

4 

defined  in  the  program  are  related  to  quantities  in  the  above  equations 

as 

D  neutral  atom  density,  N 

o 

TIME  t i me  after  collision,  t 

G  statistical  weight,  g 

U  partition  function,  U 

E  energy  of  lower  level,  E 

TC  temperature  at  Impact  center,  Tc 

WAVE  wavelength,  \ 

o 

F  oscillator  strength,  f  k 


A  collision  damping  constant,  Y 

CONST  absorption  slab  thickness,  r.  -  r.  , 

J  J-l 

<  A 

0  is  in  units  of  10  cm“a ,  TIME  in  usee,  E  in  ev,  TC  in  °K,  WAVE 
and  A  are  in  angstrom  units. 

The  calculations  proceed  in  the  following  manner ,  The  initial 
position  ft  at  which  the  neutral  atom  density  is  0  is  first  cal¬ 
culated  in  accord  with  eg.  (1-3).  The  intensity  profile  of  the  radia¬ 
tion  incident  on  the  first  absorbing  slab  RELIN(J)  is  then  determined 
at  fixed  wavenumber  increments  DNU(J) .  This  profile  has  a  Lorentzian 
shape  with  a  halfwidth  determined  by  the  neutral  atom  density.  DD(J) 
is  the  wavelength  equivalent  of  DNU(J).  The  radial  positions  of  the 
faces  of  the  absorbing  slab  are  called  X  ,  and  XH  ,  and  are  equivalent 
to  x_  and  x+  in  the  above  equations.  The  average  collision  damp¬ 
ing  constant  A!  and  temperature  TEMP  in  the  absorbing  slab  are  then 
calculated.  The  quantities  ALPHA  and  B(J)  which  are  the  equivalent 
of  cr  and  P  in  the  above  equations  are  obtained  and  the  opacity  of 
the  slab  BRIT(J)  at  the  wavelength  increment  00(J)  is  then  calculated 
as  in  eq.  (1-7).  The  intensity  profile  after  passage  through  the 
first  slab  is  found  and  re  lab  led  BRIT(J).  This  profile  becomes  the 
incident  profile  for  the  next  slab,  RELIN(J).  The  radial  positions 
X,  XM  of  the  next  adjacent  absorbing  slab  are  then  obtained  and  simi¬ 
lar  calculations  leading  to  the  intensity  profile  after  passage  through 
the  second  slab  are  done.  These  calculations  are  further  repeated 
until  the  position  of  the  front  of  the  expanding  collision-produced 
sheet  is  reached. 


The  program  output  consist*  of  the  initial  intensity  profile 
before  absorption,  the  total  opacity  OPAC  of  the  absorbing  regions  as 
given  by  eq.  (2-36),  and  the  intensity  profile  of  the  emerging  radia¬ 
tion.  This  latter  profile  Is  also  normalized  to  its  peak  computed 
value  to  facilitate  comparison  with  profiles  of  observed  spectral  tines. 

A  number  of  self-absorbed  spectral  profiles  for  the  Cut  4651 
and  5153  transitions  have  been  computed  the  results  of  which  ere  shown 
In  the  following  figures.  Values  of  the  neutral  atom  density  used  in 
the  calculations  were  obtained  from  the  measured  density  variations 
illustrated  in  Fig.  26.  Estimates  of  the  Initial  collision  damping 
haifwidth  were  obtained  from  these  density  values  and  known  oscillator 
strengths  be  means  of  eq.  (2-2!a).  Tha  temperature  of  the  material 
at  the  entrance  face  of  the  absorbing  layers  mis  treated  as  a  parameter 
in  the  computations  and  is  indicated  in  the  figures. 
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I  READ 
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D 
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TIME 

READ 
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G.U.E 

‘ 

READ 
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T 

READ 
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WAVE  «F  .  A 

i 

READ 
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DELaM.LTERM 

5 

1 

REAO 
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CONST 

4 

$ 
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PUNCH  99*  WAVE  »  F  *  A 

PUNCH  99*  CONST  «  G  *  U.  t 

PUNCH  96 

TT»TIME*T1ME 

TTR=1  •O/'TT 

u 1  ST  =D*TT* 1 .0996E-Ob 
OiSTsSORTF  <DicsT  ) 
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RF=1 .O/RF 
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FIG.I-4.  SELF-ABSORBED  Cu  I  5153  PROFILE 
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